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1 . 1 Differential cross section angular distributions for the deuteron photo-disintegration 
process for selected incident photon energies below 0.8 GeV. Data points 

are from the following experiments: Bonn [5] (open squares), Frascati [S] 
(open crosses), MIT [TU] (open triangles), LEGS T2] (open dots), and 
Mainz/DAPHNE fT^ (solid dots). The curves belong to the following the- 
oretical models: the Impulse Approximation (dashed line), the Coupled 
Channel Approach [S] (up to 0.4 GeV) (continuous line), the Diagrammatic 
Approach [B] (dot dashed line), and the Modified Current Conservation [7] 
(dotted line). As can be seen from plots, the predictions of the theoretical 
models have problems in reproducing the data, starting from relatively low 
incident photon energies. (The plots are a courtesy of P. Pedroni - INFN, 
Pavia) i 

1.2 Legendre coefficients for the differential cross section. At the highest pho- 
ton energy reached at Mainz (around 0.8 GeV), the Ai and A3 coefficients 
exhibit a small but significant deviation from the trend for the lowest en- 
ergies. Data points are from the following experiments: Bonn [8] (open 
squares), Frascati [9] (open crosses), MIT [10] (open triangles), LEGS [12] 
(open dots), and Mainz/DAPHNE [13] (solid dots). The curves belong to 
the following theoretical models: the Impulse Approximation (dashed 
line), the Coupled Channel Approach [5] (up to 0.4 GeV) (continuous line), 
the Diagrammatic Approach [B] (dot dashed line), and the Modified Current 
Conservation [7] (dotted line). (Courtesy of P. Pedroni - INFN, Pavia) . . 

1.3 Total photo-disintegration cross section. In the energy region around 0.8 GeV 
the total cross section behavior is not a smooth decreasing function of the 
incident photon energy: this could be explained by the intermediate exci- 
tation of baryon resonances in the second resonance region which includes 
Pii (1440), 1)13(1520), and S'ii(1535). Data points are from the following 
experiments: Bonn [5] (open squares), Frascati [5] (open crosses), MIT [TU] 
(open triangles), LEGS [l^ (open dots), and Mainz/DAPHNE [13] (solid 
dots). The curves belong to the following theoretical models: the Impulse 
Approximation [4] (dashed line), the Coupled Channel Approach [5] (up to 
0.4 GeV) (continuous line), the Diagrammatic Approach [6] (dot dashed 
line), and the Modified Current Conservation 7 (dotted line). (Courtesy of 

P. Pedroni - INFN, Pavia) i 

1.4 Results form the NE8 experiment: the cross section shows scaling for inci- 
dent photon energies larger then 1 GeV. The solid line represents a meson- 
nucleon description which do not account for the data starting from 0.5 GeV 

of incident photon energy 
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1.5 Deuteron photo-disintegration differential cross section multiplied by s^^. 
Experimental data are from Mainz [13] (open/magenta circles), SLAC [17l 
fT8] (open/green squares), Jlab Hall C [24l[25] (open/blue triangles), and Hall 
A [26] (open/black diamonds). Points are multiplied by s^^, and arrows 
indicate the expected threshold for the onset of the CCR scaling, where 
the momentum transfer to the proton, t is of the order of 1 GeV^ [22] . 
Theoretical curves are described in the text 

1.6 Angular distributions for three incident photon energies: 1.6 (top plot), 
1.9 (middle plot), and 2.4 GeV (bottom plot). The experimental data are 
from Refs [TH] [24l [25l [26] while the curves represent the theoretical predic- 
tions of the QGSM [31] and the HRM [35 . Prediction from other models are 
not reported since the were not provided or could not be calculated. From 
the plots it is seen that the QGSM well describes the data. This model also 
predicts a forward to backward asymmetry. More data especially at very 
forward and backward angles are necessary to verify this prediction 

1.7 Proton polarization transfers Cx (transverse in-plane, left panel) and Cz 
(longitudinal in-plane, right panel) in the deuteron photo-disintegration at 

= 90°. The y-axis is perpendicular to the scattering plane, the z-axis is 
parallel to the scattered proton momentum and the x-axis is perpendicular 
to y and z. Experimental data are from |41] and have been corrected for 
spin rotation due to the LAB — ^ CM transformation. The solid/red curve 
represents the QGSM prediction for Cz [13 ■ The QGSM calculation for Cx 
is underway [l 

1.8 Asymmetry S for linearly polarized photons in the deuteron photo-disintegration 
at = 90° as a function of the incident photon energy. Experimental data 

are from Ref. [47]. The solid/red curve is the prediction of the QGSM [42]. [T 

1.9 Left panel: planar diagram for the reaction ah — > cd. Right panel: the same 

for the reaction ab ^ X. Full lines denotes quarks, wavy lines, gluons. . . [T 

1.10 Left panel: Cylinder type diagram. Right panel: cutting of these diagrams 
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in the s-channel 17 



1.11 Topological diagram with one handle corresponding to the exchange of two 
Ponierons in the t-channel [17 

1.12 Hadron as a "bubble" in the vacuum [T 

1.13 Interaction and decay of a color tube [T 

1.14 Simple linear baryon Regge trajectories separated by A J = 2 [70] [2 

1.15 Color strings: (a) massless string, (b) homogeneous string connecting a qq 
couple, (c) inhomogeneous string with variable tension due to color vacuum 



polarization [21 

1.16 Total cross sections of the reactions pp — > (i7r+ and pd — > p-K^ as a function 

of the CM energy, and predictions of the QGS model [2 

1.17 Planar diagrams describing the binary reactions tt^tt" — >■ tt+tt" (a: left panel), 
NN TT+vr" (b: central panel), and NN NN (c: right panel) [i 

1.18 Diagrams describing three valence quark exchanges in the i-channel (left 



panel) and u-channel (right panel) [27 

2.1 Schematic layout of the JLab Continuous Electron Beam Accelerator. . . . [s^ 



2.2 Left panel: Schematic view of the Hah B enclosure, including the CLAS 
metal supporting structures. The tagging spectrometer is shown at the 
entrance on the Hall, upstream whit respect to the electron beam direction. 
Right panel: The CLAS detector magnified view showing the schematics for 
the detector subsystems 

2.3 Hall B beam-line, the CLAS detector, and associated equipment 

2.4 Hall B photon tagging system. Photons produced in the radiator continue 
toward the CLAS target following the same beam-line through the magnet 
yoke that is traversed by the electron beam when the tagger is not in use. 
The photon-tagging system can tag photons with energies between 20% and 
95% of the incident electron energy, and is capable of operating with electron 
beam energies up to 6.1 GeV 

2.5 Tagger spectrometer electronics scheme. See the text for details 

2.6 HaU B Total Absorption Counter (TAC). The TAC is essentially 100% ef- 
ficient allowing for the efficiency of the diff'erent elements of the tagger ho- 
doscope to be determined. Due to counting pile-up problems, the TAC can 
only be operated at beam currents up to 100 pA, and must be retracted 
from the beam-line under normal running conditions 

2.7 A schematic top view of the CLAS detector cut along the beam line. Typical 
photon, electron and proton tracks (from top to bottom) from an interaction 
in the target are shown (the photon track through the CLAS drift chamber 
system is drawn only to guide the eye) 

2.8 Unpolarized cryogenic target cooling system (from Ref. [97]) 

2.9 The start counter assembly is a set of scintillators surrounding the target 
with the same geometric acceptance of the CLAS returning the time at 
which a charged particle has left the target region 

2.10 Left panel: Contours of constant magnetic field for the main CLAS toroid in 
the mid-plane between two coils (see also Fig. 12. 7p . Right panel: Magnetic 
field vectors for the CLAS toroid transverse to the beam direction in a plane 
centered on the target. The six coils are shown in cross-section 

2.11 Left panel: Portion of a region three drift chamber, showing the layout of its 
two super-layers. The high-lighted drift cells have fired at a charged particle 
passing by. Right panel: Section of the whole region 3 drift chamber. . . . 

2.12 Coverage of one CLAS sector by the time-of-flight scintillators 

2.13 Scheme of the electronic setup for a TOP counter. Each counter consists 
of single scintillator bar with a PMT at each end. The timing resolution of 
the scintillators varies from 80 ps for the shorter counters to 160 ps for the 
longer ones 

2.14 Upper panel: schematic view of one module of the forward electromagnetic 
calorimeter. The expanded view shows the three layers of lead and scintilla- 
tor (conventionally referred as u, v, and w) making up the calorimeter sand- 
wich structure. Lower panel: the large angle calorimeter tapered geometry 
(up- left), its lead-scintillator sandwich structure (up-right), and granularity 
(bottom-center) 

2.15 Event rate and integrated events accumulated by the CLAS data system vs. 
time for a three year period 



2.16 Data flow schematic for CLAS, with the primary stream using the ETl 
shared memory. Also shown are a variety of monitoring processes accessing 
data on other shared memory (ET2, ET3) managed by the ET system. See 
text for additional details 

3.1 Top pictures: Schematic view of the experimental setup used during the g2 
run period. From left to right: electron beam, radiator, the tagger spec- 
trometer, the photon beam-line, the deuteron target (TGT), the CLAS de- 
tector, and the down-stream photon flux normalization components. For 
the present experiment, the hadronic events generated in the target volume 
are detected in CLAS using the time-of-flight system, in coincidence with 
a tagged photon trigger and a start counter signal. Bottom pictures: the 
Level 1 trigger and the data acquisition system. Only the charged trigger 
scheme is indicated since the neutral trigger is not used in the current analysis. 

3.2 Block diagram of the off-line data processing. The raw data files (simulated 
-left panel- or measured -right panel-) are converted into reconstructed files 
and processed by CLAS applications or user written code to obtain standard 
n-tuple files 

3.3 The tracking procedure in the CLAS drift chambers system. Left panel: a 
charged particle passing through a super-layer. The field wires are indicated 
by the open dots while crosses indicate the sensitive wires which collect the 
released charge. The drift times are associated to distances of closest ap- 
proach (DOCAs) represented by the dotted lines. Right panel: two different 
values for the distance of closest approach of the particle from the sense wire 
can be defined, giving origin to a loss of spatial resolution. See text for details. 

3.4 Typical missing mass spectra calculated for the six CLAS sectors from 
Eq. 13.51 for an incident photon energy of E.y =1.05 GeV and a forward 
proton scattering angle 9^^^ = 25°. The peak corresponding to the neu- 
tron rest mass can be clearly identified, and the shape of the distributions 
is well reproduced by a Gaussian plus an exponential fitting function. The 
red/dotted arrows shown in the panels represent the (— Bcr, 3cr) interval defin- 
ing the width of the neutron peak, being a the Gaussian function standard 
deviation 

3.5 Example of missing mass spectra calculated for the six CLAS sectors from 
Eq. 13.51 for an incident photon energy of ==1.05 GeV and a backward 
proton scattering angle 9^^^ = 115°. The distribution is less pronounced 
compared to the one obtained at forward angles since the momentum res- 
olution of the CLAS is worse in the backward direction. Nevertheless the 
neutron signal is clearly evident and the shape of the distributions is well 
reproduced by a Gaussian plus an exponential fitting function. Again, the 
red/dotted arrows represent the (—3(7, 3(t) interval defining the width of the 
neutron peak, being a the Gaussian function standard deviation 

3.6 Behavior of the missing mass cuts as a function of the proton scattering angle 
9p^ for four incident photon energies around 1 GeV. The solid horizontal 
line represents the value of the neutron rest mass, while the different dotted 
lines represents the missing mass cuts for the six CLAS sectors 
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3.7 Behavior of the missing mass cuts as a function of the proton scattering angle 
9p^ for four incident photon energies around 2.95 GeV. The solid horizontal 
line represents the value of the neutron rest mass, while the different dotted 
lines represents the missing mass cuts for the six CLAS sectors. In this 
case, the cuts are wider than those calculated for lower photon energies 
nevertheless the missing mass distributions are still peaked enough for the 



3.8 Distribution of the vertex z-component for proton events originating in the 
target volume. The distribution shape reflects the target dimensions (10 
cm length). The small peak on the left is the contribution of a target 
support structure and is eliminated by a cut on the vertex z-coordinate. 
The information on the interaction vertex and the x and y positions of the 
beam is contained in the MVRT bank in the case of multi-track events. The 
solid line represents a fit to the distribution defining the effective target 



3.9 TDC timing spectrum of a single T-counter. The peak which caused the 
trigger and the accidental background can be clearly identified. Utilizing 
the accidental hits in the tagger, the normalization method determines the 
number of tagged photons from each T-counter. This is done by measuring 
the number of accidental hits of each T-counter in a given time interval t, 



3.10 Schematic representation of a normalization run: the electron beam current 
is lowered to 100 pA and the total absorption counter is inserted into the 
beam-line, downstream of the CLAS. The trigger is formed only by the 
tagger master OR (MOR). This setup allows for a direct measurement of 
the photon flux from the total absorption counter (composed by a set of four 



3.11 Behavior of the tagging efficiency for the tagger spectrometer T-counters for 
different runs (indicated by the different symbols). The tagging efficiency is 
defined as the number of tagged photons per energy interval divided by the 
total number of counts observed in the tagging counter defining that interval. 
The knowledge of this quantity is very important because it directly affects 
the determination of the photo-reaction cross sections and strongly depends 



3.12 Variation of the tagging efficiency with respect to the mean value for some 
tagger spectrometer T-counters. As can be seen the variation is around 



fitting procedure to be applied |67 



length [68 



which is called the out-of-time window 68 



lead-glass blocks with a detection efficiency close to 100%) |71 



on the beam setup |72 
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3.13 Distribution of the z-component of the interaction vertex for the proton 
events (normalized to the number of incident photons) between a full (blue 
line) and empty target (black line) run. The difference between the two dis- 
tributions (red line) is seen to be less than 1% in the interval (— 15 < z < —10) cm, 
confirming the validity of the photon normalization procedure [14 
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3.14 Behavior of the normahzed number of protons as a function of time for 
the "August" run period characterized by the tagger spectrometer in the 
prescale ON configuration. Top plot: for each run several files are analyzed 
(represented by the single points) and the normalized proton yield is seen 
to be constant in time. Bottom Plot: Overall distribution of the normalized 
proton yield fitted with a Gaussian function. The ratio between the standard 
deviation a and the mean /i is less than 3% 

3.15 Distribution of the A values calculated for the normalized photo-disintegration 
proton yield N\. at an incident photon energy of 0.95 - 1.05 GeV (top plots) 
and 2.15 - 2.25 GeV (bottom plots). The horizontal lines represent the 
fiducial limits derived from the Chauvenet's criterion. The distribution is 
practically flat and most of the runs yield values of A* whitin the limits 
of the fiducial interval (— Aq, Aq) represented by the two horizontal lines. 
Only two files have been discarded from the overall dataset 

3.16 Ratio of the normalized photo-disintegration proton yields from sectors 1 to 
5 over sector 6 as a function of the run number. As it can be seen, all the 
distributions are flat and no time dependence on the run period is shown. 
The average values of uniformity are reported in table 13.31 

3.17 Normalized total proton yield per T-counter as a function of incident photon 
energy: the different symbols refer to different running conditions (tagger 
prescale ON or OFF and incident electron energy Eq). The behavior of 
the normalized proton yield per tagger T-counter does not show evident 
dependence on the experimental settings 

3.18 Distribution of the event time difference At between the tagger time and the 
start counter time for all outgoing charged particles (top plot), for proton 
events (middle plot), and for photo-disintegration events (bottom plot). 

3.19 Level of accidental coincidences i? shown as a function of the run number. 
The relatively flat behavior has two different levels: one is relative to the 
tagger prescale OFF configuration and has an average value around 1% while 
the other is relative to the tagger prescale ON configuration and is less then 
2% 

3.20 Region of the kinematic plane where the momentum corrections could be 
calculated using two body kinematics 

3.21 Behavior of the ratio < > as a function of ^ in CLAS sector 1 for 
three values of the measured proton momentum P^: 0.55 GeV (top plot), 
0.65 GeV (middle plot), and 0.75 GeV (bottom plot). In each (/)-bin the 
ratio i?0 — P^l Pq is calculated and its distribution is fitted using a Gaussian 
function to obtain the average value < i?^ > 

3.22 missing mass distributions for — >■ pX events for sector 2 (top plot) and 
sector 3 (bottom plot) for two different proton momenta and scattering 
angles: P=0.65 GeV, 9 = 115° and P=1.15 GeV, 9 = 15°, respectively. The 
missing mass distributions are sharper after the application of the correction 
procedure and their peak positions are closer to the neutron rest mass value. 
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3.23 Before the momentum correction. Ratio Mx/Mn (being A/„ the nominal 
value for the neutron rest mass) plotted against the proton scattering angle 
in the LAB system separately for the six CLAS sectors and for the interme- 
diate proton momentum of 1.05 GeV. The deviation from the unity affects 
practically all six CLAS sectors, and shows a dependence on the proton scat- 
tering angle. The maximum deviation is of the order of 3% in some sectors 
(specifically sector 4 and 6) 

3.24 After the momentum correction. Ratio Mx/Mn (being M„ the nominal 
value for the neutron rest mass) plotted against the proton scattering angle 
in the LAB system separately for the six CLAS sectors and for the intermedi- 
ate proton momentum 1.05 GeV. The situation is clearly improved compared 
to the one shown in Fig. 13.231 since the values obtained for Mx/Mn are of 
the order of 1 and the dependence on the proton scattering angle is removed 
by the correction 

3.25 Before the momentum correction. Distributions of the Mx values obtained 
from the missing mass distributions for 7^ — )■ pX events. The top plot shows 
the overall Mx distribution for the six CLAS sectors, while the remaining 
six plots show the contributions of the different sectors. It is seen that, the 
Mx distributions exhibit large widths which are different in the six CLAS 
sectors 

3.26 After the momentum correction. Distributions of the Mx values obtained 
from the missing mass distributions for jd — > pX events. The top plot shows 
the overall Mx distribution for the six CLAS sectors, while the remaining six 
plots show the contributions of the different sectors. After the application 
of the momentum correction procedure, the widths of the Mx distributions 
are reduced and are more uniform in the six sectors 

3.27 Position of the missing mass distribution peaks (associated with their RMS) 
before (solid/black squares) and after the correction (solid/red dots). The 
horizontal line represents the nominal value of the neutron rest mass. . . 

3.28 The CLAS proton detection efficiency evaluated using GSim as a function 
of the azimuthal angle (j) for GeV, 9^^^ = 45°, and with a <j) 
binning of 5°. The shapes of the six CLAS sectors are reflected by the 
behavior of the proton detection efhciency and are clearly identified by the 
sudden decrease of the efficiency near the sector borders 

3.29 Behavior of the proton detection efficiency obtained from GSim shown as 
a function of the azimuthal and polar angles (on the abscissa and on the 
ordinates, respectively) of the generated protons for P^^^=0.95 GeV. The 
quota represent the proton detection efficiency values and are shown in dif- 
ferent colors/grays according to the scale indicated on the right side of the 
plot. The CLAS sectors are numbered from left to right. In each sector, a 
region of practically constant high efficiency (close to 95%) is clearly iden- 
tified. The vertical low efficiency zones among sectors correspond the torus 
magnet coils. The horizontal gaps in the efficiency function (along 6) corre- 
spond to to out of work time-of-flight paddles or drift chambers wires taken 
into account in the simulation using the GPP tool 



3.30 Missing mass distribution for the event 7c? — > pn^ X at an intermediate pro- 
ton momentum P^^^ = 0.95 GeV, an azimuthal angle (j) = 10° (correspond- 
ing to the center of sector 1), and for three different proton scattering angles 
6*^^^: 10° (top), 40° (middle), and 70° (bottom). The curves superimposed 
to the plots are Gaussian plus straight line fitting functions calculated in 
the interval (0.6 — 1.05) GeV. The background contamination is evaluated 
according to Eq. 13.231 

3.31 CLAS proton detection efficiency in CLAS evaluated from GSIM (black/full 
dots) compared with the results measured from the 7c? pn^p reaction 
(red/open dots) for P^^^=0.95 GeV and 0^^^ = 45°. The resuhing values 
for the proton detection efficiency are in good agreement. Note the different 
(f) binning. The datum corresponding to a sector border azimuthal angle 
(j)B obtained from the jd — > pTT^p data is roughly the average of the GSim 
results at ipB ~ 5° and 0^-1-5° 

3.32 Proton detection efficiency in CLAS measured from the — > pir^p reac- 
tion shown as a function of the proton azimuthal and polar angles for 
Pp^^—0-Q5 GeV. The quota represent the values of the proton detection 
efficiency in different colors/grays according to the scale indicated on the 
right side of the plot. The CLAS sectors are again numbered from left to 
right. The regions of constant high efficiency (close to 95%) can be clearly 
identified, even if the sampling rate Acf) = 10° corresponds to a lower res- 
olution. The vertical low efficiency regions among sectors are due to the 
presence the CLAS torus magnet coils. The horizontal gaps in the efficiency 
function (along 6) correspond to to out of work time-of-fiight paddles and 
drift-chamber wires 

3.33 Final comparison between the results obtained from the GSim evaluation 
of the proton detection efficiency using the 7c? — ?■ pn reaction (black/full 
dots) and the measurement from the data using the jd — > pn^p reaction 
(red/empty dots). The points represent the weighted averages of the de- 
tection efficiency over proton momenta in the interval (0.5— 1.0) GeV for 
the sectors central region (defined by the interval — ±10° assuming the 
sector azimuthal angle between (—30°, 30°)) and are shown as a function of 
the proton scattering angle 9^^^. It is clearly seen that the GSim result 
is well matched by the ^d — >■ pn^p result. In addition, the actual detector 
status is well reproduced by the simulations: this is shown by the dips in 
the proton detection efficiency found at the same scattering angle 9^^^ in 
both results 

3.34 Distribution of the ratio R = ^data between the proton detection efficiency 
measured from 7c? — > pn^p and simulated from ^d — > pn (defined only in 
the kinematic bins APp^^ and A9p^^ populated by both kind of data) for 
the sectors central region (defined by the interval = ±10° assuming the 
sector azimuthal angle between (—30°, 30°)). The parameters representing 
the Gaussian fit are shown in the legend. It can be seen that the mean value 
of the distribution is of the order of the unity while the spread is about 8%. 



3.35 An example of fiducial cuts in CLAS sector 2 for P^^^^1.15 GeV. Upper 
panel: for 9^^^ = 15°. Lower panel: for 9^^^ = 65°. The horizontal shaded 
area shows the allowed variation band eo ± 6eo while the vertical lines de- 
limit the fiducial region according to the behavior of the proton detection 
efficiency as a function of the azimuthal angles (f) 

3.36 Fiducial cuts in the CLAS sector 2 for P^^^=1.15 GeV. The boundaries 
of the fiducial (p regions (along the abscissa) in the sector are shown for 
each proton scattering angle (along the ordinates). The mean values of the 
efficiency within the region are indicated by the color/gray code shown in the 
legend. The dip found at 9^^^ = 55° is due to an out of work time-of-flight 
paddle, as shown also in Fig. 13.331 

3.37 Result of the calculation of the efficiency-acceptance function for the six 
CLAS sectors. The abscissa report the proton momentum while the ordi- 
nates run on the proton scattering angle (both considered in the LAB sys- 
tem). Quotas identify regions of constant values of the efficiency-acceptance 
functions (numeric values can be obtained from the color/gray scale shown 
on the left). As can be seen from Fig. 13.371 the behavior of the six CLAS 
sectors is pretty uniform, except for the holes introduced by out of work 
channels in some sectors 

3.38 Missing mass distributions for the six CLAS sectors, for the intermediate 
energy £.^=0.95 GeV, and a forward scattering angle 9p^ — 15°: the peak 
corresponding to the neutron rest mass is clearly identified. Each distri- 
bution is fitted with a Gaussian function plus exponential curve to take 
into account the continuous background. The red/dotted arrows identify a 
(—3(7, 3cr) interval around the peak position, being a the standard deviation 
of the Gaussian function 

3.39 Missing mass distributions for the six CLAS sectors for 9p^ — 125° and for 
an energy incident photon energy E-y=0.95 GeV. The distribution are wider 
and the neutron signal is less pronounced since at backward scattering angles 
the resolution of the CLAS spectrometer is lower. Nevertheless, the missing 
mass distributions are still well described by a Gaussian plus exponential 
function fit. The red/dotted arrows identify a (—3(7, 3a) interval around the 
peak position, being a the standard deviation of the Gaussian function. . . 

3.40 The behavior of the background contamination fc as a function of the proton 
scattering angle for the six CLAS sectors and for the intermediate incident 
photon energy E^—0.95 GeV. The contribution is around 4% and is rela- 
tively flat for scattering angles less than 90° 

3.41 Behavior of the background contamination fc as a function of the incoming 
photon energy for the six CLAS sectors and for = 35°. The last three 
data points, corresponding to the energy interval 1.8 < < 2.1 GeV, are 
obtained averaging the contributions from all sectors and are identical in 
the plots shown. It is clearly seen that the data are well reproduced using a 
linear fit 



3.42 Behavior of the background contamination /c as a function of the incoming 
photon energy for the six CLAS sectors and for = 135°. The last three 
data points, corresponding to the energy interval 1.8 < < 2.1 GeV, are 
obtained averaging the contributions from all sectors and are identical in the 
plots shown. At this proton scattering angles the background contribution 
is larger but it can still be extrapolated using a linear fit of the available 
points 



ni 



4.1 Distributions of the differences between the value of the cross section ob- 
tained in the individual sector and its average over the six sectors, divided 
by the average itself, as defined in Eq. 14.21 for each bin of proton scattering 
angle and incident photon energy. The result is shown separately for the 
six CLAS sectors (numbered from left to right and top to bottom). In each 
panel the statistical parameters of the distributions are shown and it can be 
seen that a systematic shift of the order of 2 5% is present in all sectors. 111? 



4.2 The CLAS results on the differential cross section ^ as a function of the 

at 

total energy s (where t and s are the momentum transfer to the proton and 
the total energy in the Mandelstam notation) for proton scattering angles 
in the range 15° < < 85°. As can be seen, the cross section decrease 
as a function of s is extremely steep and this behavior is enhanced as the 
scattering angle approaches 90°. The blue/solid line superimposed to the 
experimental results represents a power law fitting function proportional to 
s" 



4.3 The CLAS results on the differential cross section ^ as a function of the 
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total energy s (where t and s are the momentum transfer to the proton and 
the total energy in the Mandelstam notation) for proton scattering angles in 
the range 95° < 9™ < 155°. It can be seen that the cross section tendency 
to a fast decrease as a function of s, slows down going toward backward 
scattering angles. The blue/solid line superimposed to the experimental 
results represents a power law fitting function proportional to s" 120l 



4.4 Exponents n obtained by the angular survey on the cross section. The data 
behavior is of the type ^ ~ (const)s" as a function of the total energy s. It 
can be seen that the cross section decrease as a function of s follows different 
exponents n for different scattering angles. At very forward and backward 
angles (corresponding to the dominance of the t or u channel, and where u 
is the momentum transfer to the neutron) the exponent is of the order of 
—9 and —7 respectively. At more intermediate angles (around 90°), where 
both the t-dominance and it-dominance are suppressed, n drops around —13. 
A clear asymmetry exist between forward (n = —9) and backward angles 
{n = —7). The errors on the data points are those given by MINUIT on the 
fitting parameter n 1121 
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4.5 Angular distributions of the deuteron photo-disintegration cross section mea- 
sured by the CLAS (sohd/red dots) in the incident photon energy range 
0.55 — 1.25 GeV. It can be seen that the CLAS data cover a broad range 
in the proton scattering angle: 15° < < 155°. Results from Mainz [13] 
(open dots), SLAC [T71 [TB] (open squares), JLab Hall A (open dia- 
monds) and Hall C [24l [25] (open triangles) are included and are generally 
in agreement with the CLAS data. The angular dependence of the photo- 
disintegration cross section shows a pronounced asymmetry between the for- 
ward and backward proton scattering angles. For incident photon energies 
higher than 1 GeV, the solid line represents the non perturbative calculation 
of the QGSM ^31^ which well reproduces the data. In particular, the QGSM 
calculation accounts for the forward to backward asymmetry shown by the 
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4.6 Angular distributions of the deuteron photo-disintegration cross section mea- 
sured by the CLAS (solid/red dots) in the incident photon energy range 
1.35 — 2.05 GeV. Results from already published experiments such as SLAC [TT] 
[T8] (open squares), JLab Hall A [26] (open diamonds) and Hall C [24l [25] 
(open triangles) are included, where available. The cross section asymmetry 
between the forward and backward angles persists in this photon energy 
range. The QGSM [3T] calculation (shown as a solid line) well reproduces 



the data behavior also in this energy range [124 



4 . 7 Angular distributions of the deuteron photo-disintegration cross section mea- 
sured by the CLAS (solid/red dots) in the incident photon energy range 
2.15-2.85 GeV. Resuhs from SLAC [H] [H] (open squares), JLab Hah 
A [26] (open diamonds) and HaU C [24l [25] (open triangles) are included, 
where available. The angular distributions are well reproduced by the 
QGSM [31] calculation (shown as a solid line) which accounts for the data 



asymmetry between forward and backward scattering directions [125 



4.8 Deuteron photo-disintegration differential cross section expressed as da/dt 
multiplied by s^"'^ (being t the momentum transfer and s the total energy 
in the Mandelstam notation) in the angular range from = 15° up to 
= 85°. The CLAS results are represented by the sohd dots. Also shown 
are previously available data: SLAC [TT] [18] data are represented by open 
squares, JLab Hall C [23] [5S] data by open triangles, and JLab Hall A [5S] 
data by open diamonds. The arrows in each panel, indicate the photon 
energies where the onset of the s~^^ scaling is foreseen. This limit correspond 
to proton transverse momenta of 1 GeV/c. For the jd — > pn reaction the 
proton transverse momentum can be expressed — ^E^M d sin^ (9^^). 
According to the CCR, a constant behavior of the cross section as a function 
of the incident photon energy is expected when the limit for Pt is reached. 
The solid line represents the QGSM [31] non perturbative calculation for the 
cross section. This model, as discussed in the first Chapter of this thesis, 
gives a non perturbative description of the photo-disintegration cross section 
assuming that the scattering amplitude at very high energy is dominated by 



the exchange of three valence quarks in the t-channel [127 
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Deuteron photo-disintegration differential cross section expressed as da/dt 
multiplied by s^-^ (being t the momentum transfer and s the total energy 
in the Mandelstam notation) in the angular range from 9™ = 95° up to 
6»CM = 1550^ rpj^g CLAS resuhs are represented by the solid dots. Also 
shown are previously available data: SLAC [T71 [TB] data are represented by 
open squares, JLab Hall C [211 ES] data by open triangles, and JLab Hall 
A [2 6) data by open diamonds. The arrows in each panel, indicate the photon 
energies where the onset of the s~^^ scaling is foreseen. This limit correspond 
to proton transverse momenta of 1 GeV/c. For the 'yd — > pn reaction the 
proton transverse momentum can be expressed by = \E^Md [9^^). 
According to the CCR, a constant behavior of the cross section as a function 
of the incident photon energy is expected when the limit for Pt is reached. 
It can be seen that the asymptotic behavior for the cross section is well 
estabhshed from 6*™ = 75° up to 6*™ = 125° while at backward angles, the 
cross section dependence on the photon energy follows powers of s different 
from —11 as shown in Fig. 14.41 The solid line represents the QGSM 31 , a 
non perturbative calculation for the cross section 
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Abstract 



The study of two body deuteron photo-disintegration is ideal for investigating the 
transition region from meson exchange to the Quantum Chromo Dynamics (QCD) 
regime because of its simphcity and amenabihty to calculation [Ij. 

In the past years, the theoretical efforts have been focused in two different di- 
rections: extending at higher energies the models based on meson exchange, and 
improving the models inspired by QCD principles to extrapolate their predictions 
at lower energies, in order to find out the approach that best describes the experi- 
mental data and in which energy region the transition between the two descriptions 
takes place. 

Discriminating among the available models is still difficult because data are rel- 
atively scarce and mainly cover a limited angular range. 

The possibilities for a complete experimental investigation of the 'yd — )• pn pro- 
cess are much improved since the availability of an intense (200 ^A) electron beam 
of 6 GeV energy at the Jefferson Lab Continuous Electron Beam Accelerator Facility 
(CEBAF). In this new experimental context, a comprehensive study of the 'jd — t- pn 
reaction was started in recent years, using complementary equipments: 

• very high energy, high-intensity polarized and unpolarized bremsstrahlung pho- 
ton beams in conjunction with high resolution spectrometers of moderate an- 
gular acceptance for measurements at few scattering angles (JLab Halls "A" 
and "C"); 

• high energy (unpolarized) tagged photon beam in conjunction with the large 
angular acceptance (nearly Att) and good momentum resolution of the CLAS 
detector (JLab Hall "B"). 

This thesis reports on the first, wide-ranging survey of the jd — > pn differential 
cross section, for incident photon energies between 0.5 and 3.2 GeV, characterized 
by the high statistics and large kinematic coverage obtained with the CEBAF Large 
Acceptance Spectrometer at Hall B (CLAS). 

The subject is introduced from the physical motivation for the measurement 
together with an overview of the theoretical efforts made to describe the photo- 
disintegration process in the few GeV region. Then, the CLAS detector characteris- 
tics are briefly described in order to better illustrate the full data analysis procedure 
which, together with the final results and their theoretical interpretation, constitute 
the very original part of this work. 



Chapter 1 

Physical Motivations 



1.1 Introduction 



The study of high-energy two body deuteron photo-disintegration has received 
a renewed interest in recent years in the attempt to evidence quark related effects 
in nuclei. In fact, this reaction is well suited for studying the interplay of nuclear 
and particle physics while allowing the best possible separation between the known 
electromagnetic reaction mechanism and the nuclear system structure. 

The physics interest deals with how and at what energy the transition takes 
place from the hadronic picture of the deuteron, established in the reaction at low 
energies (below 1 GeV), to the quark-gluon picture, which is expected to be correct 
for energies much higher then 1 GeV, where distances of the order of few tenth of 
Fermi play a role. 

From the experimental point of view, the history of — )• pn measurements 
has been characterized by large discrepancies among different datasets. The situa- 
tion was complicated by the fact that the majority of the older experiments used 
bremsstrahlung photons and were subject to uncertainties due the poor knowledge 
of the photon flux and energy. Moreover, large solid angle detectors were not widely 
available, leading to relatively small coverage in the proton scattering angle. 

1.2 The Low Energy Region 

The total photo-disintegration cross section data at incident photon energy below 
1 GeV show a discrepancy up to 40% in the A(1232) resonance region. However, 
considering data from experiments performed with quasi mono-energetic photons or 
improved untagged bremsstrahlung techniques, a fair agreement (of order 8-^10%) is 
found especially above 300 MeV of photon energy [21 [3]. 

In this energy region the cross section data can be described with a number of 
theoretical models all based on the "traditional" nuclear picture of meson exchanges 
such as: the Impulse Approximation the Coupled Channel Approach the 
Diagrammatic Approach by [6], and the Modified Current Conservation [7]. 

The measured differential cross section is reproduced using a 4*^ order Legendre 
polynomial expansion: 



1 



4 
1=0 



(1.1) 



and the total cross section can be computed using the relation: atot = ^ttAq. 

A selection from the measured angular distributions for the 'yd — )• pn cross sec- 
tion in the photon energy range between 0.1 and 0.8 GeV is shown in Fig. 1 1.1 1 where 
the results are from the following experiments: Bonn f8], Frascati [9j, MIT [10] . 
LEGS [12j, and Mainz/DAPHNE |13j. The overall dataset cover a relatively wide 
range in the proton scattering angle 9p^, that is from 25° up to 160°. 

The most significant contribution to the cross section data in this incident photon 
energy region is provided by the Mainz experiment, performed with a tagged photon 
beam in the energy range from 0.1 to 0.8 GeV in conjunction with the DAPHNE 
large-angle spectrometer. Results from this experiment are shown in Fig. 11.11 to- 
gether with data points from Refs [141 [TSl I16j . It can be seen that below 0.4 GeV 
the DAPHNE data better match the results of experiments carried out using a 
tagged photon beam such as Bonn and LEGS. This peculiarity is also reflected by 
the Legendre coefficients Ai — A4 plotted in the left panel of Fig. 11.21 

The consistency between the angular distributions provided by different exper- 
iments is better than that of the absolute cross section shown in Fig. 11.31 The 
DAPHNE results show a noticeable change in slope of the cross section, moving 
toward higher photon energies. The effects starts around 0.5 GeV and it is reflected 
in the shape of the Legendre polynomial coefficients shown in Fig. 11.21 

In this energy region, also the Ai and A3 coefficients exhibit a small but signifi- 
cant deviation from the trend for the lowest energies: this could be explained by the 
intermediate excitation of baryon resonances in the second resonance region which 
includes Pii(1440), 1)13(1520), and S'ii(1535). 

1.2.1 Theoretical Interpretation 

Included in Figs 11.11 11.21 and 11.31 are the results of four theoretical models based 
on meson-nucleon degrees of freedom. 

The Impulse Approximation and the Coupled Channel Approach consider the 
nucleon internal degrees of freedom explicitly in the hamiltonian and the deuteron 
wave function is a superposition of different intrinsic configurations (NN, NA, AA) 
containing also excited nucleons. These configurations are mutually coupled by non- 
diagonal transition potentials. The two models need some kind of phenomenological 
input and start on the same assumptions, but handle the resulting set of coupled 
equations differently: the Impulse Approximation is a perturbative solution in which 
only the A^A^ ground configuration is retained while the isobar degree of freedom are 
only implicitly taken into account via the dispersive part of the N — N potential. 
On the other hand, the Coupled Channel Approach gives the full solution for the 
A^A^ — A A coupled channel using custom models for the different potentials involved 
and for the A width. Its results includes only multipoles up to L = 4 while for 
scattering waves with j > 4 it still relies on the Impulse Approximation. 

Figs 11.11 and 11.21 show that the Impulse Approximation reproduces the angular 
distributions up to 360 MeV even though it overestimates the total cross section. 
On the other hand, the Coupled Channel Approach gives a good description of 
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Figure 1.1: Differential cross section angular distributions for the deuteron photo- 
disintegration process for selected incident photon energies below 0.8 GeV. Data points 
are from the following experiments: Bonn [8] (open squares), Frascati [9] (open crosses), 
MIT [TU] (open triangles), LEGS IH] (open dots), and Mainz/DAPHNE [H] (solid dots). 
The curves belong to the following theoretical models: the Impulse Approximation [i] 
(dashed line), the Coupled Channel Approach [S] (up to 0.4 GeV) (continuous line), the 
Diagrammatic Approach [6] (dot dashed line), and the Modified Current Conservation [7] 
(dotted line). As can be seen from plots, the predictions of the theoretical models have 
problems in reproducing the data, starting from relatively low incident photon energies. 
(The plots are a courtesy of P. Pedroni - INFN, Pavia) 




Figure 1.2: Legendre coefficients for the differential cross section. At the highest photon 
energy reached at Mainz (around 0.8 GeV), the Ai and A3 coefficients exhibit a small 
but significant deviation from the trend for the lowest energies. Data points are from 
the following experiments: Bonn [8] (open squares), Frascati [9] (open crosses), MIT [10] 
(open triangles), LEGS [I^ (open dots), and Mainz/DAPHNE p3] (solid dots). The curves 
belong to the following theoretical models: the Impulse Approximation [4] (dashed line), 
the Coupled Channel Approach [S] (up to 0.4 GeV) (continuous line), the Diagrammatic 
Approach [6] (dot dashed line), and the Modified Current Conservation [7] (dotted line). 
(Courtesy of P. Pedroni - INFN, Pavia) 
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Figure 1.3: Total photo-disintegration cross section. In the energy region around 0.8 GeV 
the total cross section behavior is not a smooth decreasing function of the incident photon 
energy: this could be explained by the intermediate excitation of baryon resonances in 
the second resonance region which includes Pii(1440), Di3(1520), and 5*11(1535). Data 
points are from the following experiments: Bonn [5] (open squares), Frascati [5] (open 
crosses), MIT [TD] (open triangles), LEGS [12] (open dots), and Mainz/DAPHNE [H] (solid 
dots). The curves belong to the following theoretical models: the Impulse Approximation [3] 
(dashed line), the Coupled Channel Approach [S] (up to 0.4 GeV) (continuous line), the 
Diagrammatic Approach [6] (dot dashed line), and the Modified Current Conservation [7] 
(dotted line). (Courtesy of P. Pedroni - INFN, Pavia) 



the total cross section, but has some problems describing the shape of the angular 
distributions, since a dip is found at 90° which is not reflected by the data. 

This feature appears also in the behavior of the A2 coefficient shown in Fig. 11.21 
(upper left) which determines the curvature of the angular distribution around 90° 
and which is strongly overestimated above 160 MeV. This anomalous behavior is 
produced by the A^A configurations in higher partial waves [U [5] with j > 4 which 
are important for the differential cross sections, but negligible for the total cross 
section. 

In the Laget Diagrammatical Approach, the amplitude is expanded in a series 
of relevant interaction mechanisms. Each elementary transition amplitude is then 
determined starting from effective Lagrangians for the basic A^vr and couplings. 
This model is again similar to the Impulse Approximation but predicts an incorrect 
energy for the A peak excitation. This discrepancy is due to dynamical effects on 
the A mass given by the corresponding N — N and — A channels that are not 
taken accounted for, so that the resulting energy for the A peak is shifted. This 
effect would be compensated by the — A re-scattering, which is expected to shift 
back in place the A peak. 

The approach of Jaus, Bofinger and Woolcock, calculates the multipole ampli- 
tude from an effective current density which takes into account convection, spin, 
and meson exchange currents. This procedure is based on an effective one-pion ex- 
change current density which does not violate the current conservation law. The 
model allows for two adjustable parameters: the ttNN coupling constant and the 
pion cut-off mass, which can be used to fit the existing experimental data below 
200 MeV. This limit stems from the fact that only a limited number of multipole 
amplitudes have been considered and the finite A width has not been taken in to 
account. The calculations predict a total cross section which is very similar to the 
Coupled Channel result. Both methods agree with the Mainz cross section values 
up to about 120 MeV, but then between 120 and 180 MeV they both overestimate 
the experimental results. 



1.3 The High Energy Region 



For incident photon energies of the order of 1 GeV, the deuteron photo-disintegration 
differential cross section was measured for the first time by the SLAG NE8 experi- 
ment [17]. Results from the NE8 experiment are shown Fig. 11.41 for a proton scat- 
tering angle 6^^ = 90°. The data are multiplied by the s^^ factor, where s is the 
square of the total energy in the usual Mandelstam notation. 

Clearly, for > 1 GeV, the data shows a flat behavior. This represents the 
first evidence for the possibility of the so called scaling in the deuteron photo- 
disintegration cross section, suggesting that point-like constituents as quarks and 
gluons (i.e. the QGD degrees of freedom) could be the correct description for the 
process, in analogy to the Deep Inelastic Scattering (DIS). On the contrary, the tra- 
ditional meson exchange description, discussed in the previous section to interpret 
the low energy data, breaks down. 



Figure 1.4: Results form the NE8 experiment: the cross section shows scaling for incident 
photon energies larger then 1 GeV. The solid line represents a meson-nucleon description 
which do not account for the data starting from 0.5 GeV of incident photon energy. 

This striking result was supported and extended by another experiment subse- 
quently performed at SLAG for photon energies up to 2.8 GeV [IB]. 

The SLAG results have renewed the interest for the study of the deuteron photo- 
disintegration at high energy and as a consequence a vast physics program has been 
approved at Jefferson Lab. This program involves all the three JLab collaborations in 
a full investigation of the physics topics related to the deuteron photo-disintegration 
process: 

• Hall A Gollaboration: angular distributions at three medium energies and 
investigation of the polarization observables [19j ; 

• Hall B Gollaboration: nearly complete angular distributions from 0.5 to 3.0 GeV 
of photon energy pO] : 

• Hall G Gollaboration: differential cross section at few angles and very high 
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energy (up to 5.5 GeV) [21j. 



7 



1.3.1 Theoretical Interpretation 

The experimental data provided by SLAC, JLab Hall A, and Hall C are shown 
in Fig. 11.51 for the measured proton scattering angles. The cross section is again 
multiplied by s^^ in order to evidence the possible scaling behavior. 

In this energy region an asymptotic behavior could be expected on the basis of 
the simple model of the Constituent Counting Rules (CCR, [22j). This model con- 
siders the carriers of the currents within the hadrons as structureless and describable 
in terms of a renormalizable field theory, with the conditions of asymptotically scale 
invariance for the interactions among the constituents, and finiteness at the origin 
of the hadronic wave functions. Accordingly, applying dimensional counting to the 
minimum quark field component of a hadron accounts for many experimental conse- 
quences of its compositeness. The constituent counting rule prediction for exclusive 
scattering is: 

f.»<-./(i). a.2) 

where s — t- oo and t/s is fixed, where s, t (and u) are the usual relativistic invariants 
in the Mandelstam notation. Here n is the number of elementary fields in the initial 
and final states. 

Considering the 'yd — )• pn process it turns out that n = 13 since the initial state 
has one photon plus 6 quark fields while the final state has 6 quark fields which 
means that at at high incident photon energy and intermediate angles, conventionally 
90°, where both the t-dominance and u-dominance are suppressed the 'yd — )• pn 
differential cross section is expected to scale as s~^^. As shown in Fig. 11.51 this 
prediction is well confirmed by the experimental data for = 69° and 89°. 

The overall data behavior is less clear and may suggest a deviation from the 
predictions of the simple CCR scaling indicating, in turn, that pQCD itself is not 
plainly applicable. 

In fact, the cross section at the forward angles = 52° and 36° falls off more 
slowly, with an ~ scaling at lower energies, until the onset of the s~^^ behavior 
at about 3 and 4 GeV beam energy, respectively. However it results that an absolute 
pQCD calculation would not correctly predict the data. Further observations may 
be made for other photo-reactions, and it remains to be seen how this behavior arises 
and if pQCD is a genuine explanation for it. 

In this different approach have to be used and it seems very reasonable 

that a soft description for the process will have to be introduced. 

This can be done following opposite strategies, such as the inclusion of QCD 
degrees of freedom in the low energy description or the extrapolation of the conven- 
tional N — TT interaction mechanisms to the higher energy region. It is also possible 
to merge suitable parts of this two techniques to obtain an hybrid description (for 
a full review see Ref . [23] ) . The curves shown in Fig. 11.51 represent the predictions 
of several theoretical models implementing the outlined strategies. In the Reduced 
Nuclear Amplitude (RNA) model [271 [28] the binding of the quarks inside the nu- 
cleons and the deuteron is taken into account with empirical form factors and the 
elementary cross section is computed assuming CCR scaling. This approach is able 
to describe the ^d — t- pn cross section with appropriate normalization factors for 
6)CM = ggo j^j^j ggo ^-^j^ R > 2 GeV. 
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Figure 1.5: Deuteron photo-disintegration differential cross section multiplied by s^^. Ex- 
perimental data are from Mainz [13] (open/magenta circles), SLAC [171 El] (open/green 
squares), Jlab Hall C [24j[25] (open/blue triangles), and Hall A [26] (open/black diamonds). 
Points are multiplied by s^^, and arrows indicate the expected threshold for the onset of the 
CCR scaling, where the momentum transfer to the proton, t is of the order of 1 GeV^ [22] . 
Theoretical curves are described in the text. 




In the Hard quark Rescattering Model (HRM) [29], the elementary interaction 

consists of a quark exchange between the two nucleons. The incoming photon is 
absorbed by a quark of one nucleon which then gives up its momentum via a hard 
gluon exchange with a quark of the other nucleon. The model assumes that this 
rescattering mechanism is analogous, with some approximation, to the wide angle 
p — n scattering which is also dominated by quark exchange. The limits for the 
applicability of the model are > 2.5 GeV and momentum transfer t >2 GeV^, 
but, under particular assumptions for the short distance p — n interaction, they can 
be extended. However, the agreement with the data is poor, especially at higher 
energies where the uncertainty is large, due to the limited knowledge of the pn cross 
section for the actual kinematic conditions. 

Even if traditional models based on the N — n picture fail to reproduce the data 
at E"^ > 1 GeV, this approach can be extended, with the appropriate modifica- 
tions, in the few GeV region. The Asymptotic Meson Exchange Current (AMEC) 
model [30], uses form factors to describe the d — NN interaction vertex and an over- 
all normalization factor fixed by fitting the experimental data at 1 GeV. The results 
reproduce the energy dependence of the cross section only for = 89°. 

It is worth noticing that this non QCD-based model provides a scaling law for 
the cross section, with an exponent depending on the scattering angle. 

A non perturbative description of the photo-disintegration cross section data is 
given by the Quark Gluon String Model (QGSM) (Ref. |31j ) This model assumes that 
the scattering amplitude at very high energy is dominated by the exchange of three 
valence quarks in the i-channel. The duality property of scattering amplitudes [32| 
[33] allows its extension at intermediate energies. 

As shown in Fig. 11.51 the predictions of the QGSM are in good agreement with 
the available data in the full energy range at ^ ggo^ ggo_ the trend of 

the QGSM theoretical curve could be compatible with the last two points at the 
highest energies since their errors are quite large. At = 36° the situation is 
more delicate even if the QGSM curve has a better trend with respect to the other 
calculations. 

A better insight and a more effective discrimination among the predictions of 
the competing theoretical models are obtained examining the angular distributions 
of the differential cross section I34i • 



Experimental data recently obtained from [26] for three different incident photon 
energies: = 1.6, 1.9, and 2.4 GeV are shown in Fig. 11.61 

These angular distribution were measured only for the indicated energies and are 
not complete since the very forward and backward regions are not covered. Addi- 
tional data from Refs |18[ [Ml [25] are also shown in Fig. 11.61 together with the curves 
representing theoretical predictions. The shaded area is the HRM [35] calculation 
which has an uncertainty of 15-20% while the QGSM prediction is represented by a 
dashed line. Other calculations are not shown since they were not provided (AMEC) 
or the model could not be applied in this case (RNA) since the curves for different 
angles have arbitrary normalizations. 

Also in this case, the QGSM well reproduces the available distributions while 
the HRM (apart from a possible overall shift) show a different behavior at small and 
large angles. 

From the experimental side, it is clear that the knowledge of the full angular 
distributions, including the points at very forward and backward angles and in a 




in 




Figure 1.6: Angular distributions for three incident photon energies: 1.6 (top plot), 
1.9 (middle plot), and 2.4 GeV (bottom plot). The experimental data are from Refs [HI 
while the curves represent the theoretical predictions of the QGSM [3T] and the 
HRM ^35_. Prediction from other models are not reported since the were not provided or 
could not be calculated. From the plots it is seen that the QGSM well describes the data. 
This model also predicts a forward to backward asymmetry. More data especially at very 
forward and backward angles are necessary to verify this prediction. 
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wider incident photon energy range will be extremely helpful to better discriminate 
among the models and to further check the QGSM predictions. 

1.3.2 The Origin of Cross Section Scaling 

It has to be said that the applicability of the CCR at this incident photon en- 
ergies is still an open question. In fact, this approach is based on the assumption 
that soft and hard contributions factorize |36j even if this hypothesis (relying on 
pQCD) may be valid only for very high momentum transfers which are not at- 
tained in the few GeV energy region where experimental data are available (see 
Refs [371 [Ml EHl HO]). In order to clarify this point it is useful to consider the 
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Figure 1.7: Proton polarization transfers Cx (transverse in-plane, left panel) and Cz 
(longitudinal in-plane, right panel) in the deuteron photo-disintegration at — 90°. The 
y-axis is perpendicular to the scattering plane, the z-axis is parallel to the scattered proton 
momentum and the x-axis is perpendicular to y and z. Experimental data are from |41] 
and have been corrected for spin rotation due to the LAB — >■ CM transformation. The 
solid/red curve represents the QGSM prediction for [Hj- The QGSM calculation for Cx 
is underway. 



complementary information to the differential cross section given by polarization 
observables. Significant data on this topic were recently provided by the JLab Hall 
A Collaboration [41] and are shown in Fig. 11.71 which gives the proton in-plane po- 
larization transfers and the photon asymmetry as measured in the — )• pn reaction 
for the proton angle = 90°. In this case, the signature of pQCD effects rests 
in the so called Hadron Hehcity Conservation (HHC) [Ml Hll HH |45] . HHC arises 
from the fact that vector interactions such as photon-quark or gluon-quark couplings 
conserve chirality, leading, in turn, to conservation of the sums of the components 
of the hadronic spins along their respective momentum directions. 

If pQCD can be really invoked to explain the cross section scaling at = 90° 
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it is expected that at the same angle the polarized observables will satisfy the pQCD 
based prediction of HHC of Cx = Cz = 0. 

As shown in Fig. 11.71 for > 1.5 GeV the results on Cx seem not in agreement 
with the HHC prediction while the those on Cz are more difficult to interpret but 
the QGSM calculation for Cz has a reasonable trend with respect to the JLab data. 
A similar calculation for Cx is underway. 

The problem with HHC lies in the limiting hypothesis of chirality conservation 
while the QGSM, being a non perturbative approach, is far from the chiral limit 
(characteristic of pQCD) and is able to introduce a chirality violation which may be 
reflected in the data. 

For photon energies below 1.5 GeV the proton polarization transfers Cx and 
Cz show quasi-resonance structures at ^/s =2.7 and 2.9 GeV, respectively. Such 
structures are conceptually out of the QGSM range of applicability [46]. 




Figure 1.8: Asymmetry E for linearly polarized photons in the deuteron photo- 
disintegration at = 90° as a function of the incident photon energy. Experimental 
data are from Ref. [47]. The solid/red curve is the prediction of the QGSM [42] . 



Another observable carrying information on HHC predictions is the photon asym- 
metry, defined as 

^^^^ ^ - d^ii + da^ ' (^-^^ 

for which the HHC prediction at fl^'^ = 90° is 1 (see Ref. [32]). 

As can be seen from Fig. 11.81 most of the experimental data are limited to the 
incident photon energy region below 1.5 GeV where the perturbative regime may be 
not fully established, in any case the trend of the data does not seem to indicate that 
the HHC prediction is verified. In the higher energy region where > 1.5 GeV, 



the QGSM predicts a slow decrease of S(90°) Unfortunately the overlap with the 
available data is quite limited so definite conclusions cannot be drawn yet. 

Summarizing: at least for = 90° the cross section data and the in-plane 
polarization transfers results do not appear consistently described in a coherent 
perturbative picture. 

This indicates that a non perturbative approach may be needed in order to fully 
describe the deuteron photo-disintegration process (polarized and unpolarizcd). The 
QGSM gives a non perturbative description of the process which reproduces fairly 
well the available data: for this reason the QGSM foundations will be given in the 
next Section. 



14 



1.4 The Quark Gluon String Model 



Quantum chromodynamics is the fundamental theory of strong interactions. The 
asymptotic freedom of QCD imphes that its couphng constant as{Q'^) becomes small 
at small distances (large Q^) and this gives the possibility to use perturbation theory 
in processes with large momentum transfer. Experimental data for these processes 
are in good agreement with pQCD predictions. On the other hand, at large dis- 
tances (of order 1/Aqcd, where Aqcd — 0.3 GeV) the coupling constant is not 
small anymore and non perturbative effects, which are responsible for confinement, 
come into play. The large distance dynamics is necessary to understand processes 
with small momentum transfer, which give dominant contributions to high energy 
hadronic interactions. 

In this region a non perturbative approach should be used and there are not many 
ways to introduce non perturbative effects in QCD. One possibility is to consider 
lattice calculations [38], which are most directly related to the exact QCD results. 
This method, however, needs very large lattices to obtain accurate and reliable 
results and it is difficult to apply to complex phenomena as exclusive reactions at 
high energies. Another possibility is from the QCD sum rules [49] which express 
physical quantities in terms of several parameters such as vacuum correlators of 
gluonic and quark fields. As well as lattice calculations, this method has been 
applied mainly to static properties of hadrons. A third method, which will be briefly 
introduced in this overview, is the (or topological) expansion in QCD |50i I5H 
[52], where N is the number of colors or flavors. Dynamical content for this model 
will be drawn from the string (or color tube) models [53l [Ml ISHl [Ml [57j, which 
are closely related to the space-time picture of interaction in 1/N expansion. The 
string-like configuration of gluonic fields is confirmed by lattice calculations. 

These ingredients form the basis for the Quark Gluon String Model (QGSM, [58^ 
159] ) which then accounts for the confinement of quarks and can be considered a 
microscopic theory underlying the Regge phenomenology. As a result many rela- 
tions between parameters of reggeon theory (based on analyticity and unitarity) 
can be established from the QGSM [60], since a quark-gluon string can be identi- 
fied with a corresponding Regge trajectory. The QGSM has been applied to many 
different problems of strong interactions: hadronic mass spectrum |61j . widths of 
resonances [62], relations between total cross sections, residues of Regge poles [63] . 
behavior of hadronic form factors and, in recent years, to heavy ion collisions at high 
energy ^1] and to the deuteron two body photo-disintegration [31] . 

1.4.1 The Topological Expansion in QCD 

Large distance dynamics of QCD defines what is called soft process, i.e. when 
the value of the coupling constant ~ 1 and cannot be used as a good expansion 
parameter. Since QCD has not yet been solved exactly, it is necessary to find out 
another parameter, which is small enough to allow a non perturbative expansion of 
amplitudes. 

It is possible [50j to consider as a small parameter the quantity 1/Nc, where Nc 
is the number of colors with corrections to the main term of order l/N,? ~ 0.1. 
Another suitable parameter is 1/-/V/ where A'^^ is the number of fiavors [5l]. Soft 
processes involve only the light quark fiavor (u, d, s) contributions, so thatNf ~ 3. 



An intuitive justification for the use of these quantities, hes in the fact that these 
two parameters {\/Nc, are the only dimensionless quantities in massless QCD, 

being the only suitable quantities to be used as expansion parameters in that limiting 
case. The formal limit oi Nc ^ oo [50j {Nf/Nc — )• 0) shows many interesting 
properties and has been intensively studied since in this limiting case there is hope to 
obtain exact solution of the theory. However, this approximation is too far from the 
real world (for instance, all resonances are considered infinitely narrow (F ~ 1/Nc))- 
The case when the ratio {Nf/Nc — ^ 1) is fixed and the expansion of all amplitudes in 
l/Nf or l/Nc is carried out [51] is more realistic and is called topological expansion, 
since each term of the series corresponds to an infinite sets of Feynmann diagrams 
with definite topology. The 1/N expansion has a dynamical character which means 
that it can work better in some regions of the kinematical variables, while in some 
other regions its convergence can be worse. For example, in the case of a binary 
(2 —7- 2) amplitude, which is a function of two variables {s,t), the expansion 
becomes more accurate as t increases in the positive t region [65], while for high 
energies additional terms have to be considered. 

The topological classification of diagrams, in this scheme, starts from the ampli- 
tude for the binary reaction ab — )• cd. Here, the main contribution is given by the 
"planar" diagrams, shown in figure Fig. 11.91 where planar means a class of diagrams 
where all lines lie on a plane without hard interactions inside the contour delimited 
by the external quark lines. These lines of valence quarks determine the border of 
the diagram. 




Figure 1.9: Left panel: planar diagram for the reaction ab cd. Right panel: the same 
for the reaction a6 — > X. Full lines denotes quarks, wavy lines, gluons. 

All planar diagrams give a contribution of the same order in 1/A^, taking into ac- 
count that each color trace gives a factor Nc and each fermion loop a factor —Nf, 
and considering the limit of large N under the condition g^N ~ 1. It should be 
emphasized that the topological expansion can be applied only to colorless ampli- 
tudes. The unitarity content of the planar diagrams can be reveled cutting such a 
diagram in the s-channel as shown in right panel of Fig. 11.91 This operation gives 
the amplitude of multi-particle production, which is planar by itself. 

Cutting of the diagram on the left panel in Fig. II. 9^ either in the s or the t- 
channel, always contain a qq pair, which enters in the boundary of the diagram 
and determines the internal quantum numbers (charge, isospin, etc.) exchanged in 
the given channel. At high energies and small momentum transfers these planar 
diagrams are determined by exchange of Regge poles in the t-channel. Diagrams 
for elastic scattering [i.e. reactions without quantum numbers exchange in the t- 
channel) , are shown in 11.101 , where the valence quarks of the colliding hadrons are 
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conserved in the process of interaction. 




Figure 1.10: Left panel: Cylinder type diagram. Right panel: cutting of these diagrams in 
the s-channcl. 

These diagrams are of cylinder type and their s-channel cutting corresponds to 
the multi-particle production configurations. This is shown by the picture on the 
right of Fig. 11.101 which corresponds to the production of two chains of particles 
with the same structure as the one shown for the planar diagram cut of 11.91 From 
the t-channel point of view, the cylinder diagrams are due to gluon exchange in 
the i-channel, and they are usually assumed to give rise to a Pomeron pole. More 
complicated diagrams can be drawn, such the one shown in Fig. II. 11^ which has two 
Pomeron exchange contribution, and one hole in the surface. 




Figure 1.11: Topological diagram with one handle corresponding to the exchange of two 
Pomerons in the t-channel. 

Holes in a diagram defines "handles" , so each topological class of surfaces is char- 
acterized by a given number of boundaries (n^) and handles (n/i). The topological 
expansion allows a complete classification of diagrams based on the two parameters 
(rifc) and (n/j) and to determine their dependence on 1/iV. 

All the diagrams belonging to a given topological class have the following depen- 
dence on 1/A^: 

T{nb,nh)=^-j (1.4) 

so that the contribution of planar diagrams {rib = ^,nh = 0) to the scattering 
amplitude is 1/A^, the cylinder diagram (nf, = 2,n/j = 0) contribution is ~ (1/A^)^ 
and that of the diagram of Fig. 11.111 (n^ = 2, n/^ = 1) is (l/AT)^. 

This equation is valid for 4-point amplitudes; amplitudes with a larger number 
of hadrons should take into account that each external hadronic state introduces a 
factor due to the normalization of its wave function. The ratio of the cylinder 

to the planar diagram is ~ 'i^/N, however for amplitudes with vacuum quantum 
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numbers in the t-channel, the cylinder type diagrams will dominate at high energies. 
This because its relative contribution increases as s°^{o)~°fl(o) a,s energy increases 
as it is implied by the dynamical character of the 1/N expansion. In many cases, 
the type of the process (and the number of boundaries) is fixed by the quantum 
numbers in the t-channel. In this case the expansion parameter depends only by the 
number of handles and, according to Eq. 11.41 is (1/-^)^- 

I. 4.2 Color Tube or Quark Gluon String 

The topological expansion gives a useful classification of all QCD diagrams but 
needs a definite space time picture to become predictive. This picture is based on 
some assumptions on the properties of confinement which are related to the proper- 
ties of the QCD vacuum. Many reasons suggest that QCD vacuum has a complicated 
structure and contains important large scale fluctuations of gluon and quark fields. 
These fluctuations generate vacuum correlations of gluon and quark fields which are 
not described by perturbation theory. In particular, the vacuum correlator of the 
square of gluonic field C^^^ has been found from the analysis of QCD sum rules |49j: 

^(0|G;iJO)«0.012(GeV)\ (1.5) 

and it can be interpreted as a vacuum condensate connected to the energy density 
of the vacuum. 

= -4(0|^G^,G;1J0) with b=^N,- ^Nj . (1.6) 

Thus the energy density of the physical vacuum with the gluon condensate given by 

II. Sl is lower then for the "empty" vacuum. If the vacuum fluctuations will be broken 
in some region of space, then the energy of this region will be increased. Such a 
breaking of vacuum fluctuations in QCD (partial or complete) arises naturally inside 
hadrons, due to color fields induced by valence quarks. As a result, there will be 
two different states of vacuum condensate, one is the vacuum outside the hadrons 
and the other one inside the hadron. It is assumed that the chromo-electric field 
of quarks do not penetrate far into the ordinary vacuum medium since the gluons 
carry the color charge so that the flux lines attract each others. This behavior leads 
to the property of confinement, so that isolated states in the vacuum can exist only 
as "white" states with all color lines closed (see Fig. I1.12p . In this picture a hadron 
can be considered as a bubble in a "vacuum liquid". 




Figure 1.12: Hadron as a "bubble" in the vacuum. 
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The process of interaction of such bubbles at high energies (Fig. I1.13P leads to 
production of new objects: color tubes or strings (a string is a limiting case of a 
color tube with no transverse momentum distribution). 



It is necessary to consider separately two cases. 

• The processes with the annihilation of valence quarks of the colliding hadrons, 
corresponding to the planar diagrams of the expansion. 

• The processes of the diffraction type, where the valence quarks are conserved, 
connected with cylinder type diagrams. 

In the first case, after an interaction and annihilation of valence quarks the 
configuration shown in Fig. 11.131 (b), where the color lines connect the spectator 
quarks moving in opposite directions (in the CM system), is realized. As they move 
apart the region where the vacuum fluctuations are broken enlarges (Fig. 11.131 (c)) 
and the color tube or string is produced and this leads to an increase of the energy. 
At some point, the breaking of the color tube will be energetically favorable and 
the production of qq pairs from the vacuum (Fig. 11.131 (d)) will occur. This process 
repeats until many white bubbles or hadrons will be produced. The time needed 
for production of a hadron with momentum p and mass m is r ~ so in the 
CM system the fastest hadrons, which contain the spectator quarks qa and qb, will 
be produced last. Each produced q and q pair has a small relative momentum in 
their rest system which means a small rapidity difference. As a result of the Lorentz 
invariance, finally produced hadrons at high energy will be uniformly distributed in 
the rapidity and will have limited transverse momenta. 

1.4.3 Non Perturbative Approach in the QGSM 

The QGSM describes rather well the experimental data on exclusive and inclusive 
hadronic cross sections at high energy [HI |50l El [66] and gives new important 
insights into the Regge phenomenology. It can explain the quark-gluon content 
of different Reggeons, relate their residues, predict inclusive spectra of different 
particles in different kinematical regions, find the relations between exclusive and 
multiple production amplitudes through the unitary conditions. Concerning this 





Figure 1.13: Interaction and decay of a color tube. 
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last topic, multiple production at high energy, the QGSM is very similar to the dual 
parton model [67] . 

The Regge pole analysis of the exclusive hadronic reactions j68j was formulated 
much before QCD and for this reason is very often seen as a dormant topic being so 
strongly identified with the pre-quark era of the 5-matrix and dispersion relations 
approach to strong interactions. This point of view is not correct [69] since the 
description of high energy behavior in terms of singularities in the complex angu- 
lar momentum plane of two body scattering amplitudes is completely general being 
based on the general principles of unitarity, analiticity and crossing symmetry. In 
recent years it has been resurrected by all the theories which aim to a string formula- 
tion for hadrons and the conjecture of Regge that high energy behavior of two body 
scattering amplitudes might be describable in terms of moving poles in the J-plane 
is experimentally verified and all hadrons do lie on trajectories which connects the 
particle mass with its total angular momentum (Fig. I1.14p . 




Figure 1.14: Simple linear baryon Regge trajectories separated by A J = 2 [70] . 

Despite the very broad applications of such approaches to the different hadronic 
reactions and to the inclusive hadron-nucleus collisions at rather high energy they 
were poorly used in the description of exclusive reactions with nuclear targets. One 
of the main reason is apparently related to the fact that those approaches were orig- 
inally introduced to describe high energy data, but data on exclusive reactions with 
nuclear targets have been limited, up to now, mainly to the intermediate energy re- 
gion. Let us say something more on the region of validity in s and t of this approach. 
Prom the quantum mechanical point of view, high energy means that the wavelength 
of the photon is much smaller than the radius of the target nucleon which in QCD 
can be taken as the radius of confinement Rq, which is about 1 Fm. Therefore, at 
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>1 GeV, the necessary condition E^Rq ^> 1 can be considered as fairly satisfied. 
In the Regge approach, high energy means s » tchar where tchar is some character- 
istic value of t which is usually taken as ~ 1 (GeV)^. As the region of validity in t 
or u is concerned, there is a common convention that this is limited by |i| or \u\ < 1 
(GeV/c)^. This convention is based on the consideration that with increasing of \t\ 
or \u\ the secondary singularities in the complex angular momentum plane as daugh- 
ter Regge trajectories or branch points [7T] may also become important. However, 
this considerations cannot be formulated in a quantitative way because the residues 
of Regge poles and discontinuities of Regge cuts are usually parameterized and are 
different for different reactions. This means that at the phenomenological level, the 
description of data include contributions of some effective Regge trajectories which 
in turn may include the effects of secondary singularities. 

A different philosophy is used in the QGSM since each graph in this model 
is classified according to its topology and the corresponding Regge singularity is 
fixed by the quark-gluon content of the graph itself. Therefore the corresponding 
amplitude can be considered for all values of s and t as an analytic functions of those 
variables. The QGSM gives the possibility to calculate parameters that before were 
considered only at the phenomenological level in the Regge theory and to establish 
a bridge between QCD and the Regge approach to strong interactions. 

1.4.4 QCD Content of Regge Trajectories 

There is experimental evidence and theoretical reasons that both meson and 
baryon states should be connected, on the J vs plane by an essentially non 
linear relation (trajectory) instead of the simple scheme shown in Fig. 11.141 In fact, 
the starting point for Regge's work was the partial wave expansion of the scattering 
amplitudes [7l] but here we would like to reinterpret the original idea of Regge pole 
and trajectory using the modern QCD concepts and language to make evident that 
the Regge trajectories must contain some non linearity. 

Let us consider a meson state composed of a qq pair connected by the color flux 



Figure 1.15: Color strings: (a) massless string, (b) homogeneous string connecting a qq 
couple, (c) inhomogeneous string with variable tension due to color vacuum polarization. 



In a simple potential description the energy can be written as: 



tube fFig. fTTSl) . 
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where V{r) = or is a linear potential and a ~ 400 (MeV)^. 

For a sake of simplicity, (Fig. 11.151 (a)) massless quarks will be considered. In the 
CM system the energy will he H = 2p + ar where the orbital angular momentum 
't will be 7" X so that: 

H{r)=2-+(Tr . (1.8) 
r 

From the relation ^ = — ^ + cr = the minimum of H is found at tq = -^/^- If this 
position is chosen as the ground state of the system, the energy can be expressed 
as: ff(ro) = E = 2^/2^^^^ or E'^ = 8cj£. 

This relation gives the dependence of the angular momentum on the energy of 
the system: I = j^E'^. Performing the substitutions a ^ £ and E"^ =^ t: 

ait) = , (1.9) 

so that the expression for a simple linear Regge trajectory can be readily recognized. 

The choice of a massless qq system gives rise to a linear relation in the Mandel- 
stam variable t. The exact calculation yields a{t) = = a't for our string with 
massless ends where a'(0) ~ 1 GeV^. 

This scheme can be generalized to a string connecting two massive quarks (Fig. 11.151 
(b)). In this case a static tension is introduced on the color string, and an additional 
term is collected in the Regge trajectory: 

Q(t) = a(0) + a'(0)t . (1.10) 

Another correction term can be obtained considering the screening effects introduced 
by the vacuum polarization in QCD which complicates the string structure. A 
distance dependent tension (Fig. 11.151 (c)) can account for this effects since the 
formation of qq pairs from the vacuum weakens the flux tube opening "holes" in the 
string. A possible form for such a potential due to qq loops comes from lattice QCD 
calculations (see Ref. (72]) and can be parametrized in terms of the screening length 
/.-I: 

-- + ar]^ , (1.11) 

r / fj,r 

where a = 0.21 ± 0.01 Fm. 

This potential has the correct Coulomb behavior at small distances and ap- 
proaches a constant at large distances. The ratio ^ represents the splitting energy 
of the qq pair, with jj''^ = 0.90 ± 0.20 Fm. 

As was shown in the simple model described above, it is expected that com- 
posite hadronic systems such mesons and baryons give rise to (slightly) non linear 
Regge trajectories. In the QGSM, the non linearity of baryon Regge trajectories in 
incorporated according to Refs [73l [TH [75l [76] . For small momentum transfers, the 
non linear part is negligible, but when t > l(GeV/c)^ effects introduced by the non 
linear pert of the Regge trajectories become critical and must be considered. 
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1.5 Quark Hadron Reactions in the QGSM 



The QGSM was originally developed to describe hadron-hadron collisions at high 
energy. In order to better define its extension to photo-reactions a brief review of 
its application and predictions for hadronic collisions will be given. 

Reactions as pp — t- dir'^ and pd — )■ pTr~ are dominated by diagrams with three 
valence quark exchanges in the t-channel [77j in analogy to the deuteron photo- 
disintegration case (see Ref . [59] and Refs [78l [79] ) . 

Fig. 11.16] shows the total cross section data for pp — )■ dTr"*" [80]. It is seen that the 
QGSM very well reproduces the data with the single exception of the near threshold 
region. The QGSM describes quite well also the data from the pd — )• p7r~ reaction 
which have been obtained from the reversed channel Tr~p — )• pd [ST\ I82j . 
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Figure 1.16: Total cross sections of the reactions pp — > dn^ and pd — > pir^ as a function 
of the CM energy, and predictions of the QGS model. 



1.5.1 Transition Probabilities and Transition Amplitudes 

Each diagram of the topological expansion has a simple interpretation within 
the framework of the space-time pattern formulated in terms of a color string (see 
Refs [66','83j). For instance, at high energy ^/s, the binary reaction shown in Fig. II. 171 
(a) occurs due to a specific quark configuration in each pion, where (in the CM 
system) one quark (or anti-quark) takes almost the entire hadron momentum and 
plays the role of a spectator, while the valence anti-quark or quark is rather slow. 

The difference in the rapidities Ay between the quark q and the anti-quark q in 
each pion is 

Ay = y,-y^c^\ln- , (1.12) 
2 So 

with the scale sq — l(GeV)^. 
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When the two slow valence particles q and q from the two initial vr s annihilate, 
the fast spectator quark and anti-quark continue to move in the previous directions 
and form a color string in the intermediate state. After that, the string breaks to 
produce gg-pairs from the vacuum and the formation of the tt^tt" system takes 
place in the final state. The same space-time pattern holds also for the situation 
depicted in Fig. 11.171 (b) with the only difference, that the string is formed after the 
annihilation of a diquark anti-diquark pair from the NN system in the initial state. 
Correspondingly, the graph of Fig. 11.171 (c) shows the formation of the qq string due 
to annihilation of the valence diquark anti-diquark pair in the initial state and the 
production of a diquark anti-diquark pair as a consequence of the breaking of the 
string. The annihilation of the initial qq pair takes place, when a gap in rapidity of 
the valence q and q is small and both interacting partons are almost at rest in CM 
system so that the relative impact parameter b_|_ — bg^ is less than their interaction 
radius. It is possible to prove that the probability to find a valence quark with a 
rapidity yq at impact parameter inside a hadron can be written as (|58 t I6 H I83j) 



w {yq - 2/0, - bo±) = 

where c is a normalization constant, yo is the average rapidity, boj. is the transverse 
coordinate in the CM system in the impact parameter representation. Furthermore, 
it is possible to relate the parameter (3 and the effective interaction radius squared 
i2^(s) in ll.l3l that specify the quark distribution inside a hadron, to the phenomeno- 
logical parameters of a Regge trajectory ai{t) which gives the dominant contribution 
to the amplitude for the considered planar graph. In this case one gets 

R^{s) = Rl + a[{yq-yo) . /3 = 1 - a,(0) , (1.14) 

where a[ = a^(0) is the slope of the dominant Regge trajectory. 

Due to the creation of a string in the intermediate state the amplitude of a binary 
reaction ab — )■ cd has the s-channel factorization property, i.e. the probability for 
the string to produce different hadrons in the final state does not depend on the type 
of the annihilated quarks and it is only determined by the flavors of the produced 
quarks. The same independence also holds for the production of the color string in 
the intermediate state from the initial hadron configuration: it depends only on the 
type of the annihilated quarks. This s-channel factorization has been formulated in 
Refs [53 [m [83] in terms of transition probabilities as defined in Eq. (jl.l3p . 

An example of this can be seen for the reactions shown in Fig. 11.171 as well as 
for the crossed reactions 7r"'"7r~ — ?■ vr'^Tr'^, 7r"'"7r~ — )• NN , and NN — t- NN . 

For large s values and finite values of the 4- momentum transfer squared t, these 
channels are described by planar diagrams with t-channel valence-quark exchanges 
as in Fig. 11.171 where the single and double solid lines correspond to valence quarks 
and diquarks, respectively (soft gluon exchanges are not shown). Introducing the 
hadron(7)-quark and quark-hadron transition amplitudes, analytic expressions for 
the form factors have been found in the time-like region which admit continuation 
in the space-like region. 
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Figure 1.17: Planar diagrams describing the binary reactions 7r°7r° —)■ tt+tt (a: left panel), 
NN -)■ TT+TT" (b: central panel), and NN -)■ iViV (c: right panel). 



This approach can be generahzed [84] introducing the amphtudes T°-''^'^^{s,h±) 
and T'^^^'^'^{s,h±), that describe the formation and the fission of an intermediate 
color string. The amplitude for the reaction ab — )• cd described by the planar 
graph can be written, using the s-channel factorization property, in the form of a 
convolution of two amplitudes: 



1 ab—^cd 



(5,q±) = ^/ ^'kx T'^^^''5(s,kx)T^5-^^<^(s,qx -kx) (1.15) 



in momentum representation, or as the product 



2s 



(1.16) 



in the impact-parameter representation. 

The solution for the quark-hadron transition amplitudes T'^'^^'^'^ (^s, k^) and J"??-*-^^ (g^ kj_) 
at large invariant energy ^/s can be found using single Regge-pole parameterizations 
of the binary hadronic amplitudes ^'^°'^"^'^+'^" , ^JV^^^w j^nn^nn 

OLM(t) 



A^^--Ms,t) = iV^ ^--^j exp(i?^^i), ^^^^^^ 



Here aju (i), as (0 and an (t) are the dominant meson, baryon and diquark-antidiquark 
trajectories while Nm, Nm and No are normalization constants; mg = so and Roi 
is the interaction radius for the z-th trajectory. The following intercepts and slopes 
are found for the dominant Regge trajectories 

aM(0)~0.5, as {0)^-0.5, 0^(0)- -1.5 (1.18) 

and 

o'm (0) ~ a's (0) ~ a'^, (0) ~ 1.0 GeV'^ (1.19) 
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Using equations (fLT6]) and (fTTTl) the amplitudes T'?"^'''' (s, b^) and T'^i^^'^ (s, bj 
can be written as 



~- - I/O 1 / S \("m(0)+1)/2 . ,2 

r^^^^^fs, b^) = A^,r^ TT -A exp 

r^g-g^NN^ b^) = iV^ , , ( ) eXp ( ^ 

where Rm (s) and Ro (s) are the effective interaction radii given by 

RlAs) = RlM + <i{o)in (-^), 



:i.2o) 



Ulr 



Rl {s) = RId + oln (0) hi ^ 2 



(1.21) 



The substitution of the amplitudes (jl.20p into the factorization formula (11.16P gives: 

I^^^"^(s,b^) = 



, / S \ i("D(0)+°M(0)) 



AttRd (s) Rm (s) V 



exp 



1 1 

+ 



8Rl {s) 8Rl {s) 
(1.22) 



Consistency of Eqs. ()1.22p and ()1.17p requires the following relations between the 
Regge parameters and normalization constants ([HI [HU [58]): 



^ 1 1 1 



Rlis) Rl,{s) RlisY ^^23) 
2q (0)^ = an (0) + om (0) , 

{NMNn f' P , , , = NB-^y (1.24) 

Rd (s) Rm [s) R% (s) 

If only light n, d quarks are involved it can be assumed that ([Ml IZSl [84] ) 

(0) = a'B (0) = a'j, (0) = a' (0) , 
i?OM (0) = RIb (0) = i?L (0) = i?o (0) , (1.25) 
(NmNd)^^^ = Nb. 

Then the relations (|1.23p and (jl.24p can be fulfilled at all s. Otherwise, they can 
only be satisfied at sufficiently large s (Ref. 
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1.6 Deuteron Photo-disintegration in the QGSM 



It is convenient start the description of the photo-disintegration reaction consid- 
ering the amphtudes for spin-less constituents before going over to the case in which 
the photon spin is included. 



'wwwww\A«=: 




Figure 1.18: Diagrams describing three valence quark exchanges in the t-channel (left 
panel) and u-channel (right panel). 



From Eq. 11.161 the transition amplitude corresponding to each quark diagram of 
Fig. 11.181 can be written as: 

I^'^^P"(s,bx) = — f^'^^'?(5g)(s^b^) f«(55)^P"(s,b^), (1.26) 
where the amplitudes T'^'^^'?(5'?)(s, b^) and ^^(^^^-^^"(s, b_L) are given by (cf. (fr20D ) 



,1/2 1 / S \("m{0)+1)/2 / ^^2 



T^'^'^^"^\s,h^) = < Ltw=7^ ^ --2 exp 



(aB(0) + l)/2 ^ ^^2 

exp ' 



D{6q) 

(1.27) 



Here the effective interaction radii RM(&q) (s) and RD{eq) (s) are defined as 

^MiGq) (s) = ^0M(6g) + «M (0) ( ~ Z;2 ) ' 

\ ''''0 / 



:i.^ 



where -RoAf(6g) ^^'^ ^0D{6q) general, different from -Rqm ^'^'^ -^OD P-^IP - 

Since the photon possess an hadronic structure, a photo-reaction such as 'yd — )• pn 
could be assimilated to a hadron-hadron collision. In this case, the photon spin must 
be taken into account since its transversal polarization leads to a non-trivial angular 
dependence of the residue of the amplitude [59j . Spin effects in the QGSM where first 
introduced using the s-channel factorization property of amplitudes, developed in 
Ref. [84J with respect to the description of the electromagnetic nucleon form factors 
Fi and F2. 
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The resulting form factors for nucleons and pions comply well with the experi- 
mental data for both positive and negative values of (see Refs [851 ESI EZl EHl ESI 
[90]. 

In the case of the deuteron photo-disintegration, the amplitude can be written 
in the form: 

^ j d^k^ (Ap;A„|r'?(5''H?'"(s,k^)|A,;A(5,)) 



where pi, PSi and are the 4- momenta of the photon, deuteron, proton, and 
neutron, respectively, while A, is the s channel helicity of the i-th particle. Further- 
more, making the simplifying assumption that the spin of the (5g) state is 1/2, the 
amplitude ff'^^^^^i^ can be written as 



(Aq; A(5g) |r (s, k_L) I A^; A^) — 

^xM)^x. (^^^) ''^''h.M^^)) D^'^'^'''\s,k^) , (1.30) 

where and ex^ are the deuteron and photon polarization vectors, D"''^^'^^^'^\s , k_|_) 
is the scalar amplitude and niq is the quark mass. In analogy to qq — )• A^A^, which was 
analyzed in Ref. p4J , we can describe the spin structure of the amplitude T<ii^<i)^P^ 
in terms of eight invariant amplitudes 



(Ap;A„|f'?+(5^)-^*^" (s,kx)|A,;A(5,)) = 

Di(s,k_L) (^ApA^'^AnAfs,) +D2{s,k^) (fJy )aj, a/a„ A(5,) + 
1)3(5, k_L) '5apA,(0-s,)a„A(5,) +D4{s,k^) (CT^)Aj,A,(cra:)A„A(5,) + 
D5{s,k±) (0-y)ApA,(0-s/)A„A(5,) + Dfi{s,k±) (cr^)Aj,A,(fT^)A„A{5,) + 
1)7(5, k_L) (0-:^)ApA,(f^^)A„A(5,) +D8{s,k^) (f7^)Aj, A, (c^x)a„ A(5,) , (1-31) 

where the z- and x-axes are directed along the photon momentum and the momen- 
tum transfer k_|_, respectively, and the y-axis is orthogonal to the scattering plane. 
Now the experimental data on the proton form factor are in agreement with the as- 
sumption that the dominant contribution stems from the amplitude corresponding 
to the conservation of the s-channel helicities ([84j). Here the same assumption is 
used, taking into account only the amplitude Di{s,k±). Thus: 



(Ap; A„|T (5, p3_l) \Xd; A^) = 



,^7rf-^.(55)(>,^kx) I)i(5,p3x -kx). (1.32) 
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Furthermore, taking into account that at high energy ^ and finite momen- 
tum transfer i — |p3_Lp — sq the momentum k is almost transversal k = (fco, k_L, kz), 

where /cq — — O ( — — ) and J (P]s.± k_L(...) P3±, it is found the following rep- 

resentation for the spin structure of the 'yd — t- pn amplitude: 



(Ap; Xn\T {s, p3±) I Ad; A^) = 

uXp{P3)h^ (--4(s,t)p3_L + B{s,t)m) exjVx^{p4) 



:i.33) 



where 



A{s,t) 



P3. 



|P3. 



P A:2 



m„ 



Bis,t) 



I m,, 



8^2 



s m 



fc2 



(1.34) 



(1.35) 



and m is the nucleon mass. In the case of a Gaussian parametrization for (g^ k_|_) 

and Di (in Eqs. (24,25)) the ratio R = A{s,t)/B{s,t) is a smooth function of t. 
Further on, it will be assumed that it is a constant and will be considered as a free 
parameter. The differential cross section for the reaction jd — ?■ pn is then 



Jt 



[St \B{s,t)\^ + Su \B{s,u)\' 



64. ITS (p^™)2 
+(-1)^+1 2Stu Re{B{s,t)B{s,u)] 



(1.36) 



where / is the isospin of the reaction, i.e. 1 = 1 (or 0) for isovector (or isoscalar) 
photons. The kinematical functions St, Su, Stu in (11.361) are given by 

St = ^ J2 [^^7 (^(^3 -Pi) + m) ex^ {pA - m) 



X {R{P3 -pi) + m) et {ps + m) 



" ^ ^ [^^d (^(^3 -Pi) + m) ex-, {P4 - m) 



-^7 1 -^d 



s 



tu 



X ^x-, (RiPs - j5i) + m) e^^ (pa + m) 
^ ^ Tr [ex-, {RiP3 - Pi) + m) exAP^ - ' 



X el^ {R{P3 - Pi) + m) el^ (ps + m 
In order to fix the energy dependence of the amplitude B{s,t) it is required that 

2q]v(0)-2 

0CM=Q 



1.37) 



da 
dt 



(1.38) 
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Taking into account that 5j ~ s for s ^ sq it is found that 

/ „ \ Qiv(0)-l/2 

B{s,t)r^(^-j . (1.39) 

Moreover, a good approximation for the energy dependence of St{9''^^^ = 0) in the 
region = 1 7.5 GeV is 

St\ 0CM=o ~ (1.40) 

with C = (36 lb 3) GeV^. Using this approximation B{s,t) can be related to the 
Regge-pole exchange amphtude as 

\Bis,t)\^ = \MKegge{s,t)\^ , (1.41) 

where 

MRegge{s,t) = F{t) i—j exp 

Here aN{t) is the trajectory of the nucleon Regge pole and sq = 4 GeV^ ~ m^. 
1.6.1 Nonlinear Nucleon Regge Trajectories 



(1.42) 



According to the data on vrA'' backward scattering (see e.g. the review [73]) the 
nucleon Regge trajectory has a nonlinearity of the form: 

Qjv(t) = aiv(O) + a'^(O) t + ^ a%{0) Kit), (1.43) 

where aAr(O) = —0.5, a'j\f{0) = 0.9 GeV~^ are the intercept and slope of the Regge 
trajectory, and a%{0) = 0.20 0.25 GeV~^ is the coefficient of the nonlinear term. 
In ()1.43p it is introduced also a cut-off function K{t). Assuming that K{t) = 1 the 
amplitude will grow very fast with s at large t which would violate unitarity. To 
prevent this fast growth the exponential form 

K{t) = exp {-pt^) (1.44) 

is chosen with /3 = 0.008 GeV~^. The small value of /? does not destroy the param- 
eterization of a{t) for —t < 1.6 GeV^ derived from Ref. |73J. Note also that the 
phenomenological Regge trajectory (ll.43|) with a power- like or exponential cut-off is 
nonlinear only for moderate values of t; at large t the quadratic term becomes small 
and the trajectory becomes essentially linear again. 

On the other hand, the QCD motivated Regge trajectories as suggested by Brisu- 
dova, Burakovsky and Goldman (BBG) [83] show a different behavior at large t. As 
shown in the previous Section the screened quark-antiquark potential (Eq. II. lip 
leads to nonlinear meson Regge trajectories. These trajectories can be parametrized 
on the whole physical sheet as 

a(t) = a(0)+7[r"-(r-t)n (1.45) 



.an 



with < u < 1/2. The Hmiting cases u = 1/2 and — >• {■yu = const) correspond 
to the square-root trajectory 



a{t) = a{0) + 7 
and the logarithmic trajectory 



a{t) = a(0) - {jiy)ln 1 - - 



T 



:i.46) 



:i.47) 



respectively. Such trajectories arise not only for heavy quarkonia, but also for light- 
flavor hadrons. 

In order to find the possible forms of nonlinear Regge trajectories for mesons 
composed of light quarks, the analytical model shown in subsection 11.6.11 has to 
be considered. In this picture, the color flux tube is stretched between quark and 
antiquark at the tube ends. The varying string tension is introduced to simulate 
dynamical effects such as the weakening of the flux tube due to pair (qq) creation. 
Within this framework it is possible to recover the form of the underlying potential 
for a given Regge trajectory. Potentials leading to "square-root" and "logarithmic" 
Regge trajectories have been taken from Ref. [83] where it was demonstrated that 
their effect is very similar to what is obtained by a screened potential in the un- 
quenched lattice QCD. In addition, this approach is able to reproduce with very 
good accuracy all the available meson spectra. 

From experimental data it is known that the slopes of meson and baryon Regge 
trajectories are almost the same a'^ ~ ~ 0.9 — 1 GeV~^ (see Ref. ^66^)- Since 
the slope is determined by the string tension which, in turn, depends on the color 
charges at the string ends, a baryon Regge trajectory can be described using the 
form suggested by the analytical string model for the meson. In the baryon case one 
of the quarks at the ends is substituted by a diquark ( [91j ) . 



1.7 Predictions for the — > pn 
Differential Cross Section 

The dependence of the residue F(t) on t can be written in the form (see [751 EH])- 

1 



m = B 



t 



exp {Rft) + C exp (R^t) 



1.48) 



where the first term in the square brackets contains the nucleon pole and the sec- 
ond term accounts for the contribution of non-nucleonic degrees of freedom in the 
deuteron such as the D wave components and isobar contributions. 

The set of parameter for the Regge trajectory is the same used in the phenomeno- 
logical nonlinear trajectory |31j : 



B = 4.01 • lO-^kb^/2 . QgY^ (J ^ Q 7 GeV-2, 
Rl = 2 GeV-2, Rl = 0.03 GeV-^ a 



N{0) 



0.25 GeV" 
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the ratio R = A{s,t)/ B{s,t) being equal to 1. 

The QCD motivated logarithmic trajectory differs only for an overall normaliza- 
tion factor B taken as 1.8 • 10""' kb^/^- GeV and for the ratio R = A{s,t)/B{s,t) 
being equal to 2. 

The free parameter T is fixed by comparison with the available experimental 
data. The QGSM predictions, based on the QCD motivated "logarithmic" trajec- 
tory of Eq. 11.471 were shown in Fig. 11.61 for T =1.7 GeV^ taking into account the 
interference of the isovector and isoscalar components of the scattering amplitude. 
The forward-backward asymmetry arises from the interference of two amplitudes 
describing the contribution of isovector (p like) and isoscalar (a; like) photons, so 
that in this case the differential cross section can be written as 

T p+LJ 

=^ —L^ \St\BP{s,t) + B^is,t)\^+ 
Su \BP{s,u)-B^{s,u)\^ 

+2StuReiBP{s,t) + B'^{s,t)y {BP{s,u) - B'^{s,u))] , (1.50) 
having used the vector dominance model in the form: 

B'^{s,t) = BP{s,t)/V8, B'^{s,u) = BP{s,u)/V8. (1.51) 



1.8 Summary 

The study of the deuteron photo-disintegration at photon energies higher then 1 GeV 
can shed light on the reaction mechanisms at work in the transition from the nuclear 
to the QCD picture of the deuteron. 

To this aim, the world cross section data on the deuteron photo-disintegration 
have been reviewed, both at low (E^ < 1 GeV) and high (1 < < 5.5 GeV) in- 
cident photon energies. In the high energy region a new phenomenon emerges: the 
cross section scaling. According to a simple interpretation it could be a consequence 
of perturbative QCD. The scaling appears at different energies for different angles 
corresponding to a perpendicular momentum transfer of approximately 1 GeV^. 
From a closer look at the photo-disintegration data for the differential cross section 
emerges that the pQCD description is not plainly applicable especially at interme- 
diate momentum transfers. 

The complementary information derived from polarization observables suggests 
that non perturbative physics could be at work even at relatively high momentum 
transfers so that soft mechanism could mimic the scaling behavior. Some non per- 
turbative quark models well describe the experimental data. Among the various 
models, the Quark Gluon String model stands out since it gives a good description 
of both cross section and proton in-plane polarization transfers data. 

For this reason, the theoretical foundations and predictions of the QGSM have 
been fully illustrated. 

From the experimental point of view, the possibility to further check the the- 
oretical predictions and to discriminate among the different models requires the 
knowledge of the differential cross section in a wide range of photon energies and for 
very forward and backward angles. 
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Chapter 2 

Experimental Apparatus 



2.1 The Jefferson Lab Accelerator Facility 

The Continuous Electron Beam Accelerator (CEBAF) at JLab uses a high power 
electron beam with energies up to 6 GeV and 100% duty cycle. A schematic lay- 
out of the accelerator and experimental halls is shown in Fig. 12. 1[ The machine 
design, based on a race-track configuration, has the advantage of limiting the length 
and cost of the accelerating sections by recirculating the beam several times into the 
linacs while keeping low the energy loss due to bremsstrahlung in the curvature parts 
of the tracks. The injection system consists of a 45 MeV linac inserting electrons 




Hall B: The CLAS detector 



Figure 2.1: Schematic layout of the JLab Continuous Electron Beam Accelerator. 



into the main accelerator. The different beams can coexist in the linacs, but must 
occupy different arcs. At the heart of the machine are the 8 five-cell niobium super- 
conducting radio frequency (RF) cavities per cryo-module having an high average 
gradient (see Tab. 12. ip and operating at 1.5 GHz. Once into the main accelerator, 





Gradient 


Quality Factor 


Design specs 
Average performance 
Best performance 


5.0 MeV/m 
9.56 MeV/m 
18.0 MeV/m 


2.4 • 10^ 
7.0 • 10^ 
9.0 • 10^ 



Table 2.1: Performance figures for the CEBAF RF super-conducting cells. 

the beam enters the north linac and gains an additional boost of 600 MeV of kinetic 
energy and then is bent through 180° and sent to the south linac where it is boosted 
once again by 600 MeV. At this point, the beam could be sent to an experimental 
hall or recirculated through the north linac. The beam can be further recirculated 
up to 4 times, picking 1200 MeV each time is sent around the accelerator. 

The primary electron beam delivered by the accelerator can be separated to be 
sent to the experimental areas. Halls A, B, and C, for simultaneous experiments. 
The electron packet structure has a frequency of 1.5 GHz and can reach a bunch-to- 
bunch intensity modulated with a periodicity of 3, in order to deliver beam to the 
different Halls with variable intensity and energy so that beams in the three Halls 
can have either the same energy or energies multiples of 1/5 of the end-point energy. 

The experimental equipment in the Halls is complementary, addressing a wide 
range of physics issues. Hall A has two identical focusing high-resolution spectrom- 
eters with a maximum momentum of 4 GeV/c [9^. Hall C has two symmetric 
focusing spectrometers: one featuring acceptance of high-momentum particles, the 
other a short path length for the detection of decaying particles [94] . 

Hall B houses the CEBAF Large Acceptance Spectrometer, CLAS, designed for 
operation with both electron and tagged-photon beams and described in more detail 
in the next section. 

2.2 The CEBAF Hall B 

The Hall B facility is devoted to experimental studies of electro and photo- 
reactions requiring the detection of several only loosely correlated particles in the 
hadronic final state and measurements at limited (for the JLab standards) lumi- 
nosities. In order to achieve high detection efficiency for multi-particle final states a 
large-acceptance detection is required. In addition, a large-acceptance detector may 
compensate for various restrictions limiting the experimental operation to moderate 
luminosities. For example, in experiments using a tagged bremsstrahlung photon 
beam, a lower luminosity is required to keep down accidental coincidences. 

A view of the Hall B layout is shown in Figs 12. 3l and l2. 21 where the electron-photon 
beam line, the tagging spectrometer for photon beam operations, the large angle 
spectrometer CLAS and its associated equipments are represented. Also, solid-state 
polarized targets can only be operated at lower beam currents, corresponding to low 
luminosity. A high detection efficiency for multi-particle events and a useful event 
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rate at limited luminosity both require a detection system with a large acceptance. 
The instantaneous intensity of a tagged photon beam produced by the Hall B tagging 
spectrometer is limited to approximately 10^ tagged 7/s by accidental coincidences, 
typically resulting in relatively low background rates in a large-acceptance detector. 




The photon tagging spectrometer 



Figure 2.2: Left panel: Schematic view of the Hall B enclosure, including the CLAS metal 
supporting structures. The tagging spectrometer is shown at the entrance on the Hall, 
upstream whit respect to the electron beam direction. Right panel: The CLAS detector 
magnified view showing the schematics for the detector subsystems. 



2.2.1 Hall B Photon Production and Tagging 



Electrons from the CEBAF accelerator strike a thin target (the "radiator") just 
upstream from a magnetic spectrometer [95] (the "tagger" ). Photons produced in 
the radiator continue toward the CLAS target following the same beam-line through 
the magnet yoke that is traversed by the electron beam when the tagger is not in 
use. The bremsstrahlung tagging technique for direct measurement of the incident 
photon energy in photo-nuclear interactions is well established [96] . The JLab system 
is the first photon tagger in the multi GeV energy range to combine high resolution 
(~ 10~^So) with a broad tagging range (20% — 95%) of Eq. 

There are three modes of operation for photon beams in Hall B: a normal non- 
polarized mode, a circularly polarized mode, and a linearly polarized mode. These 
are governed by the polarization state of the incident electron beam, and by the 
nature of the photon radiator. 
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Figure 2.4: Hall B photon tagging system. Photons produced in the radiator continue 
toward the CLAS target following the same beam-line through the magnet yoke that is 
traversed by the electron beam when the tagger is not in use. The photon-tagging system 
can tag photons with energies between 20% and 95% of the incident electron energy, and is 
capable of operating with electron beam energies up to 6.1 GeV. 



2.2.2 Photon Beam Position and Profile 

The electron beam used in photon experiments is monitored up to a point just 
before the radiator by RF cavities providing a good determination of both the pro- 
duction point and initial direction of the photon beam. Then, downstream of the 
tagging magnet, the photon beam passes through a collimator. 

In the present experiment, the collimator has a 0.861 cm diameter hole in a 
cylindrical nickel block 25 cm long. This collimator is located approximately 14 m 
from the radiator, and restricts the diameter of the photon beam at the CLAS target 
to less than 3 cm. A magnet placed just downstream from the collimator sweeps 
aside low energy secondary charged particles created in the collimator. 



The photon beam position and size are monitored by a fixed array of crossed 
scintillator fibers located 20 m behind the CLAS target. The scintillators respond 
to the electron-positron pairs produced by photons in the CLAS target, in the at- 
mosphere between the target and the hodoscope, or within the scintillating fibers 
themselves. 

2.2.3 Photon Energy and Timing 

The photon-tagging system can tag photons with energies between 20% and 95% 
of the incident electron energy Eo, and is capable of operating with electron beam 
energies up to 6.1 GeV. The field setting of the tagger magnet is matched to the 
incident beam energy so that those electrons that do not radiate will follow a circular 
arc just inside the edge of the pole face, and will be directed into a secondary shielded 
beam dump below the floor of Hall B with a maximum capacity of 800 W. 

An electron that radiates a bremsstrahlung photon has lower momentum, and 
consequently smaller radius of curvature in the tagger dipole field, so that it emerges 
from the magnet along the open edge of the pole gap. A scintillator hodoscope along 
the fiat focal plane downstream from this straight edge detects this electron, and 
thereby allows for the determination of the energy and timing of the radiated photon. 

The overall geometry of this arrangement may be seen in Fig. 12.41 The tagger 
dipole is topologically a C-magnet with a full-energy radius of curvature of 11.80 m, 
a full-energy defiection angle of 30°, and a gap width of 5.7 cm. The required 
magnetic field in the gap is 1.13 T for a beam energy of 4 GeV. 

The focal plane hodoscope consists of two separate planes of scintillator detec- 
tors. The first detector plane (called the E-plane for energy) is composed of 384 
plastic scintillators 20-cm long and 4-mm thick. Their widths (along the dispersion 
direction) range from 6 to 18 mm in order to subtend approximately constant mo- 
mentum intervals of 0.003--Eo- Each counter optically overlaps its adjacent neighbors 
by one third of their respective widths, thus creating 767 separate photon energy 
bins that provide an energy resolution of 0.001 • Eq. 

The second detector plane (called T-plane for timing) lies 20 cm downstream of 
the E-plane and contains 61 counters, 2-cm thick, which are read out using PMTs 
attached by solid light guides at both ends (transverse to the particle direction) of 
each scintillator. The RMS timing resolution of these counters is 110 ps. 

A schematic of the electronics of the tagger focal plane can be seen in Fig. 12.51 
A scattered electron first encounters an E-counter and then passes a T-counter 
producing signals in both devices that are processed in separate electronic chains 
and participate in the construction of the tagger "event" . 

The signal coming from the E-counter single PMT is sent to an ADML board, 
amplified and discriminated to produce an ECL logic pulse. In turn, this pulse 
is then sent to an E-T mapping module which performs appropriate coincidences 
according to E and T geometries in order to reduce the production of accidental 
triggers. The same ECL pulse acts as a start for the E-counter TDC which stop will 
be given by the CLAS trigger. To produce a tagger event, an ECL signal from the hit 
T-counter must be in coincidence with the E-counter signal so that both signals can 
be interpreted in term of a definite electron trajectory. To this aim, the output of the 
E-T mapping module is elongated using a discriminator so that the discriminated 
T-counter signal determines the timing. Thus, a three-fold coincidence is required 




Figure 2.5: Tagger spectrometer electronics scheme. See the text for details. 



between the "left" and "right" signals from the two sides of a T-counter scintillator 
bar and a geometrically matched E-counter signal. For the timing measurement, the 
output of each T-counter discriminator is sent to one multi-hit FASTBUS TDC, in 
common start mode, with 50 ps per channel resolution. The signal for the TDC start 
comes from the CLAS trigger while the stop from the T-counter. Since the tagger 
has 61 T-counters, there are 61 possible pulses which are sent through a cascade of 
logic modules set in OR mode (referred as "Master OR", MOR). 

2.2.4 Photon Beam Flux 

The absolute photon flux is determined at very low flux rates by inserting a 
large lead-glass total absorption shower counter (TAG) into the photon beam. The 



Property 


Specs 


Glass composition 


55%Pb + 45%Si02 


Effective Z 


48.9 


Radiation Length 


2.36 cm 


Moliere Radius 


3.1 cm 


Refraction index 


1.67 


Gerenkov angle 


34.3° 


Density 


4.08 g/cm^ 


Effective length 


17 RL 


Energy resolution 


16% at =300 MeV 


Timing resolution 


4 ns 


Efficiency 


99% 



Table 2.2: Characteristics of the total absorption counter (from Ref. ^97)). 

TAG is essentially 100% efficient (see Tab. 12. 2p . allowing for the efficiency of the 
different elements of the tagging hodoscope to be determined. Due to counting pile- 
up problems, the TAG can only be operated at beam currents up to 100 pA, and 
must be retracted from the beam-line under normal running conditions. Thus, a 
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Figure 2.6: Hall B Total Absorption Counter (TAC). The TAG is essentially 100% efficient 
allowing for the efficiency of the different elements of the tagger hodoscope to be determined. 
Due to counting pile-up problems, the TAC can only be operated at beam currents up to 
100 pA, and must be retracted from the beam-line under normal running conditions. 



secondary monitor, linear in flux over a wide range, is cross-calibrated against the 
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TAG at low rates and then used to monitor the flux at higher intensities. More 
details on the photon flux normalization will be given in the next Chapter. 
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2.3 The CEBAF Large Angle Spectrometer (CLAS) 



The CLAS is nearly An spectrometer based on a toroidal magnetic field. The pri- 
mary requirements driving this choice were the ability to measure charged particles 
with good momentum resolution, provide geometrical coverage of charged particles 
to large angles in the laboratory, and keep a magnetic-field-free region around the 
target to allow the use of dynamically polarized targets. A view of the CLAS detec- 
tor and its ancillary equipment (such as the photon tagging spectrometer and the 
downstream photon flux normalization devices) is shown in Fig. 12. 7i A spherical co- 

Large-angle Calorimeter 



Drift Cliambers . Electromagnetic Calorimeter 




Figure 2.7: A schematic top view of the CLAS detector cut along the beam line. Typical 
photon, electron and proton tracks (from top to bottom) from an interaction in the target 
are shown (the photon track through the CLAS drift chamber system is drawn only to guide 
the eye). 

ordinate system is used in all of the descriptions. The z-axis is taken to lie along the 
beam direction, with B as the polar (scattering) angle, and (j) as the azimuthal angle. 
The X and y directions are then, respectively, horizontal and vertical in the plane 
normal to the beam. The CLAS magnetic field is generated by six superconducting 
coils arranged around the beam line to produce a field pointing primarily in the (/>- 
direction. The particle detection system consists of drift chambers (DC) [981 [991 [TOO] 
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to determine the trajectories of charged particles, gas Cerenkov counters (CC) [lOlj 
for electron identification, scintillation counters |102j for measuring time of flight 
(TOF), and electromagnetic calorimeters (EC) [1031 1104] to detect showering parti- 
cles (electrons and photons) and neutrons by induced ionization. 

The six segments are individually instrumented to form six independent magnetic 
spectrometers with a common target, trigger and data acquisition system (DAQ). A 
two-level trigger system is used to begin data conversion and readout. The Level-1 
trigger can make use of the fast information from the time-of-flight counters, the 
Cerenkov counters, and the electromagnetic calorimeters. A Level-2 trigger adds 
crude track finding using hit patterns in the drift chambers. The data acquisition 
system collects the digitized data and stores the information for off-line analysis. 

In the following sections, each of the CLAS subsystems and the Hall B ancillary 
equipment are briefly described. 

2.3.1 Target 

A variety of targets have been used to date, with dimensions adapted to the 
particular needs of either electron or photon running. The most common target used 
has been liquid H2. However, reactions have also been studied using liquid D2 (as in 
the present case), ^He, and ^He; solid ^^C, Al, Fe, Pb, and CH2; and polarized NH3 
and ND3 targets. All targets are positioned inside CLAS using support structures, 
inserted from the upstream end, and are independent of the detector itself (Fig. I2.8|) . 
Unpolarized targets have commonly been positioned inside CLAS at the center of 
curvature of the inner toroidal coil. 
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^ure 2.8: Unpolarized cryogenic target cooling system (from Ref. 



2.3.2 Start Counter 



The start counter (ST, see Ref. [105J) is a set of scintillators surrounding the 
target with the same geometric acceptance of the CLAS returning the time at which 
a charged particle has left the target region. This information, combined with time- 
of-flight information, determines the time needed to a charged particle to pass across 
the CLAS in tagged-bremsstrahlung experiments. The timing resolution relies on the 
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Figure 2.9: The start counter assembly is a set of scintillators surrounding the target with 
the same geometric acceptance of the CLAS returning the time at which a charged particle 
has left the target region. 

determination of the charged particle impact point along the start counter scintillator 
panels. For this reason, light readout is done at both ends of the panels. To save 
space, in the forward direction two sectors of the start counter panels are merged in 
one scintillator. This scintillator is bent in the forward polar angle such that both 
ends of the panel are in the upstream direction allowing for a placement of a single 
photomultiplier tube. The start counter timing resolution is of 260 ps. 

2.3.3 Torus Magnet 

The magnetic field for the momentum analysis of charged particles is generated 
by six superconducting coils arranged in a toroidal geometry around the electron 
beam line. There is no iron in the system, so the magnetic field is calculated directly 
from the current in the coils. 

The kidney-shape of the coils results in a high field integral for forward-going 
particles (typically high momentum), and a lower field integral for particles emitted 
at larger angles. At the same time, this coil geometry preserves a central field-free 
volume for the operation of a polarized target. 

At the maximum design current of 3860 A, the total number of amp-turns is 
5 • 10®. At this current the integral magnetic field ■ d I reaches 2.5 Tm in the 
forward direction, dropping to 0.6 Tm at a scattering angle of 90°. The magnetic 
field is reversible: when magnets are turned on in the normal condition, positive 
particles are bent outward from the beam axis. 

The main field component is in the (/(-direction, however, there are significant 
deviations from a pure (/>- field close to the coils. The effect of these deviations on the 
particle trajectories is minimized by the circular inner shape of the coil: particles 
coming from the target do not experience a significant defiection in (p when crossing 
the inner boundary of the coil. 
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Figure 2.10: Left panel: Contours of constant magnetic field for the main CLAS toroid in 
the mid-plane between two coils (see also Fig. I2.7|) . Right panel: Magnetic field vectors for 
the CLAS toroid transverse to the beam direction in a plane centered on the target. The 
six coils are shown in cross-section. 



Routine operation has been limited to 88% (3376 A) of the maximum current to 
keep internal mechanical stresses within conservative limits. 

2.3.4 Drift Chambers 

The CLAS toroidal magnetic field bends charged particles toward or away from 
the beam axis but leaves the azimuthal angle essentially unchanged. The magnet 
coils naturally define the detector into six independent tracking areas or "sectors". 
To simplify detector design and construction, 18 separate drift chambers were built 
and located at three radial positions in each of the six sectors. These radial locations 
are referred to as "Regions". The six "Region One" chambers surround the target in 
an area of low magnetic field, the six "Region Two" chambers are situated between 
the magnet coils in an area of high field near the point of maximum track sagitta, 
while the six " Region Three" chambers are located outside of the magnet coils [98^ 

m\\m)\. 




Figure 2.11: Left panel: Portion of a region three drift chamber, showing the layout of 
its two super-layers. The high-lighted drift cells have fired at a charged particle passing by. 
Right panel: Section of the whole region 3 drift chamber. 

To optimally fill the sector volume, the chamber bodies were designed to support 
wires stretched between two end-plates, each parallel to its neighboring coil plane. 
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and thus tilted at 60° with respect to each other. This design provides maximum 
sensitivity to the track curvature since the wire direction is approximately perpen- 
dicular to the bend plane. The wire midpoints are arranged in "layers" of concentric 
(partial) circles, with the wire positions shifted by half the nominal wire spacing in 
successive layers. This pattern of wires in neighboring layers, with a repeating pat- 
tern of two field- wire layers and one sense- wire layer, results in a quasi- hexagonal 
pattern with six field wires surrounding one sense wire. The cell size increases uni- 
formly with increasing radial distance from the target. 

For pattern recognition and tracking redundancy, the wire layers in each chamber 
are grouped into two "super-layers" of six wire layers each, one axial to the magnetic 
field, and the other tilted at a 6° stereo angle to provide azimuthal information. The 
stereo super-layer of Region one is an exception to this rule, consisting of only four 
wire layers due to space constraints. A detail of the wire layout is shown in Fig. 12.111 
The total number of sense wires in the drift chamber system is 35,148. 

A high voltage system maintains the sense wires at a positive potential and the 
field wires at a negative potential whose absolute value is half that of the sense 
wires. A layer of guard wires surrounds the perimeter of each super-layer, with the 
high voltage potential adjusted to approximate the electric field configuration of an 
infinite grid. This three- voltage scheme minimizes the effects of nearby grounded 
surfaces such as the end-plate. 

2.3.5 Cerenkov counters 

The Cerenkov Counter (CC) [lOlj performs the dual function of triggering on 
electrons and separating electrons from pions. The design of the Cerenkov detector 
aims at maximizing the solid-angle coverage in each of the six sectors out to an 
angle B = 45° with the least possible amount of material (to prevent degradation 
of the energy resolution). This is achieved by placing the light collecting cones 
and photomultiplier tubes (PMTs) in the regions of (p that are already obscured 
by the magnet coils, and covering as much of the available space as possible with 
mirrors. Since charged-particle trajectories lie approximately in planes of constant 
(/>, the placement of the PMTs in the shadows of the magnet coils does not affect 
the angular coverage. The light collection optics was designed to focus the light 
only in the (p direction, which preserves information on the electron polar angle. 
The Cerenkov radiator gas used in the detector is perfluorobutane (C4F10) having 
an index of refraction of 1.00153. This results in a high photon yield and a pion 
momentum threshold of 2.5 GeV/c. 

2.3.6 Time of Flight Scintillators 

The time-of-flight (TOF) counters [102j cover the polar angular range between 
8° and 142° and the entire active range in azimuthal angle (p. The scintillators are 
located radially outside the tracking system and the Cerenkov counters but in front 
of the calorimeters. Their alignment and relative positioning with respect to the 
beam-line are shown in Fig. 12.121 

The scintillator thickness of 5.08 cm is uniform throughout, chosen to give a large 
signal for traversing minimum-ionizing particles. Each scintillator is positioned such 
that it is perpendicular to the average local particle trajectory. The width of each 
counter subtends about 1.5° of scattering angle. The forward counters (at polar 
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upstream >^ downstream 

Figure 2.12: Coverage of one CLAS sector by the time-of-flight scintillators. 



A node 

Dynode 




scintillator bar 




Anode 
Dynode 



Passive patches 



TOF Pre-trigger Board 
(16 counters) 



Passive Splitter 



Delay 



Single-ended 
son Ribbon 



Gate 



ADC 
LeCroy 
ISSIM 



Counleis 1-16 
+ ^ 



Enable 



Disc 
LeCroy 

2313 



Coumeis 17-4S V 
"Triples" 



level 1 
trigger 



Start 



TDC 
LeCroy 
1872 A 



Figure 2.13: Scheme of the electronic setup for a TOF counter. Each counter consists of 
single scintillator bar with a PMT at each end. The timing resolution of the scintillators 
varies from 80 ps for the shorter counters to 160 ps for the longer ones. 
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angles 9 less than 45°) are 15-cm wide, and the large-angle counters are 22-cm wide, 
in the polar angle direction. 

Each TOF counter consists of a scintillator with a PMT at each end. The timing 
resolution of the scintillators varies from 80 ps for the shorter counters to 160 ps for 
the longer ones. The electronic readout of a single counter is shown in Fig. 12.131 

2.3.7 Electromagnetic Calorimeters 

The forward region of the CLAS (polar angles from 8° to 45°) is equipped with 
an electromagnetic calorimeter segmented in six independent modules of triangular 
shape. Each module is composed by 13 units, and each unit consists of 3 layers 
of lead sheets 2 mm thick interleaved with scintillator bars of 10 cm length and 
thickness (top part of Fig. [27141) . Such a configuration provides the best compromise 
between energy resolution and neutron detection efficiency. 

Each layer of scintillator is rotated by 120° with respect to the following layer to 
allow the stereo readout needed to resolve multiple hits into the calorimeter. The 
five units closest to the target are referred as the inner part of the calorimeter while 
the remaining units form the outer part. Inner and outer parts have independent 
electronic readout to sample the longitudinal shape of the showers developed within 
the calorimeter. The total thickness in terms of radiation lengths is 14.8 and the en- 
ergy resolution of the forward calorimeter is ge/E = lQ%^JE{GeV). The sampling 
fraction is approximately 0.3 for electrons of 3 GeV and greater, and for smaller 
energies, there is a monotonic decrease to about 0.25 for electrons of 0.5 GeV. The 
average RMS position resolution is 2.3 cm for electron showers with more than 
0.5 GeV of energy deposited in the scintillator. The timing resolution for electrons 
averages 200 ps over the entire detector. 

For two CLAS sectors the coverage is extended to polar angles of 75° by the Large 
Angle Calorimeter (LAC). The LAC is also a sampling calorimeter composed of two 
identical modules with a multi-layer structure of lead sheets and scintillator bars, 
similar to the forward calorimeter. A single LAC module consists of 33 layers, each 
composed of a 0.20 cm thick lead foil and 1.5 cm thick plastic scintillator bar |106j . 

The calorimeter geometry is projective, with scintillators of 10-cm average width. 
To avoid optical coupling, 0.2 mm thick Teflon sheets separate neighboring scintil- 
lators. Scintillator bars in consecutive layers are rotated by 90° to form a 40 x 24 
matrix of roughly 10 x 10 cm^ cells (bottom part of Fig. I2.14p . 

To improve the e~ /tt" discrimination, the LAC modules are longitudinally di- 
vided into inner (17 layers) and outer (16 layers) regions with individual light read- 
outs. 

Each module is equipped with 256 photo-multipliers tubes [107]. The gain stabil- 
ity of each LAC photomultiplier is monitored with a precision of higher than 0.5% us- 
ing a radioactive source of a-particles (^^^ Am) sealed within a YAP scintillator crys- 
tal |108j . Each LAC module corresponds to 12.9 radiation lengths and 1.0 hadronic 
absorption lengths and has an energy resolution of a(E)/E = 7.5%/y^ E{GeV). 
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Figure 2.14: Upper panel: schematic view of one module of the forward electromagnetic 
calorimeter. The expanded view shows the three layers of lead and scintillator (convention- 
ally referred as u, v, and w) making up the calorimeter sandwich structure. Lower panel: the 
large angle calorimeter tapered geometry (up- left), its lead-scintillator sandwich structure 
(up-right), and granularity (bottom-center). 
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2.4 Trigger System and Data Acquisition 



Many of the CLAS detector subsystems use the same type of electronic modules 
for readout to minimize the complexity of the system and to simplify maintenance 
and integration. Initially, all analog signals were digitized using commercial FAST- 
BUS modules. Recently, however, more electronic modules are available in VME, 
and are being integrated into the system. 

The wire signals from the most complex CLAS subsystem, the drift chambers, 
are read out using LeCroy 1877 pipeline time-to-digital converters (TDCs) operated 
in common stop mode [98j with a least count of 0.5 ns. 

For the PMT-based detectors, the signal times are digitized using LeCroy FAST- 
BUS 1872A TDCs, and the pulse heights are converted using LeCroy FASTBUS 
1881M analog-to-digital converters (ADCs). Most TDCs are fed from leading-edge 
LeCroy 2313 CAMAC discriminators set at a low threshold, typically 20 mV. The 
high-resolution TDCs, set to a nominal least count of 50 ps, are needed to achieve 
the required timing resolution for time-of-flight measurements. The ADC informa- 
tion is used to determine the energy deposition in each detector, as well as to correct 
to the TDC measurements for time walk. 

The inputs to the Level 1 trigger are generated using custom electronics to sum 
appropriate groupings of PMT signals that can be used to generate fast signals 
with high efficiency. In order to study the trigger in software, these inputs are also 
recorded in the data stream with LeCroy 1877 FASTBUS pipeline TDCs along with 
all other trigger inputs. 

2.4.1 Trigger System 

To acquire events of interest while minimizing the dead-time, a two-level hi- 
erarchical trigger system was designed |109j . The Level 1 trigger is dead-timeless, 
processing all prompt PMT signals through a pipelined memory lookup within 90 ns. 
The resulting signals are sent to a trigger supervisor module, where they are used 
to gate the front-end electronics [llOj . This custom board [lTl\ takes the Level 1 
(and Level 2) inputs from the trigger system and produces all common signals, busy 
gates, and resets required by the detector electronics. The Level 2 trigger uses the 
fast-clear capability by clearing events that satisfied Level 1, but which have no 
tracks in the drift chambers. 

2.4.2 Data Acquisition 

The CLAS data acquisition (DAQ) system was designed for an event rate of 
2 kHz. Continued development of the DAQ over five years resulted in routine op- 
eration at event rates between 3 and 4 kHz for the 2000/01 running period as can 
be seen from Fig. 12.151 The data from the various detector components are digi- 
tized in 24 FASTBUS and VME crates within the experimental hall and collected 
by the 24 VME Readout Controllers in these crates. These data arrays, or event 
fragments, are buffered and then transfered via fast Ethernet lines to the CLAS on- 
line acquisition computer in the control room. Three primary processes, the Event 
Builder (EB), Event Transport (ET), and Event Recorder (ER), comprise the main 
data-flow elements in the acquisition computer. The EB assembles the incoming 
fragments into complete events. The completed event is then labeled by a run and 
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Figure 2.15: Event rate and integrated events accumulated by tlie CLAS data system vs. 
time for a three year period. 
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event number, an event type, and the trigger bits configuration. At this stage the 
event has the final format it needs for off-hne analysis. The data fiow is shown in 
Fig. 12.161 the EB passes the completed events to shared memory |112j on the CLAS 
on-line computer. The ET system manages this shared memory, allowing access by 
various event producer and consumer processes on the same or remote processor 
systems. The ER picks up all events for permanent storage. 

Some events are transferred to remote ET systems, for raw data checks such as 
hit maps, status, and event displays, and for on-line reconstruction, analysis, and 
monitoring. The ER writes the data in a single stream to magnetic media. The 
output files are striped across an array of local RAID disks. 

The CLAS data flow requires communication among some 100 processes running 
on a system of processors distributed around the detector in the experimental hall 
and linked to two Symmetric Multi-Processor (SMP) machines and a number of 
workstations in the control room. The framework for this organization is the CEBAF 
On-line Data Acquisition (CODA) system |113t 1114] . The data runs are set up and 
performed by the CODA "Run Control" sequence, where a single run typically lasts 
for one or two hours. 




Figure 2.16: Data flow schematic for CLAS, with, the primary stream using the ETl shared 
memory. Also shown arc a variety of monitoring processes accessing data on other shared 
memory (ET2, ET3) managed by the ET system. See text for additional details. 
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2.5 Summary 



The Large Acceptance Spectrometer (CLAS) allows the study photo- and electro- 
induced nuclear and hadronic reactions by providing efficient detection of neutral and 
charged particles over a good fraction of the full solid angle. The detector capabilities 
(summarized in Tab. I2.3P are being used in a broad experimental program to study 
the structure and interactions of mesons, nucleons, and nuclei using polarized and 
unpolarized electron and photon beams and targets. 
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momentum (9 > 30°) 
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polar angle 


ag/6 1 mrad 




azimuthal angle 


a^/(f) 4 mrad 




time (charged particles) 


at w (100-250) ps 




photon energy 


ge/E ^M/Je 


Particle ID 


vr/K separation 


p< 2 GeV/c 




7r"'"/p separation 


p< 3.5 GeV/c 




TT~ misidentified as e~ 


< 10-3 


Luminosity 


electron beam 


L ~ 10^^ nucleon cm~^s~^ 




photon beam 


L ~ 5 • 10^^ nucleon cm~^s-^ 




data acquisition 


event rate 4 kHz 




data rate 


25 MB/s 


Polarized target 


magnetic field 


-^max — 5 T 



Table 2.3: Summary of the CLAS detector characteristics. 

The CLAS, in conjunction with the associated photon tagging spectrometer rep- 
resents an ideal experimental setup to measure the deuteron photo-disintegration 
differential cross section. In fact, it provides both large coverage in the proton scat- 
tering angle (from 6^^^ = 8° to 140°) and an high energy (up to 3 GeV) tagged 
photon beam. 
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Chapter 3 



Event Reconstruction 
and Data Analysis 

3.1 Introduction 

The 'yd — > pn differential cross section has been measured using the Hah B tagged 
photon beam in the energy range between 0.5 — 2.95 GeV and the CLAS large ac- 
ceptance spectrometer. The various CLAS production runs are organized by "run 
periods" where experiments sharing common running conditions take data simulta- 
neously. In this experiment (JLab E93-017, a part of the g2 run period experiments, 
where g2 stands for "7 on a deuteron target") a high-statistics dataset has been 
collected using a loose trigger requiring only one charged hadron in the final state 
in coincidence with a tagger counter. Four weeks of beam time have been used to 
accumulate approximately 2.4 billion valid triggers in two different periods during 
the August (18 days) and September (14 days) months of the year 1999. 

3.2 Running Conditions 

The E93-017 experiment is characterized as follows: 

• Energy of the incident electron beam: Eq=2.5 GeV ("August") and £^o=3.1 GeV 
("September"). The electron beam was impinging upon a gold radiator of 
10~^ radiation lengths; 

• Tagger spectrometer settings: full focal plane (0.2 — 0.95)£^o GeV, or prescal- 
ing OFF configuration for about 70% of the running time; 

high energy part of the focal plane (0.35 — 0.95)i?o GeV {prescaling ON con- 
figuration) for the remaining time; 

• Photon energy: 0.5-2.38 GeV ("August") and 0.6-2.95 GeV ("September"); 

• Current intensity: 10 — 13 nA; 

• Photon intensity on the single tagger counter: 2± 7- 10^ tagged 7/s for produc- 
tion runs (standard current intensity reported above) and about 10^ tagged 
7/s for normalization runs (low current intensity runs 2± 100 pA); 
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Figure 3.1: Top pictures: Schematic view of the experimental setup used during the g2 run 
period. From left to right: electron beam, radiator, the tagger spectrometer, the photon 
beam-line, the deuteron target (TGT), the CLAS detector, and the down-stream photon flux 
normalization components. For the present experiment, the hadronic events generated in 
the target volume are detected in CLAS using the time-of-fiight system, in coincidence with 
a tagged photon trigger and a start counter signal. Bottom pictures: the Level 1 trigger and 
the data acquisition system. Only the charged trigger scheme is indicated since the neutral 
trigger is not used in the current analysis. 



• Cryogenic target: mylar cylinder 4 cm in diameter and 10 cm long, filled with 
unpolarized liquid 1^2 at 4 K of temperature; 

• CLAS superconducting torus magnetic field: 3375 A (88% of the full field 
intensity) with positive particles out-bending (so called "normal" field config- 
uration) ; 

• Level 1 Trigger: coincidence between a signal in the tagger (Master OR) and 
either one charged track in the CLAS (80%) or 2 neutral hits (20%) in the 
electromagnetic calorimeters (the neutral trigger configuration is not used in 
the present analysis); 



CHARGED TRIGGER RUNS 
20086^20150 
20188^20350 
20584^20705 
20741^20760 



Table 3.1: Rims collected with the charged trigger used in the — ^ pn analysis. 

The experimental setup is schematically shown in Fig. 13. li Each data acquisition 
session is labeled by a run number so that, for the present experiment, run numbers 
go from 19,948 to 20,760. In addition to the basic dataset of 2.4 billion triggers. 
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additional 22 million empty target triggers have been taken mainly for background 
studies together with 178 million triggers at very low current intensity. These last 
data have been collected to measure the absolute photon flux and to study the 
photon tagger characteristics. In the present analysis, only the charged trigger runs 
have been used (shown in Tab. 13. ip for a total of 1,171 million of events. 

3.3 Event Reconstruction 



3.3.1 Physics Data 

The event files collected during the on-line data acquisition (CODA) sessions 
are stored on tape silo system, to be retrieved in the off-line phase. These files 
contain information on the detector ADC/TDC readouts at the time the event 
was digitized in CLAS and are called raw data files. The final codification, in 
the form of standard n-tuples, of the physical quantities of interest requires a multi 
step procedure performed using both official software packages (developed by the 
CLAS collaboration) and custom (user-written) code. The logic flow chart of this 
procedure is shown in Fig. 13.21 (right side): files collected during the experimental 
run are processed by the reconstruction code (called "AlC"). The first step of the 
reconstruction process is the conversion of raw ADC/TDC hits from the detector 
readout electronics into physical quantities such as energy and time. This task is 
accomplished using appropriate constants (also called calibration constants), defined 
for each detector readout subsystem and run period. The whole set of calibration 
constants is stored into the CLAS database system which can be accessed from the 
off-line software packages and provides the information used for the reconstruction 
in the form of custom files (called map files) or via an on-line SQL-database server. 

Using the calibration constants, the charged particle 4-momentum is determined 
from the time of flight and the track radius of curvature measured by the time-of- 
flight and drift chamber systems, respectively. Neutral particle momenta are basi- 
cally reconstructed using the information coming from the electromagnetic calorime- 
ters. 

The procedure for the determination of the calibration constants is iterative: a 
preliminary reconstruction pass is performed on a relatively small fraction of the to- 
tal collected raw files uniformly spaced through all the data taking period. During 
the preliminary pass, up-to-date pedestals constants for the various CLAS subsys- 
tems are used together with a "test" set of calibrations constants usually taken from 
a previous run period with similar experimental conditions. Then, the quality of 
the reconstructed data is monitored: if the result is not satisfying the constants are 
recalculated and successive passes of reconstruction are performed until the values of 
monitored quantities are acceptable. At this point the set of calibration constants is 
"frozen" into the database so that massive data reconstruction can be started (this 
process is referred as "cooking"). The output of the cooking process, is automati- 
cally screened by several data monitoring and filtering software utilities in order to 
check for possible problems or failures. For instance, diagnostic programs produce 
the status files containing the drift chambers occupancy and time-of-fiight paddles 
status for the run period of interest. 
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Figure 3.2: Block diagram of tlie off-line data processing. The raw data files (simulated 
-left panel- or measured -right panel-) are converted into reconstructed files and processed 
by CLAS applications or user written code to obtain standard n-tuple files. 
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In order to optimize the disk space requirements for the off-hne analysis, a fil- 
tering package is run on the entire reconstructed dataset so that a set of reduced 
files (skimmed files) are produced for each (group of) physics channel(s) of inter- 
est. These reduced files only contain events belonging the particular channel being 
analyzed, and only the banks strictly needed by the off-line analysis code. 

The physical quantities characterizing the reconstructed events are organized 
in a sequence of numerical fields colloquially called banks. Banks are accessed and 
managed using a FORTRAN-based memory allocation system known as BOS (Basic 
Object System) [115j . Bank objects are addressed by capitalized 4-letters pointers. 

3.3.2 Simulated Data 

A similar reconstruction chain is used for data files resulting from the GEANT 
Simulation (GSim) of the CLAS detector. As shown on the left side of Fig. 13.21 the 
output of the GSim tool is processed by the same reconstruction program used for 
real data files, with the exception that a calibration procedure is not needed. In this 
case, the database system provides the "reverse" calibration constants needed by the 
GSim tool to produce the appropriate ADC/TDC hits for the generated particles in 
their simulated interaction with the CLAS. The complete simulation of the CLAS 
detector is performed in two steps: GEANT simulation and Post Processing. 

In the first step, raw data files are created from the GSim tool considering the 
CLAS as an "ideal" detector where all the subsystems are working properly and the 
timing and spatial resolutions values are identical to the nominal ones. 

On the contrary, the post processing phase introduces finite effects for the timing 
and spatial resolutions and takes into account possible problems affecting the actual 
detector setup. In fact, during the run period of interest, some detector subsystems 
may have out-of-work electronic channels and the drift-chambers may also have out 
of work wires. 

An software package called GSim Post Processor (GPP) applies the necessary 
corrections to the GSim files using the information on the detector status (obtained 
from the diagnostic programs run after the reconstruction process) and takes into 
account the timing and spatial resolutions finite effects. 

The introduction of the finite resolution effects is done defining a certain number 
of smearing factors for the time-of-flight and drift chambers subsystems. 

The smearing factor / for the time-of-flight timing resolution is defined as follows: 

ai = 2fRi (3.1) 

where Ri is the ideal resolution of the TDC coupled to the time-of-flight scintillator 
paddle i. 

For what concerns the spatial resolution of the drift chambers, the physical 
mechanism exploited to convert the drift time into distance is shown in the left 
panel of Fig. 13. 3i A charged particle passing through a DC cell, releases a charge 
deposition on the sense wire that yields a drift time which is converted into a distance 
(DOCA, Distance Of Closest Approach) using pre-calculated functions. 

On the other hand, the overall track, as defined by the reconstruction software, 
takes into account the time-distance information coming from all the cells crossed 
by the particle by means of a best fit procedure. In this case, the distance of closest 




Figure 3.3: The tracking procedure in the CLAS drift chambers system. Left panel: a 
charged particle passing through a super-layer. The field wires are indicated by the open 
dots while crosses indicate the sensitive wires which collect the released charge. The drift 
times are associated to distances of closest approach (DOCAs) represented by the dotted 
lines. Right panel: two different values for the distance of closest approach of the particle 
from the sense wire can be defined, giving origin to a loss of spatial resolution. See text for 
details. 



approach from the sense wire under consideration will be close to the calculated 
distance, but not exactly identical (right panel of Fig [331) • 

Inside each DC region the spatial residual is defined by the difference: 

iDOCA^J-lDOCAgl (3.2) 

where j = 1, 2, 3 runs over the three DC regions. It turns out that spatial residual as 
a function of DOCAcaic is almost constant. The resulting distribution can be fitted 
using a 5*'*-order polynomial: 

S,=j2l^^B0CA^ (3.3) 

1=0 

so that: 

(Tj = const • fj • Sj (3-4) 

where the quantities /(i,2,3) £^re the three spatial resolution smearing factors. 

In order to introduce in the simulation the finite effects for the timing and spatial 
resolutions and the actual detector status, the GPP package is ran on the GSim 
output files before they are processed by the reconstruction program. 

3.4 Event Selection 

In order to constrain the two body deuteron photo-disintegration events, the 
missing mass spectra have been calculated from the relation 13.51 

Mi = (P^ + Pa- Ppf (3.5) 
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where P^, Pd, and Pp are the 4-momenta obtained by the reconstruction software 
for the identified particles. 

The resulting distributions exhibit a clear peak corresponding to the neutron rest 
mass. An example is shown in Figs 13.41 and Fig. 13.51 where the missing mass distri- 
butions are calculated for the intermediate incident photon energy oi =1.05 GeV 
and for proton scattering angles of 9^^^ = 25° and 9p^^ = 115°, respectively. The 
different plots represent the contribution to the missing mass spectra from the six 
CLAS sectors. At forward angles, the neutron peak is quite sharp, the typical sig- 
nal to noise ratio being of the order of 40 60 in this case. The shapes of the 
missing mass distributions are well reproduced using a Gaussian plus exponential 
fitting functions. The red/dotted arrows shown in the panels of Fig. 13.41 represent 
the (— 3c7, 3a) interval (where a is the standard deviation of the Gaussian fitting 
function) defining the width of the neutron peak. At backward proton scattering 
angles, the distributions peaks are less pronounced compared to those obtained at 
forward angles, since the momentum resolution of the CLAS is worse in the backward 
direction. Nevertheless, the neutron peak is still clearly evident. 

In order to exclude jd — )■ pX events yielding mass values too far from the neutron 
rest mass the width a of the Gaussian function is used to define {—3a, 3a) cuts 
around the distribution peak. 

These cuts are calculated, for each CLAS sector, as a function of the incident 
photon energy Ej and of the proton scattering angle 9p^. They have been used to 
define which proton events have to be associated to a two body photo-disintegration 
final state using also the information on the proton detection efficiency and detector 
fiducial cuts as explained in Sec. 13.81 An example of the behavior of the missing mass 
cuts as a function of 9p^ is shown in Fig. l3.6l for four incident photon energies around 
1 GeV: the different dotted lines correspond to cuts in the six CLAS sectors. The 
missing mass cuts become wider at large scattering angles reflecting the behavior of 
the missing mass distributions shown in Fig. 13.41 and Fig. 13. 5i 

At higher incident photon energies the neutron peaks in the missing mass distri- 
butions become less defined since the resolution of the spectrometer decreases. In 
fact, the reconstruction of the momentum of straighter tracks left by protons of high 
momentum is less accurate and moreover the number of events is smaller since the 
cross section decreases very rapidly as a function of the incident photon energy. 

An example of the behavior of the missing mass cuts as a function of the proton 
scattering angle 9p^ at incident photon energies around 2.9 GeV is shown in Fig. 13.71 
In this case, for the reasons explained above, the cuts are wider than those calculated 
for lower photon energies even if the missing mass distributions are still peaked 
enough for the fitting procedure to be applied. 
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Figure 3.4: Typical missing mass spectra calculated for the six CLAS sectors from Eg. 13.51 
for an incident photon energy of =1.05 GeV and a forward proton scattering angle 
9p^^ = 25°. The peak corresponding to the neutron rest mass can be clearly identified, 
and the shape of the distributions is well reproduced by a Gaussian plus an exponential 
fitting fimction. The red/dotted arrows shown in the panels represent the (— Scr, 3a) interval 
defining the width of the neutron peak, being a the Gaussian function standard deviation. 
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Figure 3.5: Example of missing mass spectra calculated for the six CLAS sectors from 
Eg. 13.51 for an incident photon energy of Ej =1.05 GeV and a backward proton scattering 
angle 9^^^ — 115°. The distribution is less pronounced compared to the one obtained 
at forward angles since the momentum resolution of the CLAS is worse in the backward 
direction. Nevertheless the neutron signal is clearly evident and the shape of the distributions 
is well reproduced by a Gaussian plus an exponential fitting function. Again, the red/dotted 
arrows represent the (— 3ct, 3ct) interval defining the width of the neutron peak, being cr the 
Gaussian function standard deviation. 
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Figure 3.6: Behavior of tlie missing mass cuts as a function of the proton scattering angle 
9™ for four incident photon energies around 1 GeV. The solid horizontal line represents 
the value of the neutron rest mass, while the different dotted lines represents the missing 
mass cuts for the six CLAS sectors 
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Figure 3.7: Behavior of the missing mass cuts as a function of the proton scattering angle 
6™ for four incident photon energies around 2.95 GeV. The solid horizontal line represents 
the value of the neutron rest mass, while the different dotted lines represents the missing 
mass cuts for the six CLAS sectors. In this case, the cuts are wider than those calculated for 
lower photon energies nevertheless the missing mass distributions are still peaked enough 
for the fitting procedure to be applied. 
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3.4.1 Interaction Vertex 

The interaction vertex of the incident photon with a deuteron in the target 
volume is calculated from events with more than one charged particle in the final 
state using the intersection between the different tracks. 

The distribution of the vertex z variable is shown in Fig. 13.81 where z is the 
direction of the incoming photon beam. As can be seen, the events are uniformly 
distributed along the 10 cm length target with a small contribution coming from a 
support structure located about 12 cm on the left of the target center. The solid 
line represents a fit to the distribution defining the effective target length. 
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Figure 3.8: Distribution of the vertex z-component for proton events originating in the 
target volume. The distribution shape reflects the target dimensions (10 cm length). The 
small peak on the left is the contribution of a target support structure and is eliminated by 
a cut on the vertex z-coordinate. The information on the interaction vertex and the x and 
y positions of the beam is contained in the MVRT bank in the case of multi-track events. 
The solid line represents a fit to the distribution defining the effective target length. 

In order to select events originating from the target volume containing liquid D2 
the following vertex cut is applied along the z-axis: 

- 5.3 < z < 4.7 (cm) , (3.6) 

where a displacement of —0.3 cm has been introduced to compensate for the shift 
of the center of the target from the origin of the z axis. 

3.4.2 Energy Loss 

The charged particles emitted in the target volume lose energy from ionization 
processes while traversing the liquid deuterium. The same effect occurs outside 
the target region when the charged hadron traverse the start counter scintillators. 
For these reasons a momentum correction has been introduced (Ref. [116j ) in the 
reconstruction procedure to compensate for the energy lost by the outgoing charged 
particles in the target and start-counter regions. 




3.5 Photon Flux Normalization 



In order to calculate the cross section, the incident photon flux at the CLAS 
target must be known and the goal of the normalization procedure |117j is to de- 
termine it as a function of the energy. This is done using the information of the 
out-of-time photons by measuring the number of accidental hits of each T-counter 
in a given time interval r, which is called the out of time window. Fig. 13.91 shows 
the histogram of the TDC timing spectrum for a single T-counter. The peak which 
caused the trigger and the accidental background can be clearly identified. Since the 
tagger TDCs operates in "common start" mode this peak implies that it was this 
T-counter which started all the other TDCs from the CLAS trigger associated to 
the hadronic event generated by this particular photon (electron) . Most of the time 
the TDCs which were not involved with the trigger will just overflow. But other 
times an accidental hit occurs in a T-counter which did not produce the trigger. 
This occurs with probability 1 — exp (— A''i*^(s) • r) where Ne\s) is the number of 
times the TDCs for each T-counter stop counting in the out of time window r. 




Figure 3.9: TDC timing spectrum of a single T-counter. The peak which caused the 
trigger and the accidental background can be clearly identified. Utilizing the accidental hits 
in the tagger, the normalization method determines the number of tagged photons from 
each T-counter. This is done by measuring the number of accidental hits of each T-counter 
in a given time interval r, which is called the out-of-time window. 



This measured electron rate is calculated "instantaneously" between each scaler 
dump in the DAQ, which is approximately every 10 sec. This procedure is based on 
the tagger T-counters and TAC TDCs information only (scaler-less). 

The total number of measured electrons (or out-of-time photons) giving the 
photon flux Ny{E!^'^) at a the photon energy E^^ is then determined from the sum: 

N,{E^^y) = ef^^^-Y,NP{s) (3.7) 

s 
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where: 

• i is the T-counter number (1 61); 

• Ne\s) is the number of measured electrons per T-counter per scaler interval 
corresponding to tagged photons not associated with an hadronic event in 
CLAS; 

• e?™™ is the tagging efficiency for the z*'^-tagger channel measured during a 
normalization run. 

The live-time of the DAQ in this time interval, is multiplied to the measured 
electron rate to give the number of measured electrons on each T-counter. So the 
number Ne\s) can be re-written as: 

AtW (s) = (s) • K{s) ■ LT(s) (3.8) 



where the DAQ live-time, indicated by LT(s), is calculated on an event basis as: 

Accepted Triggers 

LT s) = (3.9) 

iotal iriggers 



where 

^ TRGS|0].ClockHG(s) 

^ ^ lOKHz ^ ' 

is the gated clock for the run, and where i?^*) (s) is the electron rate for each T-counter 
which can be written as: 

= ./tdc(^) (3_ii) 



In Eq. [3TT] iVff)c(s) is the number of TDC hits per T-counter and per scaler 
time interval, -^t^qj^(s) is the total number of events with successful tagger recon- 
struction from the TAGR bank and N^p{s) is the number of times per scaler interval 
where there is a TDC hit before the out-of-time window r. 

In the present calculation scheme, -/V'^'^qj^(s) -r represents the total time, during 

10 s (scaler time interval), in which the TDC was open while N^p{s) ■ r represents 
a dead time. 

In order to calculate the photon flux, the tagging efficiency, which is defined as 
the number of tagged photons per energy interval divided by the total number of 
counts observed in the tagging counter defining that interval, must be measured. The 
knowledge of this quantity is very important because it directly affects the determi- 
nation of the photo-reaction cross sections. The tagging efficiency strongly depends 
on the beam setup. For instance, it will be reduced whenever tight collimation is 
required on the photon beam. 

The measurement of the tagging efficiency is preformed during low intensity 
"normalization" runs (schematically shown if Fig. I3.10p at approximately 10% of 
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the production beam current with the total absorption counter placed in the beam- 
line. This setup allows for a direct measurement of the photon flux from the total 
absorption counter (composed by a set of four lead-glass blocks with a detection 
efficiency close to 100%). Since the total absorption counter cannot handle photon 
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Figure 3.10: Schematic representation of a normalization run: the electron beam current 
is lowered to 100 pA and the total absorption counter is inserted into the beam-line, down- 
stream of the CLAS. The trigger is formed only by the tagger master OR (MOR). This 
setup allows for a direct measurement of the photon flux from the total absorption counter 
(composed by a set of four lead-glass blocks with a detection efficiency close to 100%). 



beam rates higher than 10^ Hz without noticeable degradation it could not be in- 
serted into the beam-line during the production runs to measure the photon flux 
directly. 

During the data taking, normalization runs were performed at the end of each 
data acquisition session characterized by a group of production runs with common 
experimental conditions, such as trigger configuration and tagger prescaling settings. 

Using this procedure, a total of 178 million normalization events have been col- 
lected. The relation giving the tagging efficiency for each tagger T-counter in a 
normalization run is: 

norm ^^'^^ [Tj • TAC] 



(3.12) 



norm 



where [T/^'"]^orm counts the number of coincidences Tieft • T^ght while [Tj • TAC] is 
the number of those hits detected in coincidence with the TAC. A small correction 
is done (~ 2 — 3%) to the tagging efficiency to account for the loss of photons from 
the target to the TAC [TT8] . 

The overall tagging efficiency as a function of the tagger spectrometer T-counter 
is shown in Fig. l3.11l for different runs. As can be seen also from Fig. 13.101 the higher 
numbered T-counters correspond to high energy electrons poorly deviated from the 
incident beam direction, so that their energy release AEq to the emitted photon is 
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small. According to the relation 0.y ~ nie/AEo the photon emission angle is larger 
and for this reason the collimation of the low energy photons is poor. Thus, the loss 
in efficiency is due the larger fraction of photons cut out by the collimator. On the 
contrary, the dip in the efficiency in the very first part of the tagging range is due 
to background coming from the electron beam. 
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Figure 3.11: Behavior of the tagging efficiency for the tagger spectrometer T-counters 
for different runs (indicated by the different symbols). The tagging ef&ciency is defined as 
the number of tagged photons per energy interval divided by the total number of counts 
observed in the tagging counter defining that interval. The knowledge of this quantity is very 
important because it directly affects the determination of the photo-reaction cross sections 
and strongly depends on the beam setup. 



A more detailed example of the stability of the tagging efficiency can be seen 
in Fig 13.121 where the variation of this quantity with respect to the mean value 
calculated for each run is shown for selected tagger counters. As it can be clearly 
seen, the tagging ratio is stable at level of ~ 3 4%. 

The normalization code have been used to process the normalization runs re- 
ported in Tab. 13. 2[ Since some T and E-counters were out-of-work, the procedure 
have been skipped for some of them, more precisely: T-counters: 20, 46, 59, 60 and 
E-counters: 66, 98, 101, 142, 214, 220, 222, 224, 239, 240, 285, 335, 337, 343, 380. 
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Figure 3.12: Variation of the tagging efficiency with respect to the mean value for some 
tagger spectrometer T-counters. As can be seen the variation is around 3-4%. 
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Table 3.2: Normalization runs acquired during the "August" and "September" run periods. 
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3.5.1 Validation of the Photon Flux Normalization 

Empty target runs have been used to check the photon flux normalization results. 
The collected statistics is not enough to perform a complete check over the whole 
incident photon energy range and proton scattering angles. Nevertheless the photon 
flux normalization has been validated comparing the distribution of the z component 
of the interaction vertex for proton events (normalized to the number of incident 
photons) between a full and empty target run. 

The distributions are shown in Fig. 13.131 where the empty target contribution 
(black line) exhibits two sharp peaks at z = —5 cm and z = h cm in addition 
to the one at z = —12.5 cm. These peaks are due to the interaction of the pho- 
ton beam with the target cell entrance and exit windows. The smaller peak at 
(— 15 < z < —10) cm is due to a structure holding the target cell into place. The 
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Figure 3.13: Distribution of the z-coniponent of the interaction vertex for the proton events 
(normalized to the number of incident photons) between a full (blue line) and empty target 
(black line) run. The difference between the two distributions (red line) is seen to be less than 
1% in the interval (—15 < z < —10) cm, confirming the validity of the photon normalization 
procedure. 

event distribution from the full target run, represented by the blue line of Fig. 13.131 
has contributions from events originating in the enclosure contains the liquid I?2- 
The difference between the full target and the empty target contributions is repre- 
sented by the red line of Fig. 13.131 In order to validate the photon flux normalization 
result the number of events originating in the interval (— 15 < z < —10) cm has been 
considered in the two cases. If the photon normalization procedure is correct these 
two numbers must be very close. In fact, the difference A'fuu — A^empty between the 
number of events in the chosen z interval is less then 1% since A^fuii = 10.86 it 0.93 
and iVempty = 10.94 it 0.63, thus confirming the validity of the normalization proce- 
dure. 
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3.6 Data Quality Check 



In order to select runs collected with stable running conditions, a given number 
of data quality checks have been performed on the reconstructed files used for the 
data analysis. The first quality checks have been performed during the "cooking" 
procedure where several run-based parameters, as the overall number of particles, 
the number of detected p, K^, and vr^ normalized to the incident photon flux in 
each CLAS sector, as well as the overall z vertex distribution (average value and 
width) are required to be constant at the few % level from run to run. 

An example of the stability relative to the normalized number of protons as a 
function of time is shown in Fig. 13.141 for the "August" run period characterized by 
the tagger spectrometer in the prescale ON configuration. For each run several files 
are analyzed and the average normalized proton yield is seen to be constant in time 
at less than 3% level. 
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Figure 3.14: Behavior of the normalized number of protons as a fmiction of time for the 
"August" run period characterized by the tagger spectrometer in the prescale ON configu- 
ration. Top plot: for each run several files are analyzed (represented by the single points) 
and the normalized proton yield is seen to be constant in time. Bottom Plot: Overall dis- 
tribution of the normalized proton yield fitted with a Gaussian function. The ratio between 
the standard deviation a and the mean /i is less than 3%. 

After this first selection, a second more accurate quality check based on the 
deuteron two body photo-disintegration proton yield (normalized to the number of 
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incident photons) is performed to monitor the data quahty for the present analysis. 
Significant file-based quantities have been considered: 

• the number of events with at least one charged particle in the final state for 
each tagger T-counter, Nl^ (T = 1 61); 

• the number of photo-disintegration events per 100 MeV bin, N"^^ {E = 1 -^25); 

• the number of photo-disintegration events per sector, A'^^ (-S* = 1-^6); 

where the index i runs over the number of collected data files N and each data file is 
considered a measurement of the above quantities. In order to define a criterion for 
rejecting the data files (measurements) which yield results are not compatible with 
the statistical fiuctuations of the sample, the Chauvenet's criterion [119] has been 
applied. This criterion is based on hypothesis that one result may be rejected if it 
shows an "unusually large" deviation from the average of the results obtained from 
the remaining N measurements. The meaning of unusually large is defined using 
a Gaussian distribution function to calculate the probability that such a result will 
occur: if this probability is larger than 1/A^, the result is not considered anomalous, 
since it can be due to a statical fiuctuation of the data sample. On the other 
hand, if the probability is much smaller than l/N , then it is very unlikely that the 
deviation of the result from the average is due to a fluctuation. In this case, the 
measurement is rejected and its unusually large deviation from the average of the 
remaining measurements is considered as originating from a uncontrolled change 
in the experimental conditions or as the result of mistake which was not recorded 
in the experimental logbook. The "rule of thumb" used in the literature rejects a 
measurement if the probability of obtaining it is found to be less than Pq = ^ . 

In order to apply the Chauvenet's criterion the mean values of N^, and A'^^ 
have been calculated according to the following relations: 

KA) = J, ^ (3.13) 

where = N^, or A^^ and a'^{A^) is the related statistical error. The corre- 
sponding errors on the mean values are: 

^^(^) = ' 1 (3-14) 



where the sum index i runs over the total number of collected data files A'^. 
For each of the quantities A^, the variable A can be defined as follows: 

A\A) = (3.15) 



Under the hypothesis that the variables A^ follow a Gaussian distribution also the 
variable A^{A) is distributed according to a Gaussian curve with vanishing mean 
value and standard deviation equal to the unity. 
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The Chauvenet's limiting probability Pq defines an interval of acceptability 
(— Aq, Aq) for A* such as: 




(3.16) 



-Ao 



SO that in the present analysis data files yielding A* values not satisfying the Chau- 
venet's criterion have been discarded. 

An example of the A distribution obtained for = (i.e. number of nor- 
malized photo-disintegration events) is shown in Fig. 13.151 for two different bins of 
incident photon energies since this quantity represents the most critical yield to be 
monitored in the present analysis. It can be seen that the distribution is practically 
flat and most of the data files yield values of A* whitin the limits of the fiducial 
interval (— Aq, Aq) represented by the two horizontal lines. 

Only two data files have been found to yield A values out of the confidence 
interval defined according to the criterion illustrated above and have been excluded 
from the current analysis. 

Another important parameter to be investigated is the uniformity of the response 
of the six CLAS sectors since they can be regarded as six independent spectrometers. 
The uniformity has been monitored considering the ratio between the normalized 
photo-disintegration protons yields from sectors 1 to 5 over those from sector 6, as 
a function of time (or run) . The choice of using CLAS sector 6 as a reference sector 
has been done since from the proton detection efficiency study, sector six has shown 
the best response (more details will be given in Sec. 13. Sp . 

As can be seen from Fig. 13.161 the distributions of the uniformity for sectors 1 to 
5 with respect to sector 6 are flat and no time dependence is evident. The values of 
the overall uniformity are reported in table [331 The lower value of the uniformity for 
sector 3 can be understood taking into account that a certain number time-of-flight 
paddles were out of work, as will be shown in detail in Sec. 13.81 



Sector 


Uniformity 


S1/S6 


0.97 


S2/S6 


0.97 


S3/S6 


0.91 


S4/S6 


0.99 


S5/S6 


1.01 


Average 


0.97 ±0.01 



Table 3.3: Uniformity of the response of the CLAS sectors. 

In conclusion, data files not satisfying the data quality checks have been discarded 
and so forth excluded from the current analysis. Following the procedures discussed 
above, about 7% of the originally collected data have been excluded. 
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Figure 3.15: Distribution of the A values calculated for the normalized photo-disintegration 
proton yield N'^^ at an incident photon energy of 0.95 - 1.05 GeV (top plots) and 2.15 - 
2.25 GeV (bottom plots). The horizontal lines represent the fiducial limits derived from the 
Chauvenet's criterion. The distribution is practically flat and most of the runs yield values 

of A* whitin the limits of the fiducial interval (— Ao, Aq) represented by the two horizontal 
lines. Only two files have been discarded from the overall dataset. 
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Figure 3.16: Ratio of the normalized photo-disintegration proton yields from sectors 1 to 
5 over sector 6 as a function of the run number. As it can be seen, all the distributions are 
flat and no time dependence on the run period is shown. The average values of uniformity 
are reported in table [ 
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3.6.1 Data Stability 



The data files processed according to the quality checks described above have 
been additionally checked to make sure that the experimental conditions in place 
during their acquisition where stable. 

Nevertheless, some experimental conditions must change since some runs have 
different requirements. For example, the maximum incident energy Eq of the elec- 
trons may change as well as the tagger prescale settings. 

For this reason, the data stability has been monitored using the total proton 
yield normalized to the incident photon flux and the rate of accidental coincidences, 
with different tagger prescale settings and different electron energies. 

The typical behavior of the normalized proton yield per tagger T-counter as a 
function of the incident photon energy is shown in Fig. 13.171 for different prescale 
settings of the tagger spectrometer and for different initial electron energies Eq. 
The overall normalized proton yield is proportional to the total photo-production 
cross section and is a decreasing function of the incident photon energy. On the 
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Figure 3.17: Normalized total proton yield per T-counter as a function of incident photon 
energy: the different symbols refer to different running conditions (tagger prescale ON or 
OFF and incident electron energy Eq). The behavior of the normalized proton yield per 
tagger T-counter does not show evident dependence on the experimental settings. 



contrary, it does not show evident dependence on the experimental settings since 
the contributions of runs performed in different conditions (identified by the symbols 
shown in the legend of Fig. I3.17P all have a similar behavior as a function of the 
incident photon energy. 

For what concerns the accidental coincidences, their contribution affects the data 



taking since the total trigger for photo-production experiments in CLAS is given by 
the coincidence signals coming from the CLAS and the tagger spectrometer (Master 
OR). 

The CLAS part of the trigger is formed when a charged particle, produced by 
an hadronic event in the target volume is detected in the near-by start counter 
scintillators, which seen from the target, cover the same solid angle as the time- 
of-flight system. The start counter signal has a time fluctuation of about 3 4 ns 
depending on production point of the particle in the target volume and impact 
position on the scintillators. 

At the end of its flight path through the CLAS, the charged particle fires the 
time-of-flight scintillators so that a valid CLAS Level 1 trigger may be formed. 
The overall time fluctuation of the time-of-flight signal depends on several factors. 
The velocity of the outgoing particle gives rise to a time fluctuation which can be 
estimated considering the time difference between a fast (high momentum) and slow 
(low momentum) pions, which is ~ 30 ns. Then, an additional ~ 15 ns fluctuation 
is introduced by the impact position of the particle in the TOF scintillator. In 
order to guarantee a suitable overlap with the start counter signal which is only 12 
ns wide the width of the TOF signal is chosen to be 100 ns wide. On the tagger 
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Figure 3.18: Distribution of the event time difference At between the tagger time and the 
start counter time for all outgoing charged particles (top plot), for proton events (middle 
plot), and for photo-disintegration events (bottom plot). 

spectrometer side, each of the 7 • 10^ 7/s flowing through the spectrometer can 
possibly produce a valid "good photon" trigger in one of the 121 T-counters so that 
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all the T-counters signals have to be aligned to produce a collective Master OR 
signal which is 15 ns wide. In this arrangement accidental coincidences are possible 
if, for example, multiple hits in the tagger have timings that can all be associated 
with a "good photon" pattern. 

Another possibility is realized when a photon hits the target but its energy 
is not measured in the spectrometer being out of its tagging range (very low or 
extremely high energy photons). Such photons may produce an hadronic event 
which is detected in CLAS while this event is accidentally in coincidence with a 
"good photon" tagger pattern which did not produce it. 

An example of the time difference distribution At between the tagger time and 
the start counter time is shown in Fig. 13.181 for all outgoing charged particles (top 
plot), for proton events (middle plot), and for photo-disintegration events (bottom 
plot). The ratio R of the number of events lying under time interval (—1,1) ns 
to the average number of events under the two near by peaks gives the accidental 
coincidences level, that is: 

l(iV(-3,-l) + Ar(l,3)) 

^ ^Fm^ (^-^^^ 

The values of ii as a function of the run number are shown in Fig. 13.191 It 
can be seen that the level of accidental coincidences has a relatively flat behavior 
as a function of the run number (time). On the contrary, two different levels are 
measured: one level is relative to the tagger prescale OFF configuration and has 
an average value around 1% while the other, relative to the tagger prescale ON 
configuration, is less then 2%. 
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Figure 3.19: Level of accidental coincidences R shown as a function of the run number. 

The relatively flat behavior has two different levels: one is relative to the tagger prcscale 
OFF configuration and has an average value around 1% while the other is relative to the 
tagger prescale ON configuration and is less then 2%. 



3.7 Momentum Corrections 



The missing mass distributions from jd — )• pX events have been calculated over 
the whole range of incident photon energies and proton scattering angles and the 
resulting peak position has been found to be off with respect to the reference value 
given by the neutron rest mass. This shift is not very large (of the order of 3%, at 
maximum) nevertheless the distortion introduced in the missing mass distribution 
shape can reduce the number of identified photo-disintegration events. 

A similar problem was found also in other CLAS experiments. For example, in 
the elastic ep scattering experiment the missing mass distribution peak for ep — >• eX 
events is shifted by the order of few % from the reference value of the proton rest 
mass |120j . This effect could be due to small misalignments of the drift chambers 
positions with respect to their nominal positions or to uncertainties in the magnetic 
field maps. Another possibility concerns the knowledge of the energy of the incoming 
photon which could be incorrect. Since it is not possible to disentangle between these 
two contributions, the shift of the missing mass peak is usually corrected introducing 
an empirical function f{Po, Oq, cpo) depending on the measured proton 3-momentum 
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Figure 3.20: Region of the kineniatie plane where the momentum corrections could be 
calculated using two body kinematics. 



P' = Po f{Po 



(3.18) 



where P' is the corrected momentum and and (j)Q are referred to measured proton 
angles in the LAB system. Assuming that angles and photon energy are correctly 
measured, the momentum correction function is factorized as follows: 



P' = Pofl{P",0o)f2iPo, 



(3.19) 



where P", fi, and /2 are functions to be determined from the data according to the 
following procedure: 
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1. first, the proton momentum range is divided in 200 MeV bins; 

2. in eacli bin of proton momentum, plioto-disintegration events are selected from 
the 7^ — )■ pX distributions. In order to reduce background, a {—2a, 2cr) cut 
around the peak position is applied, where a is the standard deviation obtained 
from a Gaussian fit of the missing mass distribution; 

3. event by event the momentum expected on the base of two body kinematics is 
calculated as a function of photon energy E-y and the proton scattering angle 
^0) which are assumed correctly measured. An example of the region of the 
kinematic plane where this procedure has been applied is shown in Fig. 13.201 

4. Each bin of proton momentum is additionally divided in 5° wide bins of the 
azimuthal angle (j)] 

5. in each (/)-bin the ratio ii^ = P^/Pq is calculated and its distribution is fitted 
using a Gaussian function to obtain the average value < i?^ >; 

6. the behavior of < > as a function of (j) is reproduced using a 2^^^ order 
polynomial (/2 in Eq. I3.19p . An example is shown in Fig. 13.211 for CLAS 
sector 1 and for three different values of proton momentum: Pq: 0.55 GeV 
(top plot), 0.65 GeV (middle plot), and 0.75 GeV (bottom plot). Then, a first 
correction is applied to the measured momentum Pq as a function of <p in tlie 
form P" = Po f2{PoAo) 

7. The procedure is repeated starting from point (4) and dividing each bin of 
proton momentum in 10° wide bins in the scattering angle 0. In this way, a 
second correction factor is calculated and applied to the measured momentum 
Pq as a function of 9 (/i in Eq. 13.190 . 

An example of the effect of the momentum correction procedure is shown in 
Fig. 13.221 for CLAS sector 3 (top plot) and sector 2 (bottom plot) for two different 
proton momenta and scattering angles: P=0.65 GeV, = 115° and P=1.15 GeV, 
9 = 15°, respectively. The missing mass distributions are sharper after the correction 
procedure while the peak position (indicated by Mx) is closer to the neutron rest 
mass value. 

In particular. Fig. 13.231 shows the overall situation before the application of the 
momentum correction procedure. As can be seen, the deviation from the unity of 
the ratio Mx/Mn, (being M„ the nominal value for the neutron rest mass) plotted 
against the proton scattering angle in the LAB system separately for the six CLAS 
sectors and for the intermediate proton momentum of 1.05 GeV, affects practically 
all sectors and shows a dependence on the proton scattering angle. The maximum 
deviation is of the order of 3% in some sectors (specifically sector 4 and 6). The 
same quantity, after introducing the momentum correction procedure is shown in 
Fig. I3.24[ It can be seen that the deviation from the unity is greatly reduced in all 
sectors together with the dependence of Mx/Mn on the proton scattering angle. 

The distributions of the Mx values obtained from the missing mass distributions 
for 'jd — )• pX events before and after the application of the momentum correction 
procedure are shown in Figs. 13.251 and 13.261 respectively. The top plot shows the 
overall Mx distribution for the six CLAS sectors, while the remaining six plots 
show the contributions of the different sectors. As shown in Fig. 13.251 before the 
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Figure 3.21: Behavior of the ratio < > as a function of </> in CLAS sector 1 for three 
values of the measured proton momentum Pq: 0.55 GeV (top plot), 0.65 GeV (middle plot), 
and 0.75 GcV (bottom plot). In each (/)-bin the ratio i?0 = P<^/Po is calculated and its 
distribution is fitted using a Gaussian function to obtain the average value < >. 
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Figure 3.22: missing mass distributions for jd — > pX events for sector 2 (top pfot) and 
sector 3 (bottom pfot) for two different proton momenta and scattering angles: P=0.65 GeV, 
= 115° and P=1.15 GeV, = 15°, respectively. The missing mass distributions are sharper 
after the application of the correction procedure and their peak positions are closer to the 
neutron rest mass value. 
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Figure 3.23: Before the momentum correction. Ratio Mx/Mn (being Mn the nominal 
value for the neutron rest mass) plotted against the proton scattering angle in the LAB 
system separately for the six CLAS sectors and for the intermediate proton momentum of 
1.05 GeV. The deviation from the unity affects practically all six CLAS sectors, and shows 
a dependence on the proton scattering angle. The maximum deviation is of the order of 3% 
in some sectors (specifically sector 4 and 6). 
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Figure 3.24: After the momentum correction. Ratio Mx/Mn (being M„ tlie nominal value 
for the neutron rest mass) plotted against the proton scattering angle in the LAB system 
separately for the six CLAS sectors and for the intermediate proton momentum 1.05 GeV. 
The situation is clearly improved compared to the one shown in Fig. 13.231 since the values 
obtained for Mx/Mn are of the order of 1 and the dependence on the proton scattering 
angle is removed by the correction. 



momentum correction is applied, the Mx distributions have large widths which are 
quite different in the six CLAS sectors. On the contrary, after the application of the 
momentum correction procedure, the widths of the Mx distributions are reduced 
and are uniform in the six sectors, as shown in Fig. 13.261 
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Figure 3.25: Before the momentum correction. Distributions of the Mx values obtained 
from the missing mass distributions for 'yd — >■ pX events. The top plot shows the overall Mx 
distribution for the six CLAS sectors, while the remaining six plots show the contributions 
of the different sectors. It is seen that, the Mx distributions exhibit large widths which are 
different in the six CLAS sectors. 



A summary of the improvement in the momentum reconstruction introduced 
by the correction procedure is shown in Fig. 13.271 The plot shows the position of 
the missing mass distribution peaks (associated with their RMS) before (solid/black 
squares) and after (solid/red dots) the correction while the horizontal line represents 
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Figure 3.26: After the momentum correction. Distributions of the Mx values obtained 
from the missing mass distributions for pX events. The top plot shows the overall Mx 

distribution for the six CLAS sectors, while the remaining six plots show the contributions 
of the different sectors. After the application of the momentum correction procedure, the 
widths of the Mx distributions are reduced and are more uniform in the six sectors. 
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Table 3.4: Numerical values of the missing mass distributions peaks and RMS before and 
after for the momentum correction procedure is applied in each sector, together with the 
standard deviations obtained fitting the corrected missing mass distributions. 
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Figure 3.27: Position of the missing mass distribution peaks (associated with their RMS) 
before (solid/black squares) and after the correction (solid/red dots). The horizontal line 
represents the nominal value of the neutron rest mass. 



For a more quantitative comparison, in Tab. 13.41 are given the numerical values 
of the missing mass distributions peak positions and RMS in each sector before and 
after for the correction together with the standard deviations obtained fitting the 
distributions after the correction procedure. 
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3.8 Proton Detection Efficiency 



3.8.1 Introduction 

The extraction of the photo-disintegration differential cross section requires the 
knowledge of the number of 'yd — t- pn events produced in the selected bins of the 
outgoing proton momentum and scattering angle. In order to count these events 
correctly the proton detection efficiency in CLAS have to be evaluated. 

The CLAS proton detection efficiency cannot be extracted from the photo- 
disintegration data because it also requires the detection of the photo-disintegration 
neutron in the CLAS. This procedure is not convenient since: 

• neutrons are detected by the electromagnetic calorimeters, having a much 
smaller angular coverage {Op^^ < 45° or 70° for two CLAS sectors) than the 
time-of- flight scintillators which are used to detect protons; 

• the neutron detection efficiency of the calorimeters is less then 50% (see Ref . |121j ) . 

For the above reasons, to evaluate the single proton detection efficiency a GEANT 
simulation of the CLAS detector has been used. 



3.8.2 GEANT Simulation 

In order to simulate — )• pn reaction in the CLAS, an appropriate number of 
events are produced using a simple Monte Carlo code. The events are generated 
uniformly in the CM system using an flat photon energy spectrum between 0.5 and 
2.95 GeV and then Lorentz boosted in the LAB frame. Since the proton detec- 
tion efficiency does not depend on the production mechanism, no cross section or 
background information is taken into account by the event generator. 

In order to have statistical errors smaller that those introduced from other quan- 
tities entering in the final cross section calculation, about 20-10^ photo-disintegration 
events have been generated. 



3.8.3 Proton Detection Efficiency Evaluation using GSim 

The proton detection efficiency calculated from the simulation is defined as the 
ratio between the number of generated events (N^J^p^) in a certain kinematic bin 
{APp^^, A9p^^, A(j)) over the number of reconstructed events (N^Ji^p^) in the same 
kinematic bin: 

_ N^^_^pJAP^^^,A9l;^^,A<P) 
^^^^^ - iV^NjKTOTA^p^ • 

A binning in proton momentum and polar scattering angle of AP^^^^ = 100 MeV 
and A6p^^ = 10° is chosen to match the one used in the final cross-section calcula- 
tion. On the other hand, a binning of A4> = 5° is chosen to better investigate the 
azimuthal behavior of the CLAS proton detection efficiency. In fact, according to 
the CLAS toroidal geometry, the reconstruction is worse on the borders of the six 
sectors, since these regions are adjacent to the superconducting magnet coils. 



The number of simulated photo-disintegration protons identified in CLAS (the 
numerator of Eq. I3.20p by the particle ID is sensitive to the quality of momentum 
reconstruction. In order to discard protons whose momentum is not correctly recon- 
structed, the differences between the angles of generated and reconstructed protons 



A6 = 0GEN(pLAB) _ ^REC(pLAB) ^ <A^™(pLAB) _ <^REC(pLAB) (g ^i) 

have been obtained and the standard deviations a/^g and cja^ of the distributions 
have been calculated from a Gaussian fit. 

The momentum dependence of UAe and o"a</> has been fitted using a polynomial 
function in the range (0.5 < Pp^^ < 3.2) GeV to obtain a variable cut used to dis- 
card reconstructed protons having polar and azimuthal angles differences outside 
the {—3(7, 3(t) interval, for 9 and 4> separately. 

The result of the GSim evaluation of the proton detection efficiency in CLAS is 
shown if Fig. 13.28] as a function of the CLAS azimuthal angle (/>, for ^^^"^^=0.95 GeV 
and 9p^^ = 45°. Each CLAS sector has a extension of 60° and it is sampled with 
a (p binning of 5° using the procedure described above. From Fig. 13.281 it is clear 
that the proton detection efficiency is almost constant (with an average value around 
95%) in the central regions of each CLAS sector while it drops vertically near the 
sector borders. 
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Figure 3.28: The CLAS proton detection efficiency evaluated using GSim as a function 
of the azimuthal angle (j) for GeV, 0^^^ = 45°, and with a binning of 5°. 

The shapes of the six CLAS sectors are reflected by the behavior of the proton detection 
eflticiency and are clearly identified by the sudden decrease of the efiiciency near the sector 
borders. 



The overall behavior of the proton detection efficiency as obtained from the GSim 
evaluation is shown in the contour plot of Fig. 13.291 as a function of the generated 
proton azimuthal and polar angles (running on the abscissa and on the ordinates, 
respectively) for a momentum Pp^^=0.95 GeV. The quota represent the proton 
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detection efficiency values and are shown in different colors/grays according to the 
scale indicated on the right side of the plot. The CLAS sectors are numbered in 
the plot from left to right and in each sector, a region of practically constant high 
efficiency (close to 95%) is clearly identified. The vertical low efficiency zones among 
sectors correspond to the decrease due to the presence of the torus magnet coils. The 
horizontal gaps in the efficiency function (along 6) correspond to out of work time- 
of-flight paddles or drift chambers wires introduced in the simulation by the GPP 
tool. 




Figure 3.29: Behavior of the proton detection efficiency obtained from GSim shown as 
a function of the azimuthal and polar angles (on the abscissa and on the ordinates, re- 
spectively) of the generated protons for Pj^^^=0.95 GeV. The quota represent the proton 
detection efficiency values and are shown in different colors/grays according to the scale 
indicated on the right side of the plot. The CLAS sectors are numbered from left to right. 
In each sector, a region of practically constant high efficiency (close to 95%) is clearly identi- 
fied. The vertical low efficiency zones among sectors correspond the torus magnet coils. The 
horizontal gaps in the efficiency function (along 9) correspond to to out of work time-of-flight 
paddles or drift chambers wires taken into account in the simulation using the GPP tool. 



3.8.4 Check of the GSim Evaluation Results 

The proton detection efficiency in CLAS can also be measured from the g2 data 
using the — t- pTT~p reaction. The definition of the efficiency in this case is different 
from the one given in Eq. l3.20l for the simulated events since the data does not convey 
any information on the "generated" particles. For this reason, the correct definition 
for the proton detection efficiency in CLAS in this case is the following: 

The quantity -^^^p^-x(p) ™ ^'i- 13-221 represents the number of protons identified 
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searching for a n^p pair from the particle ID and then calculating the event missing 
mass according to the relation X'^{p) = (P^ + — Pp — -Ptt-)^, where ^^p ,^- 
indicates the particles 4-momenta. 

On the contrary, ^^^^p^-p is the number of pir^p events identified directly by 
the particle ID. Clearly, in this last case, the particle ID identifies both protons, and 
there is no reason to count one proton as "missing" and the other as "identified" so 
the calculation is repeated exchanging the proton labels. 

The momenta of the identified particles in the final state have been compensated 
for the energy loss into the liquid D2 target and the start counter scintillators but 
no momentum corrections such those described in Sec. 13.71 have been applied in this 
case. 

The bins of the kinematic variables (AP^ A6'p'^^, A0) in Eg. 13.221 are always 
referred to the "missing" proton. The binning in proton momentum and scattering 
angle is identical to the one used in the simulation case. A different choice is done 
for the binning in the azimuthal angle (j) chosen to be /S.(f) = 10° since the proton 
detection efficiency in CLAS measured from the — >• pir'p reaction is mainly 
used to validate the simulation result and a full investigation of the CLAS proton 
detection capability in the (j) coordinate is not necessary in this case. 

Counting the number of protons identified by the event jd — t- p'7r~X (the first 
term in the denominator of Eq. I3.22p requires the knowledge of the event missing 
mass distributions in each bin of proton momentum and angles. These distributions 
have been calculated and fitted using a Gaussian curve and a straight line to iden- 
tify the peak position (corresponding to the proton rest mass) and the background 
contribution K (second term in the denominator of Eq. I3.22p as: 

J-3ct 

where the fit parameters p4, identify the straight line and a = p^ represents 
standard deviation of the Gaussian fitting function. In order to count the number of 
protons associated with the exclusive event ^d — )• pTr~X{p) a (— 3cr, 3a) cut around 
the peak of the missing mass distribution is applied. 

Examples of the obtained jd — )• pTr~X missing mass distributions are shown in 
Fig. 13.301 for an intermediate proton momentum Pp^^ = 0.95 GeV, an azimuthal 
angle 4> = 10°, and for three different proton scattering angles 9p^^: 10° (top plot), 
40° (middle plot), and 70° (bottom plot). The shape of the distributions is well 
reproduced by the Gaussian plus linear fitting function calculated in the interval 
(0.6 - 1.05) GeV. 

As in the reconstruction of the simulated proton tracks, the momentum deter- 
mination is worse at the border of the six CLAS sectors since these regions are 
closer to the torus magnet coils. In order to discard protons which momentum is 
not correctly reconstructed, the differences between the angles of "identified" and 
"missing" protons 

AO = ^^'^(PLAB) _ 0PID(pLAB) ^ ^ <^^(PLAB) _ 0PID(pLAB) (3.24) 

have been obtained and the standard deviations a^g and a^^ of the distributions 
have been calculated from a Gaussian fit. 

The momentum dependence of a^g and a^^ has been calculated using a poly- 
nomial fit in the range (0.4 < pLAB < 2.3) GeV to obtain a variable cut used to 
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Figure 3.30: Missing mass distribution for the event ^ pir X at an intermediate proton 
momentum P^^^ = 0.95 GeV, an azimuthal angle = 10° (corresponding to the center of 
sector 1), and for three different proton scattering angles 9^^^: 10° (top), 40° (middle), and 
70° (bottom). The curves superimposed to the plots are Gaussian plus straight line fitting 
functions calculated in the interval (0.6— 1.05) GeV. The background contamination is 
evaluated according to Eq. 13.231 
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discard reconstructed protons having polar and azimuthal angles differences outside 
the {—3a, 3a) interval, for 9 and (j) separately. 

An example of the proton detection efficiency in CLAS measured from the 
jd — )• pTT~p reaction is shown in Fig. 13.311 at the intermediate proton momentum 
pLAB^Q_95 ghAB ^ 450 rj.^^ black/fuh dots represent the GSim result for 

the proton detection efficiency while the red/open dots are the current result from 
the jd — )■ pTr~p data. As can be seen, the resulting efficiency values are in good 
agreement. It should be noted the different (p binning. The datum corresponding to 
a sector border azimuthal angle (pB obtained from 7c? — >• pTr~p data is roughly the av- 
erage of the GSim results obtained at (pB — 5° and </)b + 5°. The comparison between 
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Figure 3.31: CLAS proton detection efHciency in CLAS evaluated from GSIM (black/full 
dots) compared with the results measured from the jd — >■ p'iT~p reaction (red/open dots) for 
GeV and 9p^^ = 45°. The resulting values for the proton detection efficiency 
are in good agreement. Note the different (p binning. The datum corresponding to a sector 
border azimuthal angle (pB obtained from the jd — >■ p'K~p data is roughly the average of the 
GSim results at (pB — 5° and (pB +5°. 



the results obtained using the two methods cannot be extended to all the momentum 
bins covered by the simulated 70? — >• pn events since the phase space available to the 
jd —7- p-K^p reaction is smaller compared to the photo-disintegration case. This is 
clearly seen in Fig. 13.321 where the proton detection efficiency, measured from the 
jd — )• pTT~p reaction, is shown as a function of the proton azimuthal and polar an- 
gles (on the abscissa and on the ordinates, respectively) for Pp^^=0.95 GeV. The 
comparison of this result with Fig. 13. 29^ evidently shows that in this case the proton 
is detected only up to 9p^^ = 90° while in the photo-disintegration case, the angu- 
lar coverage extends up to 0p^^ = 125°, for this value of proton momentum. The 
quota represent the values of the proton detection efficiency in different colors/grays 
according to the scale indicated on the right side of the plot. The CLAS sectors are 
numbered from left to right. The regions of constant high efficiency (close to 95%) 



can be clearly identified. Note that the sampling rate A(j) = 10° yields a lower reso- 
lution with respect to the GSim result (Fig. 13.291) . The vertical low efficiency regions 
among sectors are due to the presence of the torus magnet coils. The horizontal gaps 
in the efficiency function (along 9) correspond to out of work time-of-flight paddles 
and drift chambers wires. 
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Figure 3.32: Proton detection efficiency in CLAS measured from the —>■ p-K~p reaction 
shown as a function of the proton azimuthal and polar angles for Pj^^^=0.95 GeV. The 
quota represent the values of the proton detection efficiency in different colors/grays ac- 
cording to the scale indicated on the right side of the plot. The CLAS sectors are again 
numbered from left to right. The regions of constant high efficiency (close to 95%) can be 
clearly identified, even if the sampling rate A0 — 10° corresponds to a lower resolution. The 
vertical low efficiency regions among sectors are due to the presence the CLAS torus magnet 
coils. The horizontal gaps in the efficiency function (along 9) correspond to to out of work 
time-of-flight paddles and drift-chamber wires. 

The result on the proton detection efficiency in CLAS measured from the — )• pTVp 
data is used, in the present analysis, to validate the GSim result in the appropriate 
kinematic range. In fact, the limitation in phase space coverage does not allow a full 
comparison between the final results nor the usage of the proton detection efficiency 
measured from the jd — )• p7r~p data in the calculation the final photo-disintegration 
cross section. 

The information obtained from the jd — >■ p'K~p data is used to define the final 
timing and spatial resolutions smearing factors to be introduced in the simulation 
in order to match the characteristics of the simulated detector with those of the real 
apparatus. The comparison is done iteratively, running partial reconstructions on 
the simulated files until the required resolution smearing factors are determined. 

After this procedure, spatial and timing resolutions smearing factors equal to 
1 have been applied to the simulation results for all proton scattering angles. An 
exception is done for very forward proton scattering angles 8° < 6^^^ < 10°, where 
a spatial resolution smearing factor equal to 2 has been used. 
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A final comparison between the results obtained from the GSim evaluation of 
the proton detection efficiency using the — t- pn reaction (black/full dots) and 
the measurement from the data using the ^ pn'p reaction (red/empty dots) 
is shown Fig. 13.331 separately for the six CLAS sectors. The points represent the 
weighted averages of the proton detection efficiency in CLAS over proton momenta 
in the interval (0.5 — 1.0) GeV for the sectors central region (defined by the interval 
A(p = ±10° assuming the sector azimuthal between (—30°, 30°)) and are shown as a 
function of the proton scattering angle 9p^^. It is clearly seen how the GSim result 
is well matched by the ^d — )• pTr~p result. In addition, the actual detector status is 
well reproduced by the simulations: this is shown by the dips in the proton detection 
efficiency found for the same scattering angle 6p^^ in both results (see, for instance, 
sector 3 for 9^^^ = 75° and 95°). 

In order to evaluate the uniformity between the two results, the distribution for 
the values of ratio R = ''^^'^^ (defined only for the kinematic bins AP}'^^ and AOl;^^ 
populated by both kind of data) between the efficiency calculated with the two 
techniques is shown in Fig. 13.341 for proton momenta in the range (0.5 — 1.0) GeV. 
As can be seen from the legend, the mean value of the distribution is of the order of 
the unity while the spread is about 8%. The best agreement between the simulation 
and data is found for sector 6. This finding can be used to proper normalize cross 
section data from other sectors. 
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Figure 3.33: Final comparison between the results obtained from the GSim evaluation of 
the proton detection efficiency using the jd — > pn reaction (black/full dots) and the mea- 
surement from the data using the -fd — )■ pir^p reaction (red/empty dots). The points repre- 
sent the weighted averages of the detection efficiency over proton momenta in the interval 
(0.5 — 1.0) GeV for the sectors central region (defined by the interval A(f> = ±10° assum- 
ing the sector azimuthal angle between (—30°, 30°)) and are shown as a function of the 
proton scattering angle 6p^^. It is clearly seen that the GSim result is well matched by 
the 7d — > p'K~p result. In addition, the actual detector status is well reproduced by the 
simulations: this is shown by the dips in the proton detection efficiency found at the same 
scattering angle 9l^^^ in both results. 




Figure 3.34: Distribution of the ratio R — ^aAm between the proton detection efhciency 

" «GSIM ^ •' 

measured from — >■ p'K~p and simulated from — >■ pn (defined only in the kinematic bins 
APp^^ and A^^^^ populated by both kind of data) for the sectors central region (defined 
by the interval A0 = ±10° assuming the sector azimuthal angle between (—30°, 30°)). The 
parameters representing the Gaussian fit are shown in the legend. It can be seen that the 
mean value of the distribution is of the order of the unity while the spread is about 8%. 
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3.8.5 Fiducial Cuts and Detector Acceptance 



The information obtained from the simulated 'yd — )• pn reaction in CLAS allows 
the determination of the so-called "fiducial" regions in the azimuthal angle (j). In 
fact, a fiducial region can be defined, for each sector S and each bin of proton 
momentum AP^^^ and scattering angle A^^"^^, as the section of the CLAS 
detector characterized by a constant efficiency. 

The procedure used to identify the fiducial regions requires the calculation of a 
reference value for the proton detection efficiency (indicated as eo) in each CLAS 
sector and for each bin of proton momentum and scattering angle. This reference 
value is obtained averaging the contributions of the most central regions of each 
sector. The corresponding indetermination associated to eo is indicated by 6eo. 

Then, in each CLAS sector (for each bin of proton momentum and scattering 
angle) the proton detection efficiency, ej, in the remaining i^^ cp-hin on the left and 
on the right of the central region is investigated and the values of whose error 
bars 6ei do not overlap with the band defined by the relation cq it Seo are rejected. 

An example of this procedure is shown in Fig. 13.351 for sector 2, for a proton 
momentum Pp^^=1.15 GeV and the scattering angles dp^^ = 15° (top plot) and 
(9LAB ^ ggo (bottom plot). 

The horizontal shaded area corresponds to the allowed variation band eg ± (5eo 
while the vertical lines mark the boundaries used to exclude efficiency data consid- 
ered out from the allowed variation band according to the criterion outlined above. 
It can be seen that the procedure retains only the contributions from to the (p region 
where the efficiency behavior is constant. 

An example of the resulting fiducial cuts in the CLAS sector 2 for Pp^^=1.15 GeV 
is shown in the contour plot of Fig. 13.361 The boundaries of the fiducial (p regions 
(along the abscissa) are shown for each proton scattering angle (along the ordinates). 
The mean values of the efficiency within the region are indicated by the color / gray 
code. The dip found at 9p^^ = 55° is due to an out of work time-of-flight paddle, 
as shown also in Fig. 13.331 

The procedure outlined above is applied in each CLAS sector and for each bin 
of proton momentum and scattering angle in order to derive the full set of fiducial 
cuts. This is needed to properly weight the number of photo-disintegration events 
entering in the final differential cross section calculation since it will be expressed 
as a function of the outgoing proton momentum, scattering angle, and sector. 

For this reason, the (p information contained in the original proton detection 
efficiency ecsiM must be averaged over the CLAS sectors fiducial regions so to obtain 
the sector-based efficiency-acceptance function: 

where 1/wi = (^ecsiM {^Pp^^, A9p^^, A0j))^ is the weight assigned to each datum 
from the original efficiency function, S stands for the CLAS sector, and the index 
i runs over the (;/)-bins of sector S for the given proton momentum and scattering 
angle. 

The result of this calculation for the efficiency-acceptance function is shown in 
the contour plot of Fig. 13.371 for the six CLAS sectors. The abscissa report the 
proton momentum while the ordinates run on the proton scattering angle (both 
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Figure 3.35: An example of fiducial cuts in CLAS sector 2 for P^^ =i.l5 GeV. Upper 
panel: for 6^^^ = 15°. Lower panel: for 9^^^ = 65°. The horizontal shaded area shows the 
allowed variation band eo ± Seo while the vertical lines delimit the fiducial region according 
to the behavior of the proton detection efficiency as a function of the azimuthal angles (j). 
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Figure 3.36: Fiducial cuts in the CLAS sector 2 for GeV. The boundaries of 

the fiducial (j> regions (along the abscissa) in the sector are shown for each proton scattering 
angle (along the ordinates). The mean values of the efficiency within the region are indicated 
by the color/gray code shown in the legend. The dip found at 6^^^ = 55° is due to an out 
of work time-of-flight paddle, as shown also in Fig. 13.331 



considered in the LAB system). Quotas identify regions of constant values of the 
efficiency-acceptance functions (numeric values can be obtained from the color /gray 
scale shown on the left). As can be seen from Fig. 13.371 the behavior of the six CLAS 
sectors is pretty uniform, except for the holes introduced by out of work channels in 
some sectors. 
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Figure 3.37: Result of the calculation of the efhciency-acceptance function for the six 
CLAS sectors. The abscissa report the proton momentum while the ordinates run on the 
proton scattering angle (both considered in the LAB system). Quotas identify regions of 
constant values of the efficiency-acceptance functions (numeric values can be obtained from 
the color/gray scale shown on the left). As can be seen from Fig. 13.371 the behavior of the 
six CLAS sectors is pretty uniform, except for the holes introduced by out of work channels 
in some sectors. 
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3.9 Background Subtraction 



The determination of the number of "true" photo-disintegration events requires 
the knowledge of the contamination from background events. In fact, the missing 
mass distribution from — )■ pX events usually exhibits a peak in correspondence 
to the neutron rest mass value sitting over a continuous background mainly due to 
photo-production events originating from the target cell walls. 

In the present analysis, the background contribution is evaluated from the data. 
Only the events detected in the fiducial regions of the CLAS, defined using the 
fiducial cuts obtained in the previous section, are taken into account. For each 
event, the missing mass is calculated and distributions are obtained for each bin of 
proton scattering angle 0p ^ and incident photon energy E-^. 

Examples of the obtained missing mass distributions are shown in Fig. 13.381 for 
the six CLAS sectors, for the intermediate energy £J^=0.95 GeV, and a forward scat- 
tering angle = 15°: the peak corresponding to the neutron rest mass is clearly 
identified and each distribution is fitted with a Gaussian function plus exponential 
curve to take into account the continuous background. 

At larger angles (for the same energy S^=0.95 GeV), the distribution are wider 
and the neutron signal is less pronounced, as shown in Fig. 13.391 Nevertheless, they 
are still well described by a Gaussian plus exponential function fit. 

The contamination from background events under the missing mass distribution 
peak is calculated according to the relation: 



r3- eP4+P5M^^ 

J —ia 

*^peak 



^^^'^7) = '-^^ — (3.26) 



where a is the standard deviation of the Gaussian function, and are the pa- 
rameters of the exponential function, and N^^,^'^ is the number of events under the 
(— 3cr, 3(t) interval around the peak position. 

The behavior of the background contamination A; as a function of the proton 
scattering angle is shown in Fig. 13.401 for the six CLAS sectors and for the interme- 
diate incident photon energy i?^=0.95 GeV. The contribution is around 5% and is 
relatively flat for scattering angles less than 90°. At larger angles the contamina- 
tion is higher since the missing mass distributions get wider since the momentum 
resolution of the detector decreases. 

A similar situation occurs at high incident energies (more than 2.1 GeV) where 
the background subtraction cannot be directly performed since the statistics is not 
enough to obtain a prominent neutron peak in the missing mass distributions. To 
overcome this limitation, the behavior of the background contamination has been 
studied as a function of the incoming photon energy and extrapolated up to 3 GeV 
using a linear fit. An example of the result is shown in Fig. 13.411 for the six CLAS 
sectors and for = 35°. The last three data points, corresponding to the energy 
interval 1.8 < < 2.1 GeV, are obtained averaging the contributions coming from 
all sectors since the statistics is not enough to have a clear neutron signal in the 
missing mass spectra, and are identical in the plots shown. It is clearly seen that 
the data are well reproduced using a linear fit. 

At large proton scattering angles the background contamination k is larger but 
it can still be extrapolated using a linear fit of the available points. This can be 
seen from Fig. 13.421 for the six CLAS sectors and for = 135°. In analogy to 
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Figure 3.38: Missing mass distributions for the six CLAS sectors, for the intermediate 
energy £'^=0.95 GeV, and a forward scattering angle = 15°: the peak corresponding to 
the neutron rest mass is clearly identified. Each distribution is fitted with a Gaussian func- 
tion plus exponential curve to take into account the continuous background. The red/dotted 
arrows identify a (— Scr, Scr) interval around the peak position, being a the standard deviation 
of the Gaussian function. 
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Figure 3.39: Missing mass distributions for the six CLAS sectors for 9^ = 125° and for an 
energy incident photon energy £^.^=0.95 GeV. The distribution are wider and the neutron 
signal is less pronounced since at backward scattering angles the resolution of the CLAS 
spectrometer is lower. Nevertheless, the missing mass distributions are still well described 
by a Gaussian plus exponential function fit. The red/dotted arrows identify a (— Str, Str) 
interval around the peak position, being a the standard deviation of the Gaussian function. 
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Figure 3.40: The behavior of the background contamination fc as a func-tion of the proton 
scattering angle for the six CLAS sectors and for the intermediate incident photon energy 
Ej=0.95 GeV. The contribution is around 4% and is relatively flat for scattering angles less 
than 90°. 
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Figure 3.41: Behavior of the background contamination A; as a function of the incoming 
photon energy for the six CLAS sectors and for = 35°. The last three data points, 
corresponding to the energy interval 1.8 < < 2.1 GeV, are obtained averaging the 
contributions from all sectors and are identical in the plots shown. It is clearly seen that 
the data are well reproduced using a linear fit. 



Ill 



the situation described in Fig. 13.41] the last three points are obtained averaging the 
contribution of all sectors. 





Figure 3.42: Behavior of the background contamination fc as a function of the incoming 
photon energy for the six CLAS sectors and for 6™ = 135°. The last three data points, 
corresponding to the energy interval 1.8 < < 2.1 GeV, are obtained averaging the 
contributions from all sectors and are identical in the plots shown. At this proton scattering 
angles the background contribution is larger but it can still be extrapolated using a linear 
fit of the available points. 
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3.10 Evaluation of the Systematic Errors 



Systematic errors may result from faulty calibration of the experimental equip- 
ment or from biased procedures used in the analysis of the data. These errors must 
be estimated from reconsidering the status of the experimental conditions and the 
techniques used to analyze the data. 

In the present analysis, the main contributions to the overall systematic uncer- 
tainty come from the following procedures: 

• determination of the number of incident photons; 

• background subtraction; 

• evaluation of the proton detection efficiency. 

The evaluation of the systematic error on the photon flux normalization, is done 
using the information on the tagging efficiency as described in Sec. 13.51 As a conse- 
quence, it can be seen that the systematic uncertainties are estimated to be around 
3 - 4%. 

The background subtraction procedure introduces a systematic error of the order 
of 10% when repeated with different cuts while the proton detection efficiency in 
CLAS has been evaluated using two different techniques (as shown in Sec. 13. 8p . 

The comparison between these results, in the full overlap of the kinematics of 
the two reactions introduces a systematic error of the order of 5%. 

The different contributions are summarized in Tab. 13.51 



Procedure 


systematic error 


Photon flux normalization 

Background subtraction 
Proton detection efficiency 


3-4% 
8-10% 
5-6% 


Overall 


10-12% 



Table 3.5: Summary of the main contributions to the overall systematic uncertainty. 

It should be underlined anyhow that a more accurate evaluation of the systematic 
uncertainties is still underway. 
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3.11 Summary 



In this Chapter the various phases of the data analysis chain are fully described. 
First, the CLAS event format and the reconstruction procedure are introduced since 
these structures have been used to access and process the g2 data to perform all 
subsequent analysis. 

Then, the data calibration procedure is briefly described followed by the defini- 
tion of the photo-disintegration event. At this point, the photon flux normalization 
procedure has to be introduced since its needed to normalize the number of photo 
disintegration events used to calculate the final cross-section. 

In order to obtain accurate results, data quality checks are performed using 
appropriate criteria and several corrections are introduced to improve the quality of 
the reconstruction. 

The CLAS proton detection efficiency is evaluated both from a GEANT sim- 
ulation of the detector and from the data itself in order to verify the consistency 
between the results. After this step, the fiducial regions of the detector are defined. 

Then, the background contamination is evaluated from the data and the sub- 
traction procedure is described. 

Finally, an estimation of the overall systematic error has been given. 
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Chapter 4 

Results 



4.1 Introduction 

This Chapter is devoted to the illustration of the final results on the deuteron 
photo-disintegration cross section obtained from the analysis of the full statistics of 
the data collected with the CLAS detector during the g2 data production period. 

The CLAS data represent the first wide ranging survey of the deuteron photo- 
disintegration cross section in the few GeV region (0.55—3.0 GeV) of incident photon 
energy. The statistical uncertainties are of the order of 2 5% below 1.5 GeV, and 
of the order of 10% for aroud 2.5 GeV. 

The CLAS experimental results are discussed and compared to the published 
data. They give for the first time the possibility to study the photo-disintegration 
cross section as a function of the incident photon energy over the full range of 
proton scattering angles, thus expanding the view limited to four scattering angles 
= 36°^ 52°, 69°, and 89°) investigated in the past years (as discussed in the 
first Chapter of the present thesis) . 

Then, the angular distributions of the cross sections, integrated over the CLAS 
azimuthal coordinate, are shown for incident photon energies from 0.5 to 3.0 GeV 
in 100 MeV bins. These data provide for the first time, the cross section at very for- 
ward {9p^ = 15°) and backward proton scattering angles {9p^ = 155°) allowing to 
test the available theoretical predictions in regions where the cross section behavior 
becomes steeper. 

In addition, both the new CLAS results and the published angular distributions 
are compared to the non perturbative calculation of the Quark Gluon String Model 
(QGSM) since, as shown in the first Chapter of this thesis, this model describes 
rather well the data. 

4.2 The Differential Cross Section 

The photo-disintegration cross section is calculated weighting the number of 
protons (Njd^pn) detected in the CLAS from photo-disintegration events in each 
bin of proton scattering angle and incident photon energy, by the effective number of 
deuteron targets, the absolute number of incident photons A'^(l — -F), where F is the 
photon beam attenuation through the experimental target, the efficiency- acceptance 
function r?GSiM (Sec. 13. 8p and taking into account the background contamination k 



(as defined in Eq. I3.26j) . i.e. 



da 




pxNA{l-F)N^An r/GSiM 



)) 



(4.1) 



where A is the target molecular weight, Na is the Avogadro's number, and p and x 
(as defined in Sec. 13.4. ip the target density and efi'ective length, respectively. 

The CLAS cross section data are fully reported in Sec. 14.41 (at the end of the 
present Chapter) tabulated as a function of the incident photon energy (Tab. 14. ip . 
For each energy bin, the values of the differential cross section as da/dfl, obtained 
summing the individual contributions of the six CLAS sectors, are reported as a func- 
tion of the proton scattering angle in the CM system, together with their statistical 
uncertainties. The tables report also the values of da/dt as a function of the proton 
scattering angle in the CM system, with their associated statistical uncertainties. 

4.2.1 Uniformity of the Six CLAS Sectors 

As stated in Chap. 2, the CLAS is composed of six independent sectors that 
can be considered as six similar spectrometers. This allows for the investigation 
of the uniformity of the final cross section results since the measurements are per- 
formed independently in each spectrometer (sector), and provides an estimate of the 
systematic uncertainty due to the CLAS detector efficiency. 

The uniformity of the results has been evaluated considering the differences be- 
tween the values of the cross section obtained in each bin of incident photon energy 
and proton scattering angle for each sector S (with 5 = 1-^6) and the related 
average over the six sectors, divided by the average itself, or: 



The resulting distributions for the values of V are shown in Fig. 14.11 separately for 
the six CLAS sectors (numbered from left to right and top to bottom). In each 
panel the statistical parameters of the distributions are shown and it can be seen 
that a systematic shift of the order of 2 5% is present in all sectors. 
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Figure 4.1: Distributions of tlic differences between the value of the cross section obtained 
in the individual sector and its average over the six sectors, divided by the average itself, as 
defined in Eq. 14.21 for each bin of proton scattering angle and incident photon energy. The 
result is shown separately for the six CLAS sectors (numbered from left to right and top to 
bottom). In each panel the statistical parameters of the distributions are shown and it can 
be seen that a systematic shift of the order of 2 5% is present in all sectors. 
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4.2.2 The Energy Dependence of the Cross Section 

The new CLAS data give for the first time the opportunity to study the energy 
dependence of the deuteron photo disintegration differential cross section over the 
wide range of proton scattering angles covered with the CLAS detector and with a 
very good statistical accuracy. 

The results on the differential cross section ^ are presented in Fig. 14.21 and 
Fig. 14.31 (solid/red dots) as a function of the total energy s (here t, u, and s are 
the momentum transfer to the proton, to the neutron and the total energy in the 
Mandelstam notation) for proton scattering angles in the range 15° < < 155°. 
As can be seen from Fig. 14.21 and Fig. I4.3l the cross section decrease as a function of 
s is extremely steep and this behavior is enhanced for scattering angles around 90°. 
The blue/solid line superimposed to the experimental results represents a power law 
fitting function proportional to s"". 

The resulting exponents n obtained from the fits of the cross section for the 
various angles are plotted in Fig. 14.41 It can be seen that the cross section decrease 
as a function of s occurs according to different exponents n for the different scattering 
angles. At very forward and backward angles, corresponding to the dominance of 
the t or u channel, the exponent is of the order of —9 and —7 respectively. At more 
intermediate angles (around 90°), where both the t-dominance and n-dominance are 
suppressed, n drops around —13. Moreover, one can also argue that an asymmetry 
exist between forward (n = —9) and backward angles {n = —7) and a local maximum 
can be identified for 9^^ = 155°. 
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Figure 4.2: The CLAS results on the differential cross section ^ as a function of the total 
energy s (where t and s are the momentum transfer to the proton and the total energy in 
the Mandelstam notation) for proton scattering angles in the range 15° < < 85°. As 
can be seen, the cross section decrease as a function of s is extremely steep and this behavior 
is enhanced as the scattering angle approaches 90°. The blue/solid line superimposed to the 
experimental results represents a power law fitting function proportional to s". 
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Figure 4.3: The CLAS results on the differential cross section ^ as a function of the total 
energy s (where t and s are the momentum transfer to the proton and the total energy 
in the Mandelstam notation) for proton scattering angles in the range 95° < < 155°. 
It can be seen that the cross section tendency to a fast decrease as a function of s, slows 
down going toward backward scattering angles. The blue/solid line superimposed to the 
experimental results represents a power law fitting function proportional to s". 
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Figure 4.4: Exponents n obtained by tlie angular survey on the cross section. The data 
behavior is of the type ^ ~ (const)s" as a function of the total energy s. It can be seen 
that the cross section decrease as a function of s follows different exponents n for different 
scattering angles. At very forward and backward angles (corresponding to the dominance of 
the t OT u channel, and where u is the momentum transfer to the neutron) the exponent is of 
the order of —9 and —7 respectively. At more intermediate angles (around 90°), where both 
the ^-dominance and u-dominance are suppressed, n drops around —13. A clear asymmetry 
exist between forward {n = —9) and backward angles {n — —7). The errors on the data 
points are those given by MINUIT on the fitting parameter n . 
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4.2.3 Angular Distributions 

The angular distributions derived using Eq. 14.11 are shown in Figs. 14.51 Fig. 14. 6| 
and Fig. 14.71 for incident photon energy bins of 100 MeV, in the range from 0.55 up 
to 2.9 GeV. As seen from the figures, the data cover a broad range in the proton 
scattering angle: 15° < < 155° allowing for the first time to investigate the 
behavior of the cross section at very forward and backward scattering angles. 

Results from published data are also included in the plots. More specifically, 
results in the highest energy bins obtained by the Mainz/DAPHNE [13] experiment 
overlap the CLAS data in the photon energy region 0.55 < < 0.75 GeV. At higher 
photon energies, results from the SLAG [171118) experiments are shown together with 
the recent data from the JLab Hall G [ 24l 125) and Hall A |26) collaborations. As can 
be seen from Fig. 14.51 in the energy region below 1 GeV, the angular distributions 
still exhibit a trace of the excitation of baryon resonances in the second resonance 
region, which comprises Pii(1440), L>i3(1520), and S'ii(1535). At = 0.55 GeV, 
the Mainz data are slightly higher then the GLAS data, mostly at intermediate 
scattering angles. This discrepancy is greatly reduced at forward and backward 
angles and also becomes less evident at higher photon energies where the two datasets 
are well compatible to each other within the total statistical uncertainties. 

At E^ = 0.65 GeV, both the angular dependence and the absolute value of the 
differential cross section of the GLAS data are in good agreement with the Mainz 
result. Starting from E^ = 0.75 GeV the comparison can be extended also to the 
SLAG and JLab Hall G data. It turns our that, in general, the GLAS results well 
agree with data from these experiments. 

For incident photon energies higher than 1 GeV the angular dependence of 
the photo-disintegration cross section changes, and a more pronounced asymme- 
try emerges between the forward and backward angles, as can be seen from the 
bottom plots of Fig. 14.51 

The angular distributions obtained in the energy range 1.35 < E^ < 2.85 GeV 
are shown in Figs [33] and [4Jl At = 1.6 GeV and E^ = 1.9 GeV the GLAS 
results are compared with the angular distributions measured by the JLab Hall A 
collaboration which cover a few scattering angles in the range 26° < < 143°. 
The GLAS results confirms the Hall A data and extends that result at the very 
forward and backward angular regions where the cross section increases. In general, 
the clear angular asymmetry shown by the GLAS data persists in the whole range 
of explored photon energies, as is confirmed by Fig. 14.71 which shows the results up 
to E^ = 2.85 GeV. 

At photon energies higher that 1 GeV, where the resonance contributions starts 
to fade, the predictions of the QGSM [31] are included in the plots. As shown in 
Figs 14. 5[ 14. 6[ and Fig. 14. 7[ the angular distributions are well reproduced by the 
QGSM calculation (shown as a solid line). This model accounts very well for the 
persistent forward to backward asymmetry seen in the data by invoking the interfer- 
ence of the iso- vector and iso-scalar components of the photon, which is constructive 
at forward angles and destructive in the backward direction. 
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Figure 4.5: Angular distributions of the deuteron photo-disintegration cross section mea- 
sured by the CLAS (solid/red dots) in the incident photon energy range 0.55 — 1.25 GeV. 
It can be seen that the CLAS data cover a broad range in the proton scattering angle: 
15° < < 155°. Results from Mainz [13] (open dots), SLAC [171111] (open squares), JLab 
Hall A [26j (open diamonds) and Hall C [531 Hi] (open triangles) are included and are gener- 
ally in agreement with the CLAS data. The angular dependence of the photo-disintegration 
cross section shows a pronounced asymmetry between the forward and backward proton 
scattering angles. For incident photon energies higher than 1 GeV, the solid line represents 
the non perturbative calculation of the QGSM [31] which well reproduces the data. In par- 
ticular, the QGSM calculation accounts for the forward to backward asymmetry shown by 
the data. 
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Figure 4.6: Angular distributions of the deuteron photo-disintegration cross section mea- 
sured by the CLAS (solid/red dots) in the incident photon energy range 1.35 — 2.05 GeV. 
Results from already published experiments such as SLAG [171 HI] (open squares), JLab 
Hall A [26] (open diamonds) and Hall C [24] [25] (open triangles) are included, where avail- 
able. The cross section asymmetry between the forward and backward angles persists in this 
photon energy range. The QGSM [ST calculation (shown as a solid line) well reproduces 
the data behavior also in this energy range. 
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Figure 4.7: Angular distributions of the deuteron photo-disintegration cross section mea- 
sured by the CLAS (sohd/red dots) in the incident photon energy range 2.15 — 2.85 GeV. 
Results from SLAC [171 HI] (open squares), JLab Hall A [26l (open diamonds) and Hall 
C [24l [25] (open triangles) are included, where available. The angular distributions are well 
reproduced by the QGSM pT calculation (shown as a solid line) which accounts for the 
data asymmetry between forward and backward scattering directions. 



4.2.4 The Cross Section Scaling 



As fully discussed in the opening Chapter of this thesis, the dependence of the 
differential cross section as a function of the incident photon energy has been often 
investigated in the past years to evidence the scaling behavior. 

Such an effect is believed to show up in the cross section data when the proton 
transverse momentum 

p2 = ii?^MBsin2(0CM) ^ (4 3) 

being Mn the deuteron mass, is of the order of 1 (GeV/c)^. In this case, the pQCD 
assumption that soft and hard contributions factorize may be realized. 

In the case of deuteron photo-disintegration process, the cross section is expected 
to scale with a unique simple power law of the square of the total energy s 

f = / (^r) (4-4) 

according to the prediction of the Constituent Counting Rules (CCR). As can be 
seen from Eq. 14.41 the scaling law does not depend, apart from a constant factor, on 
the proton scattering angle. 

As seen from the results shown in the previous paragraphs, the CCR prediction 
seems not in agreement, at the photon energies of interest here, with the obtained 
values for the exponent of the scaling law n as shown in Fig. 14.41 

This situation persists even at intermediate scattering angles (around 90° ) where 
momentum transfer to the outgoing proton should be the highest possible and the 
perturbative regime is believed to be better established. 

On the contrary, the present results indicate that pQCD is not plainly applicable 
to the description of the photo-disintegration process in the region of incident photon 
energy explored here. 

In order to better underline this result, the CLAS differential cross section data, 
da/dt, multiplied by the s^^ scaling factor, are shown in Figs l4.8l and l4.9l as a function 
of the incident photon energy in the range 0.55 — 3.0 GeV. 

The s^^ factor has been used to evidence possible deviations of the experimental 
data from the constant behavior expected on the basis of plain pQCD. In addition, 
it should be noticed that the CLAS data cover the incident photon energy range in 
fine steps of 100 MeV thus allowing for an unprecedented precision in sampling the 
cross section behavior in this energy region. 

The panels of Fig. 14.81 and Fig. 14.91 show, from left to right and from top to 
bottom the calculated values for s^^^ in the angular range from = 15° up 
to = 155°. The CLAS data are represented by the solid/red dots, while the 
symbols relative to the data available in the literature are indicated in the legend 
found on the top of Fig. 14. 8[ 

The arrows in each panel indicate the photon energies where Pt = 1 GeV/c. 
According to the CCR, a constant behavior of the s^^ scaled cross section as a 
function of the incident photon energy is expected beyond that limit. As can be seen 
from Fig. 14.81 and Fig. 14.91 this situation is not perfectly realized by the experimental 
data. It is clear that, for forward proton angles, = 15° H-45°, the data scale 
with a power of the total energy s lower then —11 according to the plot shown in 
Fig. 14.41 The data show a trend compatible with the CCR s"^^ scaling roughly in 
the range of scattering angles = 55° 125°. Again, at backward angles, the 
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Figure 4.8: Deuteron photo-disintegration differential cross section expressed as da/dt 
multiplied by s^^ (being t the momentum transfer and s the total energy in the Mandelstam 
notation) in the angular range from 6'™ = 15° up to = 85°. The CLAS results are 
represented by the solid dots. Also shown are previously available data: SLAC [T71[TH] data 
are represented by open squares, JLab Hall C [Ml Hi] data by open triangles, and JLab 
Hall A [26] data by open diamonds. The arrows in each panel, indicate the photon energies 
where the onset of the s~^^ scaling is foreseen. This limit correspond to proton transverse 
momenta of 1 GeV/c. For the jd — >■ pn reaction the proton transverse momentum can 
be expressed = ^S^M^isin^ i&p^)- According to the CCR, a constant behavior of the 
cross section as a function of the incident photon energy is expected when the limit for Pt is 
reached. The solid line represents the QGSM [31] non perturbative calculation for the cross 
section. This model, as discussed in the first Chapter of this thesis, gives a non perturbative 
description of the photo-disintegration cross section assuming that the scattering amplitude 
at very high energy is dominated by the exchange of three valence quarks in the i-channel. 
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Figure 4.9: Deuteron photo-disintegration differential cross section expressed as da/dt 
multiplied by s^^ (being t the momentum transfer and s the total energy in the Mandelstam 
notation) in the angular range from 6*^^ = 95° up to 6^^ = 155°. The CLAS results are 
represented by the solid dots. Also shown are previously available data: SLAG [T71 [TB] data 
are represented by open squares, JLab Hall C [24l [25] data by open triangles, and JLab 
Hall A |26| data by open diamonds. The arrows in each panel, indicate the photon energies 
where the onset of the s~^^ scaling is foreseen. This limit correspond to proton transverse 
momenta of 1 GeV/c. For the — >■ pn reaction the proton transverse momentum can be 
expressed by = ^E-^MdSit? {9^^). According to the CCR, a constant behavior of the 
cross section as a function of the incident photon energy is expected when the limit for Pt is 
reached. It can be seen that the asymptotic behavior for the cross section is well established 
from 



75° up to = 125° while at backward angles, the cross section dependence 



on the photon energy follows powers of s different from —11 as shown in Fig. 14.41 The solid 
line represents the QGSM [21], a non perturbative calculation for the cross section. 



cross section dependence on the photon energy foUows powers of s quite lower than 
-11. 

The observed deviation from the exact n = —11 prediction of plain pQCD for 
Pt ^ 1 GeV found for the exponent of the scaling power n may indicate that pQCD 
is not fully applicable in this energy region. 

On the contrary, a non perturbative approach may be more effective in repro- 
ducing the behavior of the deuteron photo-disintegration cross section. This done 
by the non perturbative QGSM calculation for the cross section, which is shown as 
a solid line in Figs. 14.81 and 14.91 Again, it is seen that a good agreement with the 
data is found for the QGSM predictions at all proton scattering angles. 



19.Q 



4.3 Summary and Conclusions 



The aim of the present research project is the complete measurement of the 
deuteron photo-disintegration differential cross section in the few GeV region of 
incident photon energy. 

This reaction is well suited for studying the interplay of nuclear and particle 
physics in order to find out how and at what energy the transition from the hadronic 
picture of the deuteron, established in the reaction at low energies (below 1 GeV), 
to the quark-gluon picture, takes place. 

With respect to this goal, the physical motivation for such an investigation has 
been discussed, showing that the already published data do not fill the gap of knowl- 
edge on the physical mechanisms underlying the photo-disintegration process. 

Different theoretical models have been developed to describe the deuteron photo- 
disintegration cross section. In particular, the QCD inspired models have been 
reviewed since there is common agreement on the fact that quark and gluons are 
more appropriate degrees of freedom for the description of the process in this energy 
region. The theoretical questioning is not on the use of QCD but, if perturbative 
QCD can be applied at this incident photon energies, or, on the contrary, non 
perturbative QCD must be necessarily used. 

In order to answer the above question, the photo-disintegration cross section 
needs to be known with great statistical accuracy and over a broad angular range 
so to investigate its behavior in the most critical kinematic regions. 

To this aim, the — t- pn cross section has been measured with the CLAS large 
angle spectrometer at the Jefferson Lab Hall B, over a wide range of proton scattering 
angles (15° < < 155°) using a photon beam in the energy range from 0.55 up 
to 3.0 GeV. 

This powerful investigation tool, combining the unique characteristic of a nearly 
47r acceptance in conjunction with a tagged photon beam of high energy and intensity 
allowed the collection of a very high statistics and high quality dataset during the 
end of year 1999. 

The analysis of this large amount of data constitutes the central part of this 
work and has been fully reported in the third Chapter of this thesis to illustrate the 
procedures and techniques used to extract the differential cross section data. 

In the last part of the thesis, the new cross section data are discussed in order 
to contribute to the answering of the questions posed in the opening Chapter. 

To this aim, the differential cross section is shown as a function of the total 
energy s. The new CLAS data add significant information to the previous scenario, 
allowing the study of the cross section decrease as s" with an very good accuracy 
and over a wide angular range. 

In addition, the different exponents n have been calculated over the complete 
range of proton scattering angles and an asymmetry is found between forward 
(n = — 9) and backward angles (n = — 7). 

This result is reflected from the angular distributions of the cross section. The 
data measured at very forward and backward angles show that the photo-disintegration 
cross section exhibits a persistent asymmetry over the whole incident photon energy 
range. 

For what concerns the theoretical interpretations, all cross section data (new 
and published) have been compared to the predictions of the available theoretical 



models based on QCD degrees of freedom. 

The conclusion is a clear indication that pQCD based models do not fully repro- 
duce the data behavior in this region of incident photon energies. On the contrary, 
the non perturbative calculation provided by the Quark Gluon String Model has 
been able to well describe the overall data. 

This result can be considered a further indication that a pQCD interpretation 
of the deuteron photo-disintegration process may not be adequate in the few GeV 
region and a complete non perturbative theory is needed. 
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Table 4.1: Differential cross section data. 
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